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Nonlinear Transfer Functions with Thyrite* 
L. D. KOVACHT anp W. COMLEYT 


ummary—Since the publication of a previous paper! the uses 
hyrite? in the synthesis of nonlinear transfer functions in analog 
puting have been greatly expanded. Moreover, new and con- 
irably improved methods of Thyrite selection and matching have 
n developed. 
The present paper shows how rational exponents can be repre- 
ted. Functions of the form y=kx”", with 1/5<n<6 are readily 
ained. Circuit configurations generating the functions 


| 


| 2x \ 1-74 2 a 1.74 
Pe 2)" 1-2-9) 
| Pid 7 2 


ch provide excellent approximations to the cosine and sine func- 
‘Ss, respectively, are given and an improved quarter-square multi- 
r using Thyrite squaring is discussed. 

‘A study was also made of the physical factors influencing the 
formance and accuracy of Thyrite and quantitative experimental 
a regarding these factors are presented. 


CHARACTERISTICS OF THYRITE 
Be izcar et resistors are characterized by a non- 


linear behavior in that the current varies as some 

power of the applied voltage. This behavior may 
expressed generally by (1), where m varies over a 
asiderable range. 


| 1 = ke. (1) 
e value of 2 is primarily dependent upon the type of 
-yrite and, to a lesser extent, upon the nonuniformity 
ind among samples of a given type. Fortunately, the 
aracteristics of any given sample will be relatively 
nstant when operated conservatively and within the 
ven power rating. 
The most suitable Thyrite types for the generation 
a given exponent are given below. 


Catalog Exponent Maximum 

Number Range Current 
8396839G1 <7 <5 5 ma 
8396832G1 DSSS 10 ma 


The behavior of the Thyrite resistor as given in (1) 
ist generally be modified in order to render the mate- 
1 useful for computing applications. The basic ex- 
nent of the nonlinear resistor must be somewhat 
her than that actually required for a given applica- 
yn. Resistance is connected in series, thereby tending 
linearize the characteristic of the combination. This, 
‘course, reduces the value of 2. The amount of re- 
tance which must be added will be determined by the 
ecise value of required. After proper adjustment, 


* Manuscript received by the PGEC, April 26, 1957; revised 
nuscript received, January 23, 1958. i 
+ Douglas Aircraft Co., Inc., El Segundo, Calif. a5 : 
LL. D. Kovach and W. Comley, “An analog multiplier using 
rite,” IRE TRANS. ON ELECTRONIC Computers, vol. EC-3, pp. 
: 954. 

1 nae is a registered trademark of the Carboloy Dept. of 
» General Electric Co. 


the voltage drop across the series resistance will bear 
the relationship with the input voltage given in (2). 


€o = ke;”. (@Z) 


In order that the exponent remain constant over a 
wide range of input voltages the original value of » for 
the Thyrite should not be significantly higher than the 
final required value. It is fortunate that the value of n 
for a given type of Thyrite resistor will vary as much 
as +20 per cent. This makes it possible to select samples 
particularly suitable for a wide variety of applications. 


Factors AFFECTING OVER-ALL ACCURACY 


An extensive program of tests was conducted to de- 
termine the accuracy obtainable with Thyrite and to 
define those factors which have significant effect. The 
basic accuracy of the squaring operation, 7.e., the con- 
formity of the Thyrite to the true square law was found 
to vary between 0.3 per cent and 1.0 per cent of full- 
scale for most samples. A surprisingly large number 
yielded accuracies better than 0.5 per cent. 

The basic accuracy was determined by means of the 
system illustrated in Fig. 1. The method is based upon 
the comparison of the simulated function with the 
function as generated by the proper number of succes- 
sive integrations of a constant voltage. In this manner a 
perfect parabola or cubic may be generated for the 
purpose of comparison. A typical error curve obtained 
in this manner is shown in Fig. 2. The importance of 
obtaining the optimum values of series resistance for a 
given Thyrite sample is illustrated in Fig. 3 where, for 
several samples, the series resistance has been plotted 
against the resulting basic accuracy. An examination 
of the curves indicates clearly that there is a unique 
value of series resistance which yields the best result. 

It was found that the optimum current range was 
between 0 and 5 ma. Larger currents gave less accurate 
results due principally to internal heating. Overloads 
of 100 per cent will not cause permanent modification 
in the characteristics of the material, but larger, pro- 
longed overloads may result in such changes as to neces- 
sitate a redetermination of the optimum series re- 
sistance value. Some means of overload protection, such 
as amplitude limiting, may be desirable in those cases 
where problem variable peak voltages are not suffi- 
ciently predictable. 

The adverse effects due to temperature have been 
considered separately by means of the circuit of Fig. 4. 
This arrangement provides the means for applying a 
constant input voltage and monitoring the output. The 
importance of restricting the current range is indicated 
in Fig. 5 which indicates the accuracy for various con- 
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Fig. 3—Squaring error vs series resistance. 


tinuous peak current values. At each point of this curve 
sufficient time was allowed for the error to stabilize at 
its maximum value. It was found that the error for 
the 5 ma peak value was about 0.4 per cent. The con- 
tinuous application of full current is, of course, a severe 
condition that will rarely be encountered in operation. 
In actual practice it is doubtful if the error due to in- 
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Fig. 6—Error due to ambient temperature variation. 


ternal heating will ever be as large as 0.2 per cent. 

Reduction of the accuracy due to ambient tempera- 
ture variation is dependent entirely on the environment 
in which the material is used. For this reason it is ad- 
visable to operate the Thyrite at a distance from any 
heat generating equipment or, if possible, in a constant 
temperature environment such as provided by the ovens 
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certain types of analog computers. Fig. 6 shows the 
l-scale error that will result for variations in the am- 
nt temperature up to 10°F. 

ortunately, most up-to-date analog installations are 
conditioned and the room temperature does not vary 
re than one or two degrees. Moreover, Thyrite is 
entially insensitive to humidity variations. 
Finally, it should be pointed out that Thyrite dis- 
rys a slightly semiconductive effect. The difference in 
tput for equal positive or negative inputs may amount 
‘as much as 0.3 per cent of full-scale. Assuming that 
|of the adverse effects are additive, it is still apparent 
t accuracies of 1 per cent of full-scale are not difficult 
(obtain. 
In using Thyrite for multiplying by means of the 
arter-square law two pieces of Thyrite are required. 
nce the error curves of the two pieces may be entirely 
ferent, it might appear that increased accuracy can 
achieved by careful selection. However, it can be 
pwn analytically that regardless of the shape of the 
ror curves the full-scale accuracy of the multiplication 
ll be the sum of the maximum errors of the individual 
eces of Thyrite. 
(It should be pointed out that the results presented in 
‘s section apply to the Thyrite having the lower ex- 
ment. Error evaluation on the type suitable for higher 
ponents yielded results similar to those previously 
tntioned, except that a maximum current of 10 ma is 
rmissible. 


SELECTION AND ADJUSTMENT OF EXPONENT 


[fhe circuit configurations for generating various 
actions are given in Table I (pp. 94-96), but first, 
jtable Thyrite resistors and series linear resist- 
ce must be properly chosen and combined in order to 
pvide the required exponents as indicated in the table. 
fe recommended procedure for this selection and ad- 
stment has been indicated in Fig. 1. This method is 
‘table only for the selection and adjustment of com- 
nations intended to provide exponents of 2 and 3. For 
nintegral powers it will become necessary to replace 
= second integrator in Fig. 1 with a sufficiently ac- 
tate function generator which has been adjusted to 
pvide the exact function required. The difference of 
= two functions is then taken in a summing amplifier 
= output of which may be observed on a suitable re- 
rding device. A repetitive setup and oscilloscope will 
pvide the most rapid and most accurate means for 
iking the necessary adjustments. Adjustment may 
w be made by varying the Thyrite series resistor, at 
2 same time maintaining the output amplitude. The 
sired condition will be indicated as the difference of 
= two functions approaches zero. It should be em- 
asized that the peak amplitude of the derived function 
ould be maintained at a level such that the current 
es not exceed the recommended value for maximum 
suracy. If variation of series resistance between the 
iits of approximately 300 to 1000 ohms does not result 
a suitable minimum difference between the two func- 
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tions, the process should be repeated with other Thyrite 
samples until the required degree of correspondence 
between the two functions is indicated. This difference 
function is the absolute error and may be measured to 
determine whether or not the required degree of ac- 
curacy has been obtained. 


APPLICATIONS 


Table I presents in a convenient form the circuit 
configurations and transfer functions which may be 
obtained by use of the Thyrite and linear resistance 
combinations providing the various exponents as de- 
termined by the methods outlined in the previous sec- 
tion. Use of these combinations in the input, feedback, 
and output of operational amplifiers provides a wide 
range of useful functions. Because of the passive nature 
of the nonlinear resistance, insertion losses will result. 
For a typical squaring operation this loss is approxi- 
mately 20. All or part of this may be compensated for 
by appropriate gains in the associated amplifiers. It is 
interesting to note that no Thyrite resistor combination 
providing an exponent greater than three has been used, 
although units capable of somewhat higher exponents 
are available. The excessively high voltages required to 
place these within their optimum range makes their 
use impractical in computing applications. It will be 
noted that in some instances diodes have been used to 
obtain an absolute value, thereby providing the sym- 
metry required of the even functions. 

It should be pointed out that in circuit configurations 
7, 8, and 9 of the table the Thyrite is in the input circuit 
of the amplifier. Due to the variable input impedance in 
these cases, these circuits should not be supplied from 
potentiometers or other high impedance devices. 

The circuit configurations for the sine and cosine 
functions use a parabola-like function which has been 
translated in a suitable manner in order to simulate the 
desired trigonometric function. The Thyrite resistor 
combination has been adjusted to provide a power of 
1.74. It can be shown that this exponent produces the 
best least squares fit to the trigonometric functions. 
Under these conditions the error resulting from the ap- 
proximation will not exceed 1.35 per cent. The angles 
which may be resolved, however, will be subject to the 
restrictions indicated in the table. 

An improved quarter-square multiplier circuit is 
shown in Fig. 7 (p. 97). Balance potentiometers have 
been added to the squaring circuits in order that the for- 
ward resistance of the diodes may be precisely matched, 
thereby preventing nonsymmetry in the squaring opera- 
tion. This type of multiplier has several attractive ad- 
vantages worthy of consideration: 1) high input im- 
pedance; 2) no requirement for special equipment; 
3) wide-band response; and 4) zero output for zero 
input, a difficult characteristic to achieve with con- 
ventional multiplier schemes. It should be noted that all 
elements comprising the multiplier are either linear or 
passive. As a result accuracy is in no way dependent 
upon active nonlinear devices. 
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tinuous peak current values. At each point of this curve 
sufficient time was allowed for the error to stabilize at 
its maximum value. It was found that the error for 
the 5 ma peak value was about 0.4 per cent. The con- 
tinuous application of full current is, of course, a severe 
condition that will rarely be encountered in operation. 
In actual practice it is doubtful if the error due to in- 
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ternal heating will ever be as large as 0.2 per cent. 
Reduction of the accuracy due to ambient tempera- 
ture variation is dependent entirely on the environment 
in which the material is used. For this reason it is ad- 
visable to operate the Thyrite at a distance from any 
heat generating equipment or, if possible, in a constant 
temperature environment such as provided by the ovens 
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ertain types of analog computers. Fig. 6 shows the 
-scale error that will result for variations in the am- 
int temperature up to 10°F. 
fortunately, most up-to-date analog installations are 
conditioned and the room temperature does not vary 
re than one or two degrees. Moreover, Thyrite is 
entially insensitive to humidity variations. 
inally, it should be pointed out that Thyrite dis- 
Lys a slightly semiconductive effect. The difference in 
tput for equal positive or negative inputs may amount 
‘as much as 0.3 per cent of full-scale. Assuming that 
(of the adverse effects are additive, it is still apparent 
at accuracies of 1 per cent of full-scale are not difficult 
lobtain. 
[In using Thyrite for multiplying by means of the 
arter-square law two pieces of Thyrite are required. 
ice the error curves of the two pieces may be entirely 
ferent, it might appear that increased accuracy can 
_achieved by careful selection. However, it can be 
»wn analytically that regardless of the shape of the 
ror curves the full-scale accuracy of the multiplication 
\l be the sum of the maximum errors of the individual 
eces of Thyrite. 
[t should be pointed out that the results presented in 
‘s section apply to the Thyrite having the lower ex- 
nent. Error evaluation on the type suitable for higher 
ponents yielded results similar to those previously 
tntioned, except that a maximum current of 10 ma is 
rmissible. 


SELECTION AND ADJUSTMENT OF EXPONENT 


ihe circuit configurations for generating various 
fictions are given in Table I (pp. 94-96), but first, 
table Thyrite resistors and series linear resist- 
fe must be properly chosen and combined in order to 
pvide the required exponents as indicated in the table. 
‘e recommended procedure for this selection and ad- 
itment has been indicated in Fig. 1. This method is 
‘table only for the selection and adjustment of com- 
sations intended to provide exponents of 2 and 3. For 
nintegral powers it will become necessary to replace 
> second integrator in Fig. 1 with a sufficiently ac- 
rate function generator which has been adjusted to 
pvide the exact function required. The difference of 
> two functions is then taken in a summing amplifier 
> output of which may be observed on a suitable re- 
‘ding device. A repetitive setup and oscilloscope will 
yvide the most rapid and most accurate means for 
king the necessary adjustments. Adjustment may 
w be made by varying the Thyrite series resistor, at 
» same time maintaining the output amplitude. The 
sired condition will be indicated as the difference of 
> two functions approaches zero. It should be em- 
asized that the peak amplitude of the derived function 
5uld be maintained at a level such that the current 
es not exceed the recommended value for maximum 
suracy. If variation of series resistance between the 
1its of approximately 300 to 1000 ohms does not result 
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tions, the process should be repeated with other Thyrite 
samples until the required degree of correspondence 
between the two functions is indicated. This difference 
function is the absolute error and may be measured to 
determine whether or not the required degree of ac- 
curacy has been obtained. 


APPLICATIONS 


Table I presents in a convenient form the circuit 
configurations and transfer functions which may be 
obtained by use of the Thyrite and linear resistance 
combinations providing the various exponents as de- 
termined by the methods outlined in the previous sec- 
tion. Use of these combinations in the input, feedback, 
and output of operational amplifiers provides a wide 
range of useful functions. Because of the passive nature 
of the nonlinear resistance, insertion losses will result. 
For a typical squaring operation this loss is approxi- 
mately 20. All or part of this may be compensated for 
by appropriate gains in the associated amplifiers. It is 
interesting to note that no Thyrite resistor combination 
providing an exponent greater than three has been used, 
although units capable of somewhat higher exponents 
are available. The excessively high voltages required to 
place these within their optimum range makes their 
use impractical in computing applications. It will be 
noted that in some instances diodes have been used to 
obtain an absolute value, thereby providing the sym- 
metry required of the even functions. 

It should be pointed out that in circuit configurations 
7, 8, and 9 of the table the Thyrite is in the input circuit 
of the amplifier. Due to the variable input impedance in 
these cases, these circuits should not be supplied from 
potentiometers or other high impedance devices. 

The circuit configurations for the sine and cosine 
functions use a parabola-like function which has been 
translated in a suitable manner in order to simulate the 
desired trigonometric function. The Thyrite resistor 
combination has been adjusted to provide a power of 
1.74. It can be shown that this exponent produces the 
best least squares fit to the trigonometric functions. 
Under these conditions the error resulting from the ap- 
proximation will not exceed 1.35 per cent. The angles 
which may be resolved, however, will be subject to the 
restrictions indicated in the table. 

An improved quarter-square multiplier circuit is 
shown in Fig. 7 (p. 97). Balance potentiometers have 
been added to the squaring circuits in order that the for- 
ward resistance of the diodes may be precisely matched, 
thereby preventing nonsymmetry in the squaring opera- 
tion. This type of multiplier has several attractive ad- 
vantages worthy of consideration: 1) high input im- 
pedance; 2) no requirement for special equipment; 
3) wide-band response; and 4) zero output for zero 
input, a difficult characteristic to achieve with con- 
ventional multiplier schemes. It should be noted that all 
elements comprising the multiplier are either linear or 
passive. As a result accuracy is in no way dependent 
upon active nonlinear devices. 
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TABLE I (Cond) 


CIRCUIT CONFIGURATION TRANSFER FUNCTION 
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'* 50 ~~ pots to balance forward resistance of diodes. 

~~ input pot to balance outputs of thyrite squaring circuits. 
[L] 10M ~ pot to adjust K. 

|All fixed resistors + 0.25 per cent. 

‘Amplifier drift should not exceed 0.005 volt referred to the input. 
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Fig. 7—Thyrite multiplier circuit. 
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| Summary—Polynomial equations occur quite often in physical 
i mathematical problems and require an exact solution. Several 
fblems in control engineering, aerodynamics, and many other 
ads require the formulation and exact solution of algebraic equa- 
ins of quite a high order. In the past, several mechanical and elec- 
eal machines for the determination of the roots of polynomials have 
tn developed, but all of them suffer from the serious disadvantage 
[t they are very expensive and are incapable of high accuracy and 
lidity of operation. 

‘This paper describes the development, design, and construction 
an inexpensive and portable isograph which is capable of locating 
roots of a polynomial with fairly high accuracy. This isograph can 
ve as a valuable aid to large-scale computing machinery, for which 
roots of a polynomial must be isolated before further calculations 
be performed. It can also serve as a very useful adjunct to other 
thematical processes such as complementary solution of differ- 
ial equations, determination of latent roots of a matrix, and similar 
er applications. 

‘It is hoped that this inexpensive and portable instrument will 
ve extremely useful in the design of servosystems and in other 
ds where the problem may be presented as a polynomial and that 
oretical studies of engineering problems will be greatly aug- 
nted by its use. 


: INTRODUCTION 


-T is well-known that polynomial equations occur 
quite often in physical and mathematical problems 
- and require an exact solution. Several problems in 


* Manuscript received by the PGEC, November 7, 1957; revised 
nuscript received, January 31, 1958. ; : 
+ Dept. of Power Eng., Indian Inst. of Sci., Banalore, India. 


There is no intention to imply that the configurations 
given in the table are inflexible or all inclusive. They are 
intended as a guide to what may be done and other con- 
figurations are possible. These may be added to the 
table as dictated by future requirements. 


CONCLUSION 


Thyrite resistors, properly chosen and adjusted by the 
addition of the proper linear resistance, provide the 
means for generating a large class of nonlinear functions 
at small expense and with very little complication. The 
technique is valuable in releasing more complicated and 
expensive function generating equipment for the more 
difficult tasks for which they were intended, thereby in- 
creasing the total function generating capabilities of any 
analog computing installation. 


A Novel Type of Isograph (Algebraic 


Equation Solver)* 
P. VENKATA RAOT 


servomechanisms, communication engineering, aero- 
dynamics, and many other fields require the formulation 
and exact solution of algebraic equations of quite a high 
order. In the field of applied mechanics, all vibration 
problems involve the solution of a frequency equation 
which turns out to be a polynomial of degree corre- 
sponding to the degree of freedom enjoyed by the sys- 
tem. General analytical methods for obtaining the roots 
of equations as high as the fourth are available in stand- 
ard textbooks on applied mathematics. However, Abel 
has proven that general solutions for polynomials of 
degree higher than the fourth cannot be obtained, and 
that particular solutions are the only ones possible in 
equations of the fifth degree and higher. Several numer- 
ical methods have been developed for the solution of 
higher degree algebraic equations, but each equation 
has to be treated as an individual case and its particular 
solution obtained. Some of the more well-known numer- 
ical methods are Ferrari's Method, Shih-Nge-Lin’s 
Method, and Graffe’s Method. Although these methods 
are quite useful for lower degree algebraic equations, 
they become extremely tedious and time consuming for 
equations of very high order, particularly if complex 
roots are involved. Consequently, several mechanical 
and electrical machines for the determination of the 
roots of polynomials have been developed, but all of 
them suffer from the disadvantage that they are very 
expensive and are incapable of high accuracy and rapid- 
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ity of operation. Some of the more well-known machines 
are the Bell Telephone Laboratories Isograph, the 
Stibitz Isograph, the University of Texas Mechanical 
Synthesizer, Schooley’s Isograph, and the Magslip 
Isograph of Parker and Williams. It is needless to point 
out that large and even medium present day digital 
computers are capable of solving very high-order alge- 
braic equations in a very short time, with extraordinary 
accuracy. 


HARMONIC SYNTHESIS 


The design of most of the mechanical and electrical 
isographs, including the one described in this paper, is 
based on the principle of harmonic synthesis. The 
method of operation can best be explained by first trans- 
forming the general algebraic equation into the usual 
form in which it is generally set on the isograph. 

A typical polynomial of high degree is: 


f(x) = do + aye + aox? + agx? + - - 


where the coefficients a, are real numbers, whose value 
is less than unity. In general, the roots of the equation 
are of the form: 


ane = 05 (1) 


r=u+tjv = re? = r(cos 6+ 7 sin #) 


and there will be as many roots as the order of the equa- 
tion. 
By De Moivre’s theorem 7 = (u?+ 7)!/? and tan 0=v/u. 
Rewriting (1): 
J («) = ao+air(cos +7 sin 6) + aer?(cos 20+7 sin 26) 
+a3r?(cos 30-+7 sin 30)+ - - 


Separating, the real and imaginary parts, and equating 
each to zero 


: dnr™(cos nO+7 sin n0=0. 


- dar” cos nf = 0 (2) 


- dar” sinn@ = 0. (3) 


ao + ayr cos 6 + aor? cos 26: - 
ay sin 8 + aor? sin 26 + agr? sin 30 - - 


Multiply (2) by sin wt and (3) by cos wt, where 
w=2nf, f being an arbitrarily chosen independent vari- 
able. After adding these two products and collecting like 
terms in powers of 7, one obtains 


ao Sin wt + ayr sin (wt + 6) + aor? sin (wt + 26) 
+ azr? sin (wt + 30) +++ aar” sin (wi + 19) = 0. (4) 


Each of the terms of the above equation could be repre- 
sented by a sine wave voltage, the modulus and phase 
angle of which can be varied over the desired range. 


Basic PRINCIPLE 


The isographs developed by Schooley, Hart and 
Travis, and others consist of m+1, identical single-phase 
ac alternators of which m are provided with rotatable 
stators. All the alternators are driven by a motor on a 
common shaft in order to ensure that they operate at 
the same frequency. The rotatable stators are connected 
to a common shaft through gear ratios of 1:1, 2:1, 3:1, 

- , +1, so that as the hand crank is rotated slowly, 
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the stator of the first alternator turns through 0, the 
second through 26, etc., down the line, to the mth stator 
which would turn through 70. q 

The output voltages of the various alternators would 
have the phase relations indicated in Fig. 1, and they 
are impressed across coefficient potentiometers set to 
the fractional values a1, @2, a3, * - + , dn. The outputs of 
these coefficient potentiometers are then impressed 
across suitably ganged modulus potentiometers with _ 
tapped windings, in order to produce the various powers 
of the modulus. The correct phase relationship of the 
alternators, with reference to the first one, is thus ad- 
justed by means of the knob on the hand crank whose 
calibration extends from zero to 360°. 

In the actual operation of this type of isograph, an 
estimated value of the modulus is set on the ganged 
modulus potentiometers and the crank handle knob is 
slowly turned while observing the indicating instrument 
for a null or a minimum deflection. The outputs of the 
various ganged potentiometers are connected in series 
together with the voltage a) and the resulting voltage is 
read by the meter. By trial and error, a satisfactory null 
point can be obtained and the setting of the modulus 
potentiometer and the position of the knob of the hand 
crank would determine the modulus and argument of 
the root, respectively. This procedure should, of course, 
be repeated until all the roots of the polynomial have 
been determined. With this type of isograph it is also 
possible to solve equations with complex coefficients by 
setting the stator of each alternator to an initial phase 
displacement corresponding to the phase angle of the 
complex coefficient instead of at zero or 180° as is the 
case with ordinary positive and negative coefficients, 
respectively. 

In the isograph developed and built by Marshall the 
different harmonic components are generated by using 
square generators of the rotating commutator type in 
conjunction with cascaded ganged potentiometers and 
unity gain isolating amplifiers. Filtering circuits shape 
the waves into sinusoidal form and potentiometers 
serve to bring in the coefficients a, of the equation under 


_ study. The outputs are split electrically into the re- 


spective sine and cosine components. summed very 
accurately, and they are then displayed on a cathode- 
ray oscilloscope much as the mechanical isograph does 
it on a plain drawing board. 

Choudhury’s isograph, also designed on the principle 
of harmonic synthesis, is a purely electronic equipment 
without any mechanical moving parts. The different 
harmonic components are generated with the help of a 
delay line, fed from a matched generator, the frequency — 
of which is swept by the time base voltage of the cathode 
ray oscillograph. The variation in frequency of the os- 
cillator corresponds to the variation of the argument 
of the harmonic components. It is thus quite evident 
that, by controlling the amount of frequency sweep of 
the oscillator, any desired interval of the argument can 
be suitably expanded, thereby increasing the accuracy 
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measurement. An interesting method of solving 
wher degree algebraic equations was suggested re- 
ntly by Barker, who employed a transformer analog 
mmputer of the Blackburn type. Bubb suggested the 
e of operational amplifiers for generating polynomials 
id determining their zeros. 
In the Magslip isograph developed and constructed 
- Parker and Williams the terms of the polynomials 
= simulated by ac voltages of varying phase obtained 
om Magslip transmitters excited by a balanced three- 
ase supply. The application of the principle of har- 
onic synthesis resulted in a fairly compact machine, 
pable of finding the roots with high accuracy. This 
soraph, like the earlier ones, could also be used for the 
lution of polynomials with complex coefficients and 
- their evaluation at any point in the complex domain. 
1e Magslip isograph designed to solve sixth degree 
yebraic equations with real or complex coefficients 
cupies a standard rack 19 inches wide and 6 feet high, 
d the three-phase supply to the equipment is gen- 
ated by an oscillator the components of which are 
used in a similar rack of half the height. 


CUMULATIVE PHASE SHIFTING 


The design of the isograph described in this paper is 
sed on the principle of cumulative phase shifting for 
e generation of the various terms of the polynomials 
1ose roots are to be determined. The application of 


Solver) 99 


Fig. 1—Schematic arrangement of the isograph built by Schooley and others. 


this simple principle results in a compact, inexpensive, 
and portable instrument, the accuracy of which com- 
pares very favorably with that obtained by the very 
best of isographs built so far. As constructed, the iso- 
graph is capable of solving sixth degree equations having 
real coefficients, but the principle could be extended to 
deal with equations of higher degrees. It occupies a 
rack 19 inches wide and 16 inches high and it could be 
operated on the mains supply. 

The block diagram of Fig. 2 illustrates the basic 
principle of working the isograph developed and built in 
the Indian Institute of Science. From the block dia- 
gram it is easily seen that the isograph essentially con- 
sists of a single type of unit which gives an output equal 
to rv times the input but shifted in time phase by a con- 
trollable angle 6. The modulus of the root is set by 
ganged decade switches selecting the correct tap ratios 
of highly accurate three-winding transformers of special 
design, and the argument of the root is set by means of 
suitably ganged potentiometers in a set of simple RC 
phase shifter circuits. Referring now to (4) and Fig. 2, it 
can be seen that the isograph generates each one of the 
terms of the polynomial under study. Terms of the type 
rei? are obtained, simply by cascading the requisite 
number of units of identical design and providing con- 
trol for 7 and @ with the aid of ganged switches. Of 
course, it is necessary to multiply these terms by the 
fractional values of the coefficients before the summa- 
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Fig. 2—Block diagram of the isograph described in this paper. 


Tens 


Units 


Fig. 3—Circuit arrangement of a typical complex multiplier unit. 


tion of all the terms of the polynomial is effected by - 


connecting the various voltages in series with a suitable 
measuring instrument. The mode of operation is thus 
to set the coefficients of the equation in the manner de- 
scribed below and to vary the modulus and phase angle 
of the variable until a zero reading is obtained. It is 
quite obvious that the isograph could also be used for 
the evaluation of a polynomial at any point in the com- 
plex domain. 


THE re”? UNIT 


Fig. 3 gives the actual circuit arrangement of a 
typical re” unit. It essentially consists of a simple RC 
phase shifting network and a three-winding, multi- 
tapped transformer coupled by a cathode follower which 
acts merely as an impedance matching device. The 
equipment provides an output voltage which is exactly 


IRE TRANSACTIONS GN ELECTRONIC COMPUTERS 


equal to r times the magnitude of the input voltage and 
at an angle 0 degrees relative to the voltage applied to 
the input terminals. The modulus ¢ is set as the turns 
ratio on the r winding of the transformer and @ is set 
as the indication of a calibrated potentiometer in a 
simple RC phase shifting network, consisting of two 
exactly identical capacitors in series with two ganged 
resistors connected in the manner shown, in Fig. 4(a). 
The vector diagram, in Fig. 4(b), shows the various 
voltages in the network; it can be seen that the output 
voltage is always equal to the input voltage but is dis- 
placed from the latter by an angle @ which can be 
changed from 0° when R is equal to zero to nearly 180° 
when R is equal to ~. 

The output voltage, therefore, can be set at any 
phase angle with reference to the input voltage merely 
by controlling the value of R with the aid of a ganged 
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(b) 


Fig. 4—(a) Circuit arrangement of the phase shifter. 
(b) Vector diagram—phase shifter. 


[tentiometer. This voltage is then applied to the input 
rminals of a cathode-follower circuit which presents 
very high input impedance so that the performance of 
= phase shifter is not adversely affected. The output 
ithe cathode follower is impressed across the primary 
‘a three-winding transformer, built out of mumetal 
impings with a fairly large area of cross section for 
» core. The primary winding is provided with a num- 
+ of taps which can be adjusted over a small range 
id the two secondaries are tapped in decade arrange- 
tnt so as to secure tap ratios from 0.001 to 0.999. The 
mmber of turns of the primary winding is adjusted so 
at the voltage at a secondary of 2000 turns (which is 
mj included for checking purposes) and the input 
itage to the phase shifter are exactly equal. These two 
itage signals are brought to a test panel at the back- 
ve of the isograph in order to check frequently for the 
ity gain setting of the primary of the transformer. As 
well known the gain of a cathode follower is about 0.9 
H minor variations of this, if any, can be compensated 
by varying the tap on the primary side of the output 
‘nsformer. This, evidently, does not affect the settings 
vy and a. To maintain the accuracy of this unit over a 
Je range of working voltages, regulation of the trans- 
mer is reduced to an absolute minimum by employing 
ecially treated high permeability stampings and 
ting the base number of turns as 2000. As the maxi- 
im load on the output transformer corresponds to the 
nacitive reactance of the two capacitors in parallel 
nen R is set at zero), their value is chosen to be as high 
possible consistent with the availability of suitable 
ues of R to sweep 6 through the range from 0 to 180°. 
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The above arrangement is, in effect, a highly stabi- 
lized amplifier of “unity gain” from the input of the 
phase shifter to the output at the 2000-turn secondary 
winding. Other values for the gain setting of the ampli- 
fier from 0.001 to 0.999 are secured independently from 
the two secondaries by means of decade switches. The 
appropriate fractional value of the coefficient is set on 
one of the secondaries. The setting on the other second- 
ary corresponds to the modulus and it is manipulated by 
suitably ganged decade switches. The unit illustrated 
above is thus capable of producing two output voltages, 
re? and are. The output of the 7 winding is applied 
to a similar arrangement consisting of a phase-shifter 
circuit and a three-winding transformer connected by a 
cathode follower. Fig. 5 isa photograph of the isograph 
designed to solve an algebraic equation of the sixth 
order, with real coefficients. 


Fig. 5—Photograph of the isograph described in this paper. 


Root FINDING WITH THE [SOGRAPH 


As in the case of all the other isographs described in 
technical literature, practical design considerations 
have made it necessary to limit the range of the modulus 
to values less than unity. The range of variation of the 
argument, however, is limited to 0-180°, and by provid- 
ing reversing switches on the primaries of the trans- 
formers it is possible to sweep @ through 0-360°. In 
order to bring all the coefficients within the compass of 
the isograph dials, it is usually necessary to divide all 
the terms of (1) by the largest coefficient or by some 
other convenient number, so that all the coefficients as- 
sume fractional values of comparable magnitude. The 
procedure for determining the roots whose modulus is 
less than unity is quite simple and straightforward. To 
take into account all those roots whose modulus is 
greater than unity, it is very necessary first to transpose 
the equation so that the reciproval of r appears in the 
transformed equation. We make the transformation 


X= ely 


By making the simple substitution and multiplying 
all the terms of the polynomial by X”, one gets 


Zp 0; aoX* + aX" 0 
Ceaaret  Taes ) a9 + a4 Ded UF 3 ANG 
This equation can be written down by inspection merely 
by inverting the order of the coefficients and introduc- 
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ing a new independent variable X which is the reciprocal 
of the actual root. This modified equation can now be 
examined for roots in the usual manner and the range 
0<|x|<1 for the modified equation corresponds to the 
range 1<|r| < in the original equation. 

In the case of all the polynomials where the ratio of 
the largest to the smallest coefficients is extraordinarily 
high, a simple substitution of the type z=kx may be 
made with advantage. 


a2? 638° i Re 


Peg aL 
The scaling factor Rk is selected in such a way as to re- 
duce the very wide disparity in magnitude of the coef- 
ficients of the given equation. The number of complex 
roots will be equal to the difference between the order 
of the equation and the number of real roots. Further- 
more, if the coefficients are real, as they are in the large 
majority of problems in servomechanisms and com- 
munication engineering, the complex roots must neces- 
sarily appear in conjugate pairs, so that only the upper 
half of the complex plane needs scanning. It is largely 
for this reason that provision is made in this isograph 
for a variation of 0-180° only, for the argument. 

The actual method of operation of the isograph is thus 
to set up the coefficients of the equation on the panel 
and then to vary the modulus and argument in a judi- 
cious manner until a zero reading is obtained on the 
indicating instrument. In actual use, 0 is set at zero and 
a search is made for positive real roots, r being varied 
until a null reading is obtained on the display meter. 
Negative real roots are obtained in an exactly similar 
fashion by setting @ at 180°. For the determination of 
complex roots, it is found in practice that it would be 
most convenient to vary one control, usually 0, until a 
minimum reading is obtained and then to vary the other 
control to reduce the minimum to zero if possible. In this 
way a root is approached fairly rapidly through suc- 
cessive approximations of the argument and modulus. 
The indicating instrument would give a sharply rising 
or falling indication with each change in either of the 


f(x) = f(2) = ao + 2 a 0. 


two controls signifying the approach to or retreat from 


a zone containing a root. 


MEASURING INSTRUMENT 


When the settings of the modulus and argument 
switches correspond to a root, the fundamental com- 
ponents of the voltage balance out; however, their har- 
monic components give a residual voltage, as the phase 
shift, introduced for different harmonics, would be dif- 
ferent. In order to be able to measure the signal strength 
to the high accuracy desired, it is absolutely necessary 
to filter out all traces of harmonic distortion from the 
output signal. This can be achieved by incorporating a 
passive discriminating filter which effectively attenuates 
the harmonics and passes the fundamental without any 
appreciable change in its magnitude. A more satisfac- 
tory discrimination, however, can be secured by a fre- 
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quency selective amplifier which gives an output with a 
very high amplification for the desired frequency and 
little or no amplification for the harmonic components. 

The output of this amplifier is fed back to the input, 
through a twin-T filter. This filter is of the rejector 
type and passes all the harmonics without much atten- 
uation except the fundamental. The negative feedback is 
operative very largely for the harmonics, while it is in- 
effective for the fundamental. In the frequency selective 
amplifier built here, for use in conjunction with a tube 
voltmeter, the discrimination is of the order of 100:1 so 
that the relative amplification of the fundamental is © 
about 100. 


COMPLEX COEFFICIENTS 


An important extension that has interesting possi- 
bilities is the solution of algebraic equations involving — 
complex coefficients. It is obviously necessary to con- — 
nect an additional phase shifting network to the output 
side of the coefficient winding of the transformer of each 
of the ve” units. A very practical use for complex coef- 
ficients is in the transformation of the point of origin — 
in the complex plane to some other convenient point, 
approximately located in the neighborhood of a cluster 
of complex roots. By treating this new point as the ~ 
origin, in exploring for the various roots located as de- 
scribed, their determination is greatly simplified. 


CONCLUSIONS 


The isograph described in this paper is capable of — 
locating the roots of a polynomial with high accuracy 
and in a short time of about one hour, provided that the 
roots are fairly well separated. In equations having 
closely spaced roots the minima obtained on the meter 
are rather broad and it is comparatively difficult to 
determine the roots to a high degree of accuracy, as it is 
even in the case of the isographs developed earlier. 
With the aid of a small amount of computational work | 
on a desk calculating machine, it is possible to find all 
the roots of an algebraic equation correct to the fourth 
decimal place. With the isograph itself a precision of | 
about 0.8 per cent can easily be achieved. The isograph © 
is capable of yielding solutions with an accuracy suffi- 
ciently great for most engineering problems and with a 
rapidity of operation far exceeding that of any other 
device developed so far. 

The isograph also lends itself to several other applica- 
tions as well. It can serve as a valuable aid to large- 
scale computing machinery, for which the roots of a 
polynomial must be isolated before further calculations 
can be performed. It can also serve as a very useful — 
adjunct to other mathematical processes, such as com- 
plementary solution of differential equations, determina- 
tion of latent roots of a matrix, and similar applications. — 

It is hoped that this inexpensive and portable instru- 
ment will prove extremely useful in the design of servo-_ 
mechanisms and in other fields where the problem — 
may be presented as a polynomial and that theoretical | 
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studies of engineering problems will be greatly aug- 
mented by its use. It will also be of great value in the 
calculation of the transient behavior of linear electrical 
and mechanical systems. 
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Logically Micro-Programmed Computers” 
JOHN V. BLANKENBAKERft 


_ Summary—Extremely simple digital computers exploiting the 
-oncepts of simulation and micro-programming are described. Logi- 
mal rather than arithmetic micro-programming operations are em- 
sloyed for generality and greater simplicity. The resulting class of 
-omputers can duplicate the behavior of any finite state digital com- 
outer except in solution time. Also, design techniques are given for 
somputers employing only multiple-bit time delays. 


INTRODUCTION 


ponents by virtue of its structure or form than an- 

other functionally equivalent computer requires, 
che first computer is potentially preferable on the basis 
of cost and reliability. In designing for an application 
with a fixed problem, it would be desirable to keep the 
somputer as simple as is consistent with the required 
solution speed. In many problems this is very low. 
) Unfortunately, it is not always clear how to achieve 
he simpler forms. In especially simple applications 
nore computing capacity may be postulated than is 
1ecessary and it is then difficult to see how the struc- 
ure may be reduced. Experience in the computer arts 
1as shown that it is generally easier to add additional 


‘Ee one digital computer design requires less com- 
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structure to a design to increase its computational abil- 
ity than it is to accomplish the converse of simplifying 
the equipment with a reduction of its computing capac- 
ity. 

One simple form of digital computer, which might be 
used as the basis for more complex designs, is described. 
In its simplest form the area of application would be 
small. With additional structure, however, the range of 
problems amenable to solution can be extended. While 
some brief suggestions are given for means of increasing 
the computational ability of the machine, these thoughts 
are not pursued in detail. 


THE Basic COMPUTER 


The first computer is described by (1). The logical 
notation employed is “+” for the nonexclusive OR, 
“.” for AND (sometimes omitted), and “~” for the 
complement of a variable or function.! Subscripts 
identify time intervals. Synchronous techniques are 
assumed in the serial computer, Fig. 1. Z is a one-bit 
time delay, R is an n-bit time delay capable of storing x 
bits of information, and similarly M and WN are each mn 


1E. C. Nelson, “An algebraic theory for use in digital computer 
design,” IRE TRANS. ON ELECTRONIC CoMPUTERS, vol. EC-3, pp. 12— 
21, September, 1954. 
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Fig. 1—Basic computer. 


bit time delays. L;, R:, M., N; are the “true” outputs of 
their respective devices, while W; is the input to R. At 
any one time M, and JN; describe one of the four micro- 
instructions which the machine can execute. The M-N 
sequences of mn “0O”s and “1”s do not change for a 
given set of computer characteristics. 


Liga = Li: My, + Ly-Re+ MiMi (1a) 
Ry => W; = M,-N,-Input + Ni R: +. M,Ni-Lt (1b) 


M, = Veen ic) 
N: a LN a (1d) 
Output — Ri, when WIN: =e (1e) 


The four micro-instructions are given in Table I. Asa 
result of the do nothing instruction, L retains its present 
state and the input to R is equal to the output bit of R. 

TABLE I 
INSTRUCTION CODE FOR Basic COMPUTER 


Code 
Name Symbol Function 
Do nothing — 0 O0O Wr=Ri, Li =Lt 
Input-Output “é 0 1 W,=input, output=R;, Lia=L; 
Logical product O 1 O Wr=Ry, Lip =H=Li Ki 
Store W 1 1 Wi=Li, lope SO 


The input bit from external sources is entered into R 
and one bit of the computer output is taken from R; in 
the input-output instruction. The logical product of L, 
and R; is formed and entered into ZL in logical product. 
In store a bit equal to ZL; is entered into R and Ly41 is 
determined as a “1”. 


PROGRAMMING 


Programming usually consists of two steps. For a 
given problem one approach might be to program the 
LMC (Logically Micro-programmed Computer) as a 
general purpose “arithmetic” computer. In general in 
this first step, an M-N program is found which causes 
the LMC to simulate another digital machine. This M-N 
program of micro-instructions constitutes a description 
of a number of Boolean algebraic equations describing 
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the simulated machine which the LMC evaluates. If the 
LMC is to have permanently the same M-N program, 
this may be entered into the machine when it is built; 
otherwise an M-N program loading mechanism is re- 
quired. 

As the second step the simulated machine is pro- 
grammed, presumably by use of a more or less conven- 
tional instruction code. The state of the simulated 
computer as given by its memory elements is repre- 
sented in the LMC by the information in R. Program- 
ming this simulated computer consists of finding the 
initial state of its memory bits. Since there are adequate 
references on this second part of programming, no dis- 
cussion will be given concerning it. ; 

Programming will be described by way of a few spe- 
cialized examples. In the way of notation the following 
is introduced. Then bit positions of R, which circulate, 
are designated as P1, P2, - - - , Pn. Time is divided into 
major and minor cycles corresponding respectively to 
the circulation times of the M-N channels and of R or 
the variable loop. Each bit time of a major cycle may 
be denoted by two numbers. From a time, which is 
taken to be coincident with the appearance of P1 at 
R;, the bit times are numbered from 1 to ” correspond- 
ing to the variable position number. This will be the 
second half of the time number. The first half of the 
time number denotes the minor cycle and runs from 1 
to m. Thus, ¢(6,4) denotes the fourth bit time of the 
sixth minor cycle. Any micro-instruction in the M-N 
program may be taken as the fiducial mark (if the P’s 
are yet unnumbered). 

The micro-instruction program can conveniently be 
written in the form of an ” by m array. The 7 columns 
correspond to the bit times in a minor cycle (and to the 
variable position appearing at R;), while rows of the array 
correspond to minor cycles. Entries are the micro-in- 
structions. 

The principal basic programming techniques are il- 
lustrated in the micro-instruction array of Table II. The 
micro-instruction occurring at (1,1) (see Table I for 
instruction symbols) has no part in the programming. 
being illustrated. As a result, however, Z is in the “1” 
state at #(1,2). At this time the logical product instruc- 
tion determines Li, as Ly-R;=1-R;; or the comple- 
ment of the variable in P2 has been entered into ZL. 


AWANIBILIS, 1h 


PROGRAM OF MiIcRO-INSTRUCTION FOR LoGicaL COMPUTATION 
of f(V2, V3, V5, V6)=V2+V3- V5: V6 


n 


OWN 
| | 
o| x39 


* The nine m numbers are the R bit positions and also define the 
variables at time ¢(1,1). 
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his value holds through a series of do nothing instruc- 
ions until ¢(2,2), when it is stored back into the P2 
osition of R. If originally the Pz position were said to 
store a variable named V2, the variable now in the P2 
peer must be called V2. A careful distinction must 
e made between the variable position and the variable 
herein. 

The two bit time computation starting at ¢(2,3) is 
imilar to the previous one but differs in. that V3 is 
stored in the P4 position. This merely illustrates that 
two variables may be stored in R which are comple- 
mentary. A slightly more complex computation of the 
AND function, V3-V5- V6, begins at ¢(3,4). After the 
(3,4) instruction, L contains V3. The logical product 
instruction at ¢(3,5) causes L;-R;, equal to V3- V5, to 
be entered into L. At ¢(3,6) the logical product of 
Lis V3-V5 and of R,=V6 is formed and entered into 
L. Thus at ¢(3, 7) the value of the quantity stored into 
R is equal to V3-V5-V6. The ability to complement 
variables and evaluate AND functions has now been 
demonstrated. 

The next two series illustrate the evaluation of OR 
Functions from the previous two operations. The de- 
sired objective is to have a variable in P9 equal to 
V2+V3-V5-V6. By Boolean manipulations 


V2 + V3-V5-V6 = (V2)-(V3-V5-V6), (2) 


which removes all “+” signs by using AND and com- 
plementation operations. P7 at the present point in the 
computation contains V3-V5- V6. The series which be- 
gins in the program matrix at ¢(3,8) and ends with ¢(4, 
8) evaluates the right side of (2), except for the upper- 
most complementation sign. This final step is completed 
with the series ending at ¢(5,9) so that the P9 variable 
is equal to the left-hand side of (2). 
- The elements of encoding, or the translation of Boo- 
lean algebraic equations into micro-instructions, have 
been illustrated by the previous examples. The set of 
equations which are to be encoded and which define the 
computer to be simulated may be obtained from an exist- 
ing computer design? or they may be obtained by orig- 
inal work.?4 

Any digital computer can be cast in the form of a set 
of logical difference equations, 


X ovi® = fi(X,}, xy, XS, a a) AP), t= Ls 2, Saal oe (3) 


where the X‘ are one-bit delay elements as denoted by 
the subscripts s+1 and s [compare with (1a)]. This 
time difference (which can be of varying periods) is the 
reason for calling these equations difference equations. 
In distinction to the use of ¢ for real time, s is used for 


2S, P. Frankel, “The logical design of a simple general purpose 
computer,” IRE Trans. ON ELECTRONIC CoMPUTERS, vol. EC-6, pp. 
5-14; March, 1957. : 

3 D. A. Huffman, “The synthesis of sequential switching circuits,” 
J. Franklin Inst., vol. 257, p. 161, March, and p. 275, April, 1954. 

4M. Phister, Jr., “Logical Design of Digital Computers,” John 
Wiley and Sons, Inc., New York, N. Y.; 1958. 
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simulated time. For the present purposes, the functions 
f' represent combinations of the logical operators AND, 
OR, and complementation. Computer memory elements 
which do not seem to represent one-bit delays, for exam- 
ple R or the n-bit delay of the LMC computer, can be 
represented as a series of one-bit delays. External inputs 
can be treated as outputs of delays which have unknown 
inputs and the outputs can be represented as inputs to 
delays whose outputs are never used. 

Once a set of equations of the form of (3) are obtained, 
they would be recast to eliminate the OR operation. 
This will not increase the number of difference equa- 
tions. Taking (1a), merely by way of example, if the OR 
operation is eliminated, it would appear as 


Lot re (L,-M,)- (Le: Rs): (M,N). (4) 


This last equation is too complex to be evaluated di- 
rectly. As implied in the example of Table II, additional 
variables are defined. In this case, they could be 


A= ln Ml By = oR, C, = tN (Sa) 
and also 
D, = A,-B,°C, (Sb) 
Then, 
bah a Ihe (6) 


To preserve the time relationship of (4), an s sub- 
script is used with the intermediate variables A, B, C, 
and D. A physical interpretation is that no delay is asso- 
ciated with formation of these variables. In micropro- 
gramming any single equation, no distinction is made 
between delay and nondelay equations. What is im- 
portant is the order in which the equations are encoded 
and evaluated. Following the time when all variables 
may be said to have their s value, all nondelay variables 
must be computed before using the nondelay variables 
in computing a delay variable. Also, a series of non- 
delay equations must be evaluated in a sequence run- 
ning from the s variables to the s+1 variables. 

Delays could also have been assigned to (5a) and 
then (5b) and (6) would be nondelay equations for 
time s+1. Also, fractional delays could be assigned to 
each equation. 

Another problem arises because all variables, in par- 
ticular the delay variables, are not advanced simul- 
taneously to their new or s+1 value. If the computation 
is at a point where all delay variables have their s value, 
these values cannot, in general, be destroyed or altered 
until all the s+1 values have been computed. Actually, 
this problem arises only if synchronous machines, 
which in effect have equal delays throughout, are being 
modeled. The design of asynchronous machines as- 
sumes that the delay times are arbitrary. Simulating an 
asynchronous machine on an LMC fixes the delays in 
ordered amounts corresponding to the order of compu- 
tation. Thus, some restrictions which must be imposed 
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in general asynchronous design may be removed if the 
design is to be simulated on an LMC. 

A simple solution exists for the problem of the un- 
equal delays in the simulation of designs based on the 
premise of equal delays. Two sets of values, old and 
new, of the delay variables are stored whenever neces- 
sary. From the old values the new values are computed 
and stored. When all new values have been computed, 
these values are transcribed to the positions formerly 
occupied by the old values and become the old values 
for the next round of computation. This process would 
be repeated once per major cycle. It is more convenient 
to temporarily store the new values as complements. 
Thus, in (4) and (5b) one variable position would store 
L,, the old value, and another position would hold D,, 
the complement of the new value. 

The variables A, B, and C represent transient infor- 
mation which may be destroyed as soon as D, or in gen- 
eral the equivalent of the new value of a variable, has 
been formed. Thus, only a few variable positions are 
needed for all of the intermediate variables. A gross 
estimate on x would be thrice the number of equivalent 
one-bit delays in a simulated synchronous machine 
(old values, complements of old values, complements of 
new values) plus the maximum number of OR terms in 
any one equation (if the equations are written in normal 
disjunctive form). If an asynchronous machine is simu- 
lated, the number of variables would be roughly equal 
to twice the number of equivalent one-bit delays (old 
values, complements of old values) plus the maximum 
number of OR terms in any of the normal disjunctive 
equations. An estimate on m is harder to form since it 
depends on the particular program used. A gross esti- 
mate, probably high, sets m equal to ” plus the number 
of repeated uses of intermediate variable positions. An 
allowance also must be made for the input-output 
routine. 

EXAMPLE 1 

An example of a simple counter is programmed to il- 
lustrate the techniques. The physical problem is to 
count the number of cycles of a waveform and indicate 
every eighth one, Fig. 2. 

A rather standard synchronous circuit for this pur- 
pose is given by (7). A, B, and C form the counter of 
eight. The count sequence is 000-001-010-011-100- 
101-110-111-000-etc., with 111 taken to be the eighth 
count. G and H form an input interlock to adapt the 
asynchronous input to the synchronous circuits. On 
GH=10 a count is made, on the other combinations of 
Gand H no count is made. 


G41 = (Input), (7a) 
Hua = 'G; (7b) 
Acai = (A, B,C; + AcB, + A;-C,)-G,-H, 

ot KGesPH, Ay (7c) 


Joye = OBS; sie BiG, -G,:-H, + (G, ate H,) “BD; (7d) 
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—> — major cycle 
nu 
Pe Sse arene 
Count 7 
Output = "1" 


Count O Count 1 


Fig. 2—Counter problem. 


Ce GH, —- (G, a His @; 
Ar Bi Cyr Geule 


(7e) 
(7f) 


The variable register (R) in the LMC uses sixteen bit 
positions with the following variable assignment: P1, 
input (G,) atid. output: P2, G,; £3, A.; P4, Av P35 eB 
Po, B,; P7,.C, P8,. Ce P9, He P10, 77 Pitsthroucn 
P16, working storage with P13, P14, P15, and P16 
sometimes storing Hei, Coy, Boys, and Avy respec- 
tively. A straightforward program matrix is given in 
‘Fable LIT. 


Gani ay 
(Output) s41 


VARIATIONS 


The basic computer can be modified in several ways. 
Instead of the logical product instruction, a logical sum 
could be used. The resulting equally simple instruction 
code would form the logical sum, Z;+R;. The store in- 
struction would leave LZ; in the “0” state. At some cost 
of convenience the Sheffer stroke function (Z;+R,) 
and the Pierce arrow function (Z;:R;) could be used as 
the basic evaluating function. In these latter two cases, 
the store instruction would leave Z in the “1” and “0” 
states, respectively. 

A longer set of instructions can be included by using 
additional channels in parallel with the M-N channels. 
With a larger set of instructions additional delays, be- 
sides L, could be useful. Additional program channels 
and delays will not increase the flexibility of the com- 
puter, but would be introduced to increase solution 
speeds. 

Another variation substitutes a multiple bit time 
delay for the one bit time delay L. In the basic com- 
puter L serves as a communication link between the 


-binary variables in R. This function can also be per- 


formed by any multiple bit delay that is relatively prime 
to n. The case discussed here is for a delay of »+1 bit 
times. Denoting this delay of x+1 bits also by L, then 
(1a) for the basic computer would be replaced by 


Linn = Li:-Mi+ Li-Ri + M,N, (8) 


and (1b) through (le) would be essentially unchanged. 
An equivalent program for a given problem will differ 
in length so that a modification to (1c) and (1d) will be 
needed by way of changing the length of the M-N de- 
lays. The wired logic of the basic computer remains un- 
changed though and only the amounts of the various 
delays are changed. The instructions are the same and, 
in general, the same programming considerations apply. 

A program written for the basic computer may be 
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TABLE ESI 
PRoGRAM Matrix For SyNcHRONOUS EIGHT CouNTER CouNTING ASYNCHRONOUS SIGNALS 
eer Ay ACN Be Chl’ C ih AH. Hy, working «++ Variable 
cle 
[! 1 3 4 5 6 7 8 DEMO eld D120 A186 alte 15°28 16 Var. position 
1 Z, — = — = — — 
— — input/output 
: 2, O W O W O W O W O W — — complements 
RE ie iO Oe ee BG: 
4 — O — O O — — = O == W = =e == — A.BG.H,. 
ae OR ye Oe es) NO ay Si eel ee ee AGH, 
meer O) i= oe Oe —- = Se Wh niceties Ged} 
| 7 _ -- — O — — — —- — O —_- — — — Wie sine FAG 
: 8 RS ee te ee a OO OOO - 0 = We Ane 
9 — O oe ae a O O nas O oa W BC.GH, 
ee Oe i Oe O10 Se B.C.GH, 
) 11 a) a) Ww — = + GB, 
12 eee ee SO ee Oe ee = 
413 GeO. COL" ion WW =e ee 
14 Bere O O a O ae W at ai hams aaa CG is 
| 15 O 0 = eS WS OS ES UE 
(16 — — = — — a — O — O = am W <= == Ja(Cr 
1:7 O O O W aa am Cos 
18 O a — — — — — Ye = ies 
19 = O a O == O — O O == 
20 ie eee eT ee ee Oe oe ee a ee 0 — — — output (G.H.A.B.C.) 
ast = == = — W O == — HAsyw-H, 
| 22 O —— Css 
23 — = — a= W == = - = = O BB; 
| 24 =—Tv = W AsnAs 
* Some of the programmed equations, (7), are redundant. Minor cycles 1 and 24 may be combined. 
| 
jlapted by the mechanical process of inserting » do Cee 


othing instructions between each pair in the original 
-ogram. This means though that only one bit in L 
ever of interest and the program is therefore very in- 
ficient. A more efficient procedure involves computing 
multaneously on up to ~+1 evaluations. This limit 
sually cannot be reached due to timing difficulties; 
e., some variables cannot be computed until others 
ave been computed. A few more detailed comments are 
ven in the next example. 

EXAMPLE 2 

‘The same problem as in the previous example is 
sived with the simulated counter taken to be asyn- 
4ronous and with ZL taken to be an u+1 bit delay. 
ased on the coding of Table IV, then the following 
juations describe the counter. 


: G.41 = (Input), (9a) 

© Daa = Dy-G,+ Ay B,-Cy-G + Ae Be CoG 

) AASB GkG, BAB. Cr Gi (9b) 

We ce G4 A: BD, 

Saadeh alias (9c) 
Bei = B,C, + B,-D, + Be Gs + As:Cy: De: G; (94) 
eee CA, Det Aw G, + B,C, Da: Ge (96) 

Jutput)e41 = A,: Bs: Ce: Ds: Gs. (9f) 


ASYNCHRONOUS CODING OF AN EIGHT COUNTER 


Next State 
Present 
State Input (G) 
0 1 
ANE Bt Gan” As BaCe DD Al Be GieD) 
Op OW ROO Ope Os Ose O55 OO 
OF POR FOr st OF =0 1 if OPO Oat 
Om 1 1 One0 1 1 ORO tonne), 
Pe) it 0 0 1 im =O OF 20 10) 
0 1 1 0 Os al ili 40) il il 1 
0 1 1 1 0 1 0 1 0 1 1 1 
0 1 0 1 0 1 Oat 0 1 OFenO 
OF alk Oe Ea Dune AU) OF M7030 
1 1 Onen0 1 ier .0 1 0 1 
1 1 0 1 1 1 it 1 1 Tao) 1 
1 1 1 1 1 1 1 1 1 1 0) 
1 1 1 0 1 0 a0) 1 1 if =O 
1 0 ih at) 1 0 1 0 1 0 1 1 
i! 0 1 1 1 0 O 1 1 0 1 i 
il OpenO, 1 1 0 O 1 1 Ve Orr 0) 
1 OO tae O. OF OO 1 OO ReO 


* Output =“1” with 4,B.C,D,G, =1 or “1” with A.41B,C:D.G, =1. 


For this counter the delay of R can be taken as 15 bits 
and that of LZ as 16 bits. The program matrix and vari- 
able assignments for these equations are given in Table 
V. It will be easier to interpret if the micro-instructions 
are read on a descending diagonal to the right. For 
example, at ¢(5,2) the computation of A,-G, begins 
and resumes again at ¢(6,3) with a logical product 
micro-instruction. This particular variable in Z is po- 
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TABLE V 
PRoGRAM MATRIX FOR AN ASYNCHRONOUS EIGHT COUNTER 
: G. G A. A B B G CEASD RED: ++ working --- Variable 
Minor 4 : rs ig ®. y - 
Cyde i 2 oe 4h tt SA) 6 7 ty 8 eo nO ei, to ee Variable posed 

1 Zi input/output 

2 O = 

3 a W. G, 

4 — — O 

5 O O O — 

6 oo — O ie Dea 

7 = a Ss 

8 — — O = = W Desi 

9 O 

10 O O = ecee a es ie O — aes | 

11 = —- a 

12 = — — — O — — O oo — Ww AC. 

13 — _- — O O 

14 O O O —_ O = W = = = AD: 

iS O = 

16 O — O - O — Ww — — G:As 

17 Cie = 5 <0 WigihO. deceg caer BC. _ 

18 a — O O Vv — GBC. 

19 — O — — O — ~- — Ww — O — IID): 
20 — 0 ee 
Dil Ce O ae O O me W ae ease GB; 
2 — oO O @) O _ We 0 = SS CD; 
DS Sr a a O oe O O Ta ares O == oe rari O W a GA.GD: 
24 as —S W = O = ae, rae oo a => W —— O = PAs G:C; 
BS _ O O — — — O O O Saw e Abs 
26 = —— O W mr O O rey W ev G-A.E.D, 
27 Se ee ee ee OO) ee Oe G.D, 
28 == = = — — -- O Ww G.A.B:D; 
29 O W —_ O — A.BGAG: 
30 = == we 
31 W aT ee O W ay a Bes A BsG.Gs 
32 O — — — O — O a 
33 W O =p W oi Bea; ALB.C-G 
34 — == ai 
35 —- — O W A.B;C;G;s 
36 O 
37 me 
38 WwW Cost 
39 O Paieer 
40 W oe W G.A-B.C:Ds, Cal 

tentially modifiable at £(7,4), ¢(8,5),--+, #(15,12) formally interpret and execute a given set of machine 


but it is not. At ¢(16,13), A,-G, is stored in P13 and the 
“corresponding position” in Z is set to “1”. During this 
process the evaluation of other variables is being made 
simultaneously. For example, during minor cycle 7 five 
variables in Z are of interest. The micro-instruction 
located at ¢(6,9) pertains to the previous major cycle of 
computation. It should be studied as the conclusion of 
the program. 


DISCUSSION 


The concept of micro-programming of Wilkes>*® and 
Mercer’ enters by way of comparison. Mercer gives as a 
definition of micro-programming: 

Micro-programming is the technique of designing the 
control circuits of an electronic digital computer to 


5M. V. Wilkes and J. B. Stringer, “Micro-programming and the 
design of control circuits in an electronic digital computer,” Proc. 
Camb. Phil. Soc., pp. 230-238; April, 1953. 

6 M. V. Wilkes, “The best way to design an automatic calculating 
machine,” Proc. Manchester University Computer Inaugural Conf.; 
July, 1951. 

7R. J. Mercer, “Micro-programming,” J. Assoc. Comp. Mach., 
vol. 4, pp. 157-171; April, 1957. 


operations as an equivalent set of sequences of micro- 
operations, elementary operations that can be executed — 
in one pulse time. 

The computers of this article are micro-programmed in 
a broader sense than is given by this definition. Both the 


_ “arithmetic and control circuits” are micro-programmed 


using micro-operations. Most of the same conclusions 
as Mercer reaches may be made for the LMC. Micro- 
programmed machines do contain extreme versatility — 
and flexibility for very little complexity in the circuits. 
The LMC are serial machines (not out of necessity, how- | 
ever) and a conclusion to the effect that they are in- 
herently fast is not forthcoming. Wilkes and Mercer 
generally envision parallel machines. 

The basic LMC is not a fast machine. Solution speed 
decreases as approximately the square of complexity of 
the problem. For comparison purposes the simulation | 
of the SEAC characteristics is estimated to involve a _ 
speed differential of about 610% No conclusion as 
regards speed should be drawn about micro-program- 
med computers in general. The basic LMC can be much 
improved by more complex micro-instructions and by 
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ternate means of treating some classes of bits, such as 
e« memory bits of the simulated machine. 

{The means of handling input and output is restrictive 
it general. Since only a single information transfer line 
used, all data must be time coded. Additional input 
id output equipment may be convenient. 

|The work of Smith with digital system simulation 
iould definitely be noted in which he uses the concept 
| storing the state and logic of the simulated machine 
recorded information in the memory.’ The particular 
‘sign of Smith is more complex and specialized in type 
‘simulation but the value of simulation to reduce com- 
feuity is recognized. 


8 W.E. Smith, “A digital system simulator,” Proc. Western Joint 
ymp. Conf., pp. 31-36; February 26-28, 1957. 
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Perhaps one of the most significant results is the 
demonstration of how an arbitrary design may be im- 
plemented without single bit delays (i.e., without flip- 
flops or toggles). Conceivably, at some frequencies gat- 
ing and amplification functions may be mechanizable 
but single bit memory loops may be unreliable. 
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Analytical Design of Resistor-Coupled 
Transistor Logical Circuits” 
M. W. MARCOVITZ} anp E. SEIFt 


Summary—tThe object of this paper is to analyze and to develop 
‘sign procedures for a resistor-coupled transistor circuit used in the 
echanization of logical operations. The basic circuit consists of 
se transistor and a number of resistors. This circuit performs the 
R function followed by the NOT function or the AND function 
lowed by the NOT function. With these compound functions 
echanized it is possible to build any logical system. 

The first requirement for operation of this circuit is that the 
ansistor must be saturated if one or more inputs are low. The 
cond requirement is that the transistor must be cut off if all of the 
puts are high. A ‘‘worst case” analysis is performed for each of 
ese requirements. 

‘Three types of solutions are described and discussed: general 
rpose, intermediate (flexible), and special purpose. The general 
tpose design uses one standard stage for the mechanization of an 
tire logical system. The special purpose design is tailored to fit a 
mplete logical system. If any logical change is made, a large part 
the circuitry must be redesigned. An intermediate (flexible) de- 
m uses different stage designs, having different numbers of inputs 
d outputs. 

Procedures are developed for the design of the general purpose 
stem in two ways: hand computation with a slide rule; computation 
a small digital computer (Burroughs E101). Procedures are de- 
loped for the design of the intermediate flexible system using hand 
mputation. Examples of the design procedures are shown and 
timized for minimum required transistor base current. Circuits 
re constructed and tested to verify these design procedures. 


* Manuscript received by the PGEC, December 17, 1957; revised 
inuscript received, March 24, 1958. This paper is an abbreviated 
‘m of a thesis submitted May, 1957, to the Faculty of the Moore 
hool of Elec. Eng., University of Pennsylvania, in partial fulfill- 
nt of the requirements for the degree of M.S, in electrical engi- 


, Burroughs Corp. Res. Center, Paoli, Pa. 


INTRODUCTION 


HE purpose of this paper is to analyze and to de- 
velop design procedures for a transistor circuit 


which is capable of performing the three basic 
logical functions AND, OR, and NOT. This circuit has 
been described by Cole, Chien, and Propster! and by 
Rowe and Royer.?~4 The circuit is shown in Fig. 1. 

The base circuit must be so designed that the tran- 
sistor is cut off if all m input voltages are zero or very 
nearly zero. If any one or more of the m input terminals 
is negative, then enough base current must be drawn to 
assure that the transistor operates in the region of 
“saturation.” Its collector voltage will then be nearly 
zero. 

The m inputs are connected to the collectors of previ- 
ous stages and the n-output resistors are connected to 
the bases of succeeding stages. The stage under discus- 
sion will be cut off (OFF) if and only if all preceding 
stages are in saturation (ON) so that all m inputs are 
nearly zero. When the stage is OFF, its collector volt- 


1C. T. Cole, Jr., K. L. Chien, and C. H. Propster, Jr., “A transis- 
torized transcribing card punch,” Proc. Eastern Joint Computer Conf., 
pp. 80, 83; December 10-12, 1956. 

2W.D. Rowe, “Transistor NOR circuit,” Trans. Paper 57-195, 
presented at Winter General Meeting AIEE; January, 1957. 

3 W. D. Rowe and G. H. Royer, “Transistor NOR circuit design,” 
Trans. Paper 57-196, presented at Winter General Meeting AIEE; 
January, 1957, 

4 W.D. Rowe, “The transistor NOR circuit,” 1957 IRE WESCON 
CONVENTION RECORD, pt. 4, pp. 231-245; August 20-23, 1957. 
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-Vec 


Fig. 1—The basic circuit. 


age will be sufficiently negative to assure that all » suc- 
ceeding stages are ON. The stage under discussion will 
be ON if any one or more of the m-input stages are OFF. 
The collector of the stage under discussion will then be 
near zero volts and will not turn on any of the ” succeed- 
ing stages. It should be noted, however, that a suc- 
ceeding stage could, nevertheless, be ON if it has an- 
other input resistor connected to an OFF collector. 

P-N-P transistors will be used for all examples in 
this paper and all definitions will be in terns of p-n-p 
transistors. In order to use n-p-n transistors, it is merely 
necessary to reverse the polarities of the power supplies. 

Many of these stages can be interconnected to per- 
form logical functions. The input resistors of any stage 
are the output resistors of preceding stages, and the 
output resistors of any stage are the input resistors of 
succeeding stages. Since every one of these resistors is 
both an input and an output resistor, it will be more 
convenient to refer to them as transfer resistors (Rr). 
(The definitions and symbols used in this paper are 
listed in Appendix I). The same circuit configuration is 
used to represent either an AND or an OR merely by 
redefining for each stage whether a conducting tran- 
sistor shall represent the presence or the absence of a 
signal (7.e., the signal or its inverse). 


DESIGN CRITERIA 


A transistor can be in one of two states: OFF or ON, 
or in transition between them. By definition a tran- 
sistor is said to be OFF if the emitter is reverse biased 


such that I¢=Ico. Ico does not change appreciably . 


when the base-emitter reverse bias exceeds some mini- 
mum value, but is less than the breakdown voltage of 
the base-emitter junction. This minimum value is a 
characteristic of the transistor type to be used for a 
particular design and is available to the designer from 
the manufacturer’s specifications or from test data. 
Ico is similarly available. The first design criterion is 
that the base voltage of an OFF transistor must be 
equal to or more positive than some minimum required 
value obtained from transistor characteristics. The base 
voltage of a cutoff transistor is said to be “up” (i.e., in 
the positive direction) and the symbol to be used is 
+ Vau (4.€., Vase, up)» A minimum limit is indicated by 
a line below the symbol, a maximum limit by a line 
above the symbol. The minimum required value of base 
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Fig. 2—Collector characteristics. 


voltage for an OFF transistor is therefore + Vau. (Vau 
merely represents the magnitude.) The first design cri- 
terion is therefore expressed by: Vgu>Vzgv, where 
Veu is known from the characteristics. 

When a transistor is ON it is required to be “satu-_ 
rated.” For the purposes of this paper and by definition, 
a transistor is said to be saturated when its collector 
voltage is equal to or less negative than a value, — Vey 
(Vooltector, Up). Lhe region of saturation is shown shaded 
in Fig. 2. The second design criterion is that the col- 
lector voltage of an ON transistor must be equal to or 
less negative than some maximum permissible value, 
Vou <Veu. 

Other criteria may be imposed by the transistor 
specifications, such as maximum allowable collector 
voltage, emitter current, power dissipation, etc. 


Worst CASE ANALYSIS 


Every design must satisfy the two design criteria, 
Vau>Vev and Veu<Veu, in order to assure satisfac- 
tory circuit operation. It is insufficient to know that the 
design center values of a particular design satisfy these 
conditions. All resistance values and voltage supplies 
may vary between maximum and minimum limits. A 
reliable design requires that each design criterion be 
satisfied even when all circuit parameters have simul- 
taneously deviated by the maximum expected amount 
in the direction making it most difficult to satisfy the 
circuit criterion. A circuit in which all parameters have 
deviated in this manner is called a “worst case.” Since - 
there are two criteria there are two worst cases: 

1) The worst case for Vgv. In this case every circuit 
parameter has deviated by the maximum expected 
amount in that direction which will tend to make 
Vu smaller. If for this worst case Vgy is set equal 
to Vauv, then Vey will be larger than Voy for all 
other circuit conditions. Therefore, Vgy>Vazv. 

2) The worst case for Vey. In this case every circuit 
parameter has deviated by the maximum expected 
amount in that direction which will tend to make 
Vou larger. If for this worst case Voy is set equal to 
Vou, then Vey will be smaller than Veu for all 
other circuit conditions. Therefore, Vou <Vev. 

If a circuit satisfies both worst cases it is reliable and 
will never fail to meet the design criteria unless one or 
more parameters deviate more than their maximum 


<pected deviation. This is classified as a component 
iilure. Circuit diagrams describing the two “worst 


ase” conditions will be derived. For each worst case it 
ly be necessary to determine the direction in which 
ch circuit parameter must deviate in order to con- 


‘ibute to circuit failure. 


Vorst Case for Vau 


The results of this analysis are shown in Fig. 7 and 
»ference may be made to that figure for the following 
aragraphs. 

The output variable is Vgu and it is necessary to 
kamine each independent variable to determine wheth- 
- that variable should be a maximum or a minimum 
» order to make Vgu; smaller. This can be done by ex- 
ressing Vgui as a function of all independent variables 
» the circuit and obtaining the partial derivative of 
ach variable with respect to Vaui. If this derivative is 
ositive for a particular variable, then that variable 
sust be a minimum to make Vy; smaller; if the deriva- 
ve is negative, then that variable must be a maximum 
>» make Vgui smaller. For a fairly elaborate circuit this 
rocedure can be very laborious. It is possible to exam- 
1e the circuit one branch at a time, starting at Vgu., 
ad thereby greatly reducethecomplexity of solution. This 
rocedure is followed here. The first variables to be ex- 
mined are Vggi and Ra. For this purpose the rest of the 
ircuit, which is connected to Vgyi through Ico, Ir, and 
m—1) Ir, Fig. 3(a), can be replaced by an equivalent 
ircuit, Fig. 3(b). EZ, may, in general, be either positive 
, negative. The equivalent circuit is independent of Ro 
nd V zai. Ia Afn( Rai, Vea), TRE Afn( Ra, Vzpi)- If Io3 
ows in the direction shown, (Vgz1—£.) >0. 


Ro 
Vaui1 = Veszi — Roce R: (Veni — E.). (1) 
OV R 
ee 21 6 (2) 
OV zai Roo Re 


‘herefore, Vgui decreases as Vga: decreases. Since for 
he worst case for Vgui all variables deviate so as to 
take Veui smaller, Vez, should decrease and Vm is the 
rorst case. (This is only the effect of Vgz acting on the 
ircuit through Ry. The effect of Vgz acting through 
’» and R23 remains to be investigated.) 
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<r 
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(b) Fig. 4—Equivalent circuits for I71 and Ir. branches. 
Fig. 3—The Vay; node. 
OVeu1 R. 
= — (Vga1 — E.) ————— <0 (3) 
OR (Roi + Re)? 


. Ry is the worst case for Vgui. 

Since the worst case deviations of Ra and Vgz: have 
been determined from an equivalent circuit which repre- 
sents the complete circuit, Ra and Vgm can be consid- 
ered as fixed from here on. 


Now: 
Veau1 = Var — Rola (4) 
OVeu1 ses 
St 8 ae Roy << 0 
OTe (5) 
Vea, Rau. 


The worst case for Vgm requires I. 


To. = Ico. + Tru + (MH — 1)T re. (6) 


Ico1 = Ico, a known transistor characteristic; 7 = m. 


The Jri1 branch is not connected to a terminal point of 
the over-all circuit and therefore requires equivalent 
circuits at both ends of Rru. 
From Fig. 4(a): 
Fez a; Ee 


Tig oe Danae SE (7) 
i pul aR ota on 


From Fig. 4(b): 


Oe (8) 


Tru implies Rri. I7i2 implies Rriz and Veur. 
Vou = Veun = Vou1; and Inu = Ira =I 71. 


At this point it may be convenient to generalize on 
the procedure. The circuit branch under investigation 
is shown unchanged, with the rest of the circuit repre- 
sented by equivalent circuits. In each case a variable is 
to be maximized or minimized as a function of devia- 
tions, not of the equivalent circuit, but of parameters 
of the branch under discussion. The variable is either 
the output variable or an intermediate variable deter- 
mined from previous considerations (such as [712 above). 
Since the circuits are very simple, the required devia- 
tions are easily determined mathematically. For exam- 
ple, if the current through the circuit is to be maximized 
(as in the [712 circuit above), then the series resistance 
must be minimized and the voltage difference between 


112 
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Fig. 5—Simplified Vgps circuit. 
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Ee Re Sie Ra = 
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Veo: a 


Fig. 6—The Vezgni node. 


terminals must be maximized. Hence Try implies Rrie 
and Vou: 

Since the worst case deviations of the circuit up to 
this point have already been determined, the rest of the 
circuit can influence the output variable only through 
its influence on J71. 


Tri = Iu + (@ — 1)T 72 — Leu. (9) 


ZT implies Ru and Veo, where Voc is the collector sup- 
ply voltage as it influences the circuit through the Ry 
branch. Veco also acts on the circuit through other 
branches and its influence on the over-all circuit there- 
fore has not yet been determined. 

Tr. implies Rr2 and Vgps, where Vgp3 is the maxi- 
mum voltage at a conducting base which has one input 
resistor connected to a conducting collector and (m—1) 
input resistors connected to cutoff collectors. A simpli- 
fied circuit for the evaluation of Vgp: is shown in Fig. 5. 


Tor = Bim. 


B,=B, a transistor characteristic. 

To evaluate the worst conditions in the J, branch, an 
equivalent circuit is used to represent all of the rest of 
the circuit as viewed from Vgp: (Fig. 6). 


For Ip; the circuit requires Rg; and Vai. 
fs: = Irs — (m — 1)Irq — Too. (11) 


Irs implies Veuve and Rra; Tx implies Vggo and Ro»; Ins 
implies Rr3 ,and with Rr; determined: 


(12) 


Icoz=Ico, a transistor characteristic. Jy: implies Voce 
and Ryo. I 76 implies Rr. 


Ir6 = Io + Tog + (m — 1)I 75. 


Irs = Ip Tco2 — (nw — 1)T 76. 


(13) 


IRE TRANSACTIONS ON ELECTRONIC COMPUTERS 


(10) 


June 


+VBB3 

I 

331) Res 

-Vev2| »/ £82 e 
Ip2 Vee) = 
At ON 
m-i 
Tim-s)t75 


ih ~Vcus 


Ite} The Tre *Vep2 | 
esos otee Rez Ta 
(n-1) It6 : Rai 

R 
% %It3 x Le whip 
-Veo1 <—=Veo1 | b / RB! e  -Vcul m= 
CERMelces = 6 
A BS 
RT3 Ig) Vel — (m-) Iq, 
c 
ON 
; Veur RT 
Rra a <a tVaul 
m-! It, 
Ri2 
I R OFF 
ON) sletiniyzah Ru Rte 
n=1 
o o J Jo-n tre 
-Vcc2 -Vcu2 
= — Ti 3 ut 
Vee! -Vep3 


Fig. 7—The worst case for Vay. 


Tg2 implies Rg and Vago. I23 implies Vez3 and Ros. Irs 
implies Vcv3 and Rr. 

In this analysis the worst case conditions for Vgg 
were found to be Vgm, Vas: and Vgzs; for Veco they 
were found to be Vee: and Vcc2. The voltage supplies 
are, in general, common to the whole circuit. If Veo and 
Voc were used simultaneously, the resultant circuit re- 
quirements would be “worse than the worst case,” which 
is inefficient, especially since the amount of extra safety 
would be unknown. More reliable circuits are designed 
if the worst case limits of the voltage supplies are an- 
alytically derived. For the worst case for Vv, the limit 
values are Vggz and Vcc. The worst case circuit for Vag 
is shown in Fig. 7. The worst case circuit for Vey can be 
derived in a similar manner and is shown in Fig. 8. The 
two circuits shown in Figs. 7 and 8 differ in only two 
circuit parameters: Rr; and V zp. 


Basic EQUATIONS 


The worst case circuits (Figs. 7 and 8) yield the two 
basic equations for a general purpose system. 

The basic equation for the worst case for Vgy is 
written about the Vgu: node of Fig. 7. 


Ico. = [(1 — 0) Ven — Vewil (4 — p)G2 


— Veur + Veui)m(1 + p)Gr (14) 
where 
I 
G Sens 
R 


Veu1 is an intermediate variable. Parametric equations 
can be written to define Vey, by completing the circuit 
solution. 


ia 
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=Vcc! -Veos 
Fig. 8—The worst case for Vev. 


The basic equation for the worst case for Vcv is 
ritten about the Vevi node of Fig. 8. It can be shown 
rat Veuun=Veviz. and will hereafter be referred to as 
| 


'CUl1- 
[1 — o)Veo — Veusl(1 + p)Gi 
+ [Viens — Veus)(1-+ p)(n — 1) 
— (Veui + Veui)(1 — p)|Gr. (15) 


‘The parametric equations required for the solution 
the above equations were derived from the two worst 
ase circuits. The following approximations were made: 
epie= Veni; Veow= Ven. Ven and Vepix are the 
oitages at the collectors of a cutoff transistor under 
orst case conditions. See Fig. 7 and Fig. 8. 

From Fig. 8: 


Za 
: 


Marcovitz and Seif: Resistor-Coupled Transistor Logical Circuits 


(1+ 0)Ven(t — p)G2 — Veui[(1 — p) + (m — 1) + p)|Gr — Toa 
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From Fig. 8: 
ler = Bly; (20) 
we ResGiGr(m — 1)(1 — o V 
Vandauiabs eee ee ee 
(LS p)Gs “pi tie Gr 
From Fig. 8: 
Veni = Vai aie TpiRai. (22) 
From Fig. 7: 
py, Bem swGrhll so coll pCa feos (23) 


[(m—1)(1+p)+(1—p)]Gr+(1—p)Gr 


Substitution of the parametric equations in the two 
basic equations would result in complicated and non- 
linear expressions. It is easier to handle the basic equa- 
tions and parametric equations separately. 


TYPES OF SYSTEMS 


There are three general types of systems: the general 
purpose system, the special purpose system and the 
flexible, intermediate system. The choice of system is 
determined by the nature of the application. 


The General Purpose System 


In this system a standard transistor stage is used. In 
some stages not all of the inputs and outputs are used; 
the unused transfer resistors are grounded. It is not 
necessarily required to use a general purpose stage with 
the maximum number of inputs and the maximum 
number of outputs to be found anywhere in the system. 
This may turn out to be unnecessarily wasteful of com- 
ponents, power, and space. It is often advantageous to 
choose a general purpose stage which will take care of 
the great majority of circuits and to use two stages in 
parallel (or emitter followers, or higher current gain 
transistors) for the few remaining circuits. 


The Special Purpose System 


In this system all of the logical circuitry must be de- 
signed by the logical designer before the electrical circuit 


(16) 


Veui = 


pe kay 
VenieVeun = (see Appendix I). 


Bl 


From Fig. 8: 

o = [Veni — Vani(1 — p) — Veil + p)(m — 1)|Gr 
a [Vep1 aid et o)Vez\(1 + p)G2. (18) 

From Fig. 7: 

a1 = [Veni — Vani) (1 — p) — Veoi(1 + p)(m — 1)|Gr 
Fi [Veni te (1 nr o)V eel (1 aie p)Go. (19) 


(1 — p)Go + [(1 — p) + (m — 1)(1 + p) Gr 


(17) 


design can begin. The output requirements must also 
be known. Each transistor stage is designed for optimum 
performance in the particular logical configuration in 
which it is employed. If any change is made in the logic 
after the circuit design has been completed, a large 
amount of circuit redesign might be required. This sys- 
tem is not covered in this paper because of its special 
purpose nature. 
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The Intermediate, Flexible System 


An intermediate, flexible system has been developed 
which incorporates some of the advantages of both the 
special purpose and general purpose systems. In any sys- 
tem, logical elements will generally have different 
numbers of inputs and outputs. Circuit designs are de- 
veloped for the various logical elements required. As- 
sembly of these designed circuits into a system may 
require some modification as a function of the logical 
configuration. These modifications require a small 
amount of simple computation. An intermediate, flexible 
system can easily be changed to comply with required 
logical changes. 


DESIGN PROCEDURES 


Two basic equations (14) and (15) have been devel- 
oped for the general purpose circuit to satisfy the two 
basic criteria. The 15 independent variables in these 
basic equations and the associated parametric equa- 
tions, (16) through (23), are: two power supply voltages 
(Voc, Vez); power supply tolerance (oc); three resistors 
(Ki, Re, Rr); resistor tolerance (p); six transistor pa- 
rameters (B, Ico, Rz, Vz, Vov, Vv); number of inputs 
(m); number of outputs (7). Solution of the two basic 
equations will determine only two unknown variables. 
Consequently, some other means of determining 13 vari- 
ables must be established if the solution is to take this 
form. 

Many of the independent variables may be assigned 
values based on practical engineering judgment and on 
the requirements of other parts of the device for which 
the logical circuitry is intended. The selection of a tran- 
sistor type gives values for six of the variables. Selection 
of power supply voltages and tolerance, and resistor 
tolerance gives four variables. A knowledge of the logic 
of the circuit gives values for m and n. If all these values 
are assigned, the three resistors remain as the only three 
unknown independent variables. 


General Purpose System 


The following are some design procedures for a general 
purpose system. 

Assume Values for Circuit Parameters: One design 
procedure is to assume values for all the circuit parame- 
ters except two. The two basic equations would then be 
solved for the two unknown parameters. If there is no 
solution new values would be assumed and a solution 
tried again. By assigning the independent variables as 
suggested above and by assigning a value to one of the 
resistors, the problem is reduced to the solution of the 
two basic equations with two unknowns. If there is no 
solution with the value of resistance assumed, a new 
value of resistance must be chosen and the solution re- 
peated. There may be an infinite number of solutions, 
or there may be no solution, depending upon what 
values were assumed for the other circuit parameters; 
e.g.. m and n too high. 
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Assume Values for Intermediate Variables: The assign- 
ment of values to intermediate variables (e.g., Im, 
Veni) gives additional basic equations. These values 
would be selected on the basis of engineering judgment. 

Optimum Solution: Since the general purpose solution | 
involves the solving of two basic equations for the de- 
termination of more than two independent variables, 
there may be an infinite number of solutions. By estab- 
lishing additional criteria it is possible to determine a_ 
more desirable or optimum solution. Some of the cri- 
teria that might be used for optimizing are: minimum 
transistor base current; minimum transistor power dis- 
sipation; minimum required transistor current gain (B); 
minimum circuit power dissipation; optimum circuit 
for switching time considerations; maximum number of 
inputs (m); maximum number of outputs (7); maxi- 
mum resistor tolerance (p); maximum power supply 
voltage tolerance (a). Illustrations in this paper use 
minimum base current as an optimizing criterion since 
it will lead to low-transistor power dissipation and low- 
power supply requirements. 


Special Purpose System 


The above design procedures can also be used for the 
design of a special purpose system. Each logical ele- 
ment requires its own design. Since the load on the 
stage is known, a third equation is supplied. 


Intermediate, Flexible System 


The above design procedures for a general purpose 
system can also be used in designing an intermediate 
system. The two basic criteria must be satisfied. The 
load on a stage is known, giving a third equation. Ven: 
should be a constant value for all circuits. Conse- 
quently, the assignment of a value for Ven gives a_ 
fourth basic equation. 


COMPUTER SOLUTION OF A GENERAL 
PURPOSE SYSTEM 


A solution of a resistor-coupled transistor logical ele- 
ment has been performed on a Burroughs E101 digital 


‘computer. An optimum solution was obtained with min- 


imum transistor base current as a criterion. 

The two basic criteria give the two basic equations 
(14) and (15). To optimize for minimum transistor base 
current (Iz1), values are systematically assigned to Ip. 
Consequently, the equation for Ig: (18) is used as a 
third basic equation for this solution. A partial solution 
is performed to determine a minimum value for I. for 
which the two basic criteria are satisfied. The three 
conductances G;, Gz, and Gr are solved as unknowns. 
All other independent variables are assigned values. 

In order to simplify the solution, values are initially 
assigned to the intermediate variables Vans, Ven, and 
Veui, since they have second order effects on the circuit 
solution. This eliminates the need for using parametric 
equations (16), (17), (19), and (21) in the solution of 
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Fig. 9—Behavior of computer solution. 


14), (15), and (18). Safe values are assumed as follows: 
"pps=2V es; Veur= Vou; Vaui= Varn. The equations 
pr Ic: (20) and Vgpi (22) are simple to calculate since 
Bi is an assigned value. However, the equation for 
“op: (23) is a function of G, and Gr. An iterative rou- 
ine is used to solve (14), (15), and (18) with the para- 
retric equation for Ven: (23). This involves an initial 
ssumption of a value for Vep; and the direct solution 
f (14), (15), and (18) for G,, G2, and Gr. Vepy is then 
alculated using (23). The computer selects a new value 
or an assumed Vep: depending on the calculated Vepy 
*hich resulted from the previous trial value of Veni. 
‘alues of Ven: are assumed such that the calculated 
alue of Veni approaches Veni. This process is con- 
inued until it is determined if a solution exists for the 
ssigned value of Jz:. The exact solution is defined by 
he intersection of a hyperbola and the line Vep1= Veov 
Fig. 9). If a solution does exist, a lower value for Ip: is 
ssigned and a new solution performed. This procedure 
s continued until a minimum value for Ig: is deter- 
ained. The exact optimum value is defined as the point 
{ tangency of the hyperbola with the line Vep1 = Vepv. 
“Using the calculated minimum J, values are cal- 
ulated for Veui, Vev1, and Vgp3 with equations (16), 
17), (19), and (21). With these values a more accurate 
Ainimum J; is determined by using the computer solu- 
ion described above. 
When the more accurate value of minimum Jz: is 
letermined, more accurate values of Vgui, Veui, and 
7ep3 are calculated. With these values the solution is 
epeated to determine exact values for the three con- 
luctances Gi, Ge, and Gr. 

A number of circuits were designed using this com- 
yuter solution and the effects on the circuit design of 
he variation of circuit parameters were studied. 


HAND SOLUTION OF A GENERAL PURPOSE SYSTEM 


This section describes the design of standard tran- 
istor stages suitable for use in a general purpose sys- 
em. The procedure will be to select a transistor type, 
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Logical Circuits 
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(i+e)Ry = Veo, 
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=(I- Veg 
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Fig. 10—Worst case circuits for manual solution. 


voltage supplies, resistor and voltage tolerances, and the 
number of inputs per stage (m). The unknowns are the 
three resistors (Ri, Re, Rr) and the number of outputs per 
stage (7). The solution will be optimized by minimizing 
the required base current (Ig). This will be done by 
assuming various values of Js: and calculating the max- 
imum number of outputs (#%) for each value of Ig. In 
this way a chart of % vs Img is generated. From this 
chart the minimum current corresponding to the de- 
sired value of 1 is selected and the final circuit is cal- 
culated. The procedure is illustrated with a sample de- 
sign. 

A transistor type is chosen so that the following be- 
come known: Jeo, B, Vau, Vev, and the worst base and 
collector characteristics. In addition, p, o, Vas, Vcc, 
and m are chosen. This leaves four unknowns: Ri, Re, 
Rr, and n. Since there are four unknowns but only two 
basic equations, two intermediate variables are chosen 
in order to increase the number of available equations 
to equal the number of unknowns. The intermediate 
variables chosen are Ig; and Venu. 

It is now possible to draw the worst ease circuits in 
simplified sections and to write equations from them. 
Fig. 10(a)-10(c) describe the worst case for Vou. Fig. 
10(d) describes the worst case for Vav. 

Eqs. (24) to (27) are written from the four circuits. 


Veni — Vapi 


(1 =f p)Rr 
K) (1+ ¢)Vee+Ven1 Vepnilm — 1) Is (24) 
(1 — p)Re (1 — p)Rr 
(1 — «)Vee — Veni 
Cp) Ri 
V epic VEpu ea Veni(n — 1) 
She ape eae ihe eee IPS 
(1+ p)Rr (1 — p)Rr 2) 
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vn (1—«o)Vee—Veu | (Vans — Veu)(n — 1) 
aa (1 — p)Ri (1 — p)Rr 
et Veu + Vou (26) 
(1+ p)Rr 
Veu + Vev)m _ (iio) Vas ey ae (27) 
(1 — p)Rr (1 + p)Re 


These equations correspond to (18), (23), (15), and (14). 
In order to facilitate computation in the preliminary 
steps of the solution, it is advisable to make two approx- 
imations: 


Vou = Vev. 
In Fig. 10(c): Ze: = In, so that (26) can be written: 
; (a == «)Vec = Veu 


_ 28 
fe Ci ip) Ki oe 


Choosing the 2N135 for the sample design and as- 
signing values for m, the voltages, tolerances, and letting 
Ip, =0.4 ma for the first trial, the following parameters 
are known: 


Vee = 19 volts; Veg = 10 volts; o = 0.05; 
[= 0.05; Ico = O05 ma, Veu = 0.15 volt; 
Veu = 0.1 volt; m = 5; Ipi = 0.4 ma; 

Ic, = 6.9 ma (using a 0.9 safety factor); Vapi = 0.26 volt. 
Only the last three values are functions of Ig; and sub- 
ject to change for succeeding trials. The other values 
remain fixed. 

The values listed above are substituted into (24), 


(25), (27), (28); Ri, Re, and Rr are eliminated to obtain 
the relation: 


38.35 


n = 14.48 — 0.682V cp: — (29) 
CD1 
To determine the maximum value of n, 
on 38.39 
= 0 = — 0.682 + (30) 
OVeni CD1 


Ven1=7.5 volts, and »=4.2. With Ip,=0.4 ma and 
m=5, the maximum value of m is 4.2. The calculations 
are repeated with various values of Ig: to obtain Table I. 


TABLE I 
m=5 
In 0.8 0.4 0.3 0.25 0.2 Oil 
n 4.23 4.2 4.13 4.06 3.91 2.99 
Veni Uy U8 7.38 HDS Ut 6.99 


The value of ” drops sharply as Ig: approaches To. 
The highest integral value of m which can be obtained 
is 4 and the lowest base current required is 0.25 ma. 
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For 
Ip; = 0.25 ma: Ici = 4.38 ma. 
Veni = 0.233 volt. 
Veni = 7.28 volts. ° 
From (28) 
1—o)Veo — V. 18 — 0.1 
ee =4.3K. (31) 
(1 — p)Le1 (0.95) (4.38) . 


From (24) and (27), Re=58.55 K; Rr=12.95 K. With 
these resistance values the approximations made before 
can now be corrected. More accurate values can be cal- 
culated for Vgui, Vou, and Vgp3, and (26) can be used 
instead of (28). Using (16) to calculate Vgu, Vaui 
= +0.195 volt. In order to find Vou, it is first necessary 
to calculate 7z;. The only difference between J»; and — 
Ip: is due to the change in Vaz. 


1 
Alg, = — R, AV ze 
Var 5 Ves 1 
Li af —————— = 0.25 + ——_ 
ae Rs 58.55 
= 0.267 ma. (32) 


From the worst collector characteristics: 
At: 


Ip = 0.267 ma 


\ Veu => 0.084 volt. 
Ig = 4.86 ma 


Using (21) to calculate Vans: 


Vzps = 0.36 volt. (Vz3 = 0.26 volt; Ra; = 0.04 Ke 


Eqs. (24) through (27) are now solved, using all 
values corresponding to an Ig; of 0.25 ma, to obtain: 


R, = 4.35 K; Re = 62.3K; Rr = 13.2 K; n = 4.09, 
The result obtained is a general purpose stage with 


m=S and n=4. It has been optimized to use the mini- 
mum possible base current. If the logical designer re- 


_ quires a standard stage with a higher number of inputs 


and outputs, then a transistor with higher current gain 
is required. 


HAND SOLUTION OF AN INTERMEDIATE SYSTEM 


There is a need for logical systems composed of stages 
having various combinations of m and n. These stages 
must be compatible, 7.e., it must be possible to connect 
them to each other. Such circuits can be designed by 
adopting a standard value for Vep:. The base circuit of 
each type of stage is determined by the number of in-_ 
puts (m) and the number of outputs (7), not by the load 
which each output is required to drive. The maximum 
load which any following stage may present to the col- 
lector of the driving stage is defined (for the whole Sys- 
tem) as the maximum allowable load current. This 
maximum load current may be required by various 
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mbinations of m and n, and any circuit which requires 
is maximum allowable current is referred to as a 
1xximum load circuit. 

It is advantageous to let the standard Vep: equal the 
‘py. Of the general purpose circuit (m=5, n=4) de- 
ned above, and to let the maximum allowable load 
rrent equal J7; of that circuit. 


Veni = 7.28 volts, Irmax = 0.508 ma. 


sing the same values as in the general purpose solution 
r Voc, Ves, p, 7, Ico, Vev, Vou, the transistor charac- 
ristics, and assigning m, a maximum load circuit is 
signed as follows: 


R, is calculated from (25). 

Ic is calculated from (26). 

Ig, and Vgpi are obtained from the transistor char- 
acteristics. 

From Fig. 10(a): 


Veni — Vapi 
ce 5 paren ee (33) 
(1 + p)Irs 


Rr is calculated. 
Re and m are calculated from (24) and (27). 


sllowing this procedure for n=2: 

‘R, = 9.0 K; Jc: = 2.1 ma; 

ei = 0.125 ma (with a 0.9 safety factor); 

0.20 volt; Rr = 13.2 K; m = 8.2 (use m = 8); 
[Re = 41.2 K. 


WwW 
ic} 
I 


nis circuit is compatible with the m=5, n=4 circuit 
id is to be used for every stage in the system where 
=8 and m=2. Other maximum load circuits can be de- 
sned for n=1,72=3, and ~=5. If a system consists 
itirely of various maximum load circuits, it can be as- 
mbled without any further calculations, for every 
age will drive its maximum load and Vep; will be the 
andard value. 

The next circuits considered are those which require 
ss drive than the maximum load circuits, e.g., a stage 
ith m=2, n=2, or a stage with m=1, n=4. Here m is 
iown. Ri, Ic1, and Vzgp: are calculated as before. Eqs. 
4) and (27) are solved for Rr and Rs. The following 
alues are obtained for m=2, n=2: 

) R; = 9.0K; Rr = 31 K; Re = 132 K. 

he load which this circuit presents to the stages driv- 
g it is 0.218 ma, from (33). Wherever this circuit is 
sed in the system, the difference between 0.218 ma 
id the maximum permissible load (0.508 ma) should 
> taken up by a collector load resistor to ground, in 
‘der to maintain Ven at design value. 

‘Similar circuits can be calculated for any combination 
‘mand n up to those corresponding to maximum load 
rcuits. 
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Test RESULTS 


General purpose circuits were designed, using both 
the computer solution and the hand computation solu- 
tion. For the same requirements and the same tran- 
sistor the two solutions gave almost identical results. 
The circuit was constructed in its worst case condition, 
using limit values of all parameters and a “worst” ac- 
ceptable transistor, in order to check the failure points. 
Circuit failure occurred exactly as predicted by the 
analysis. 

A binary adder/subtracter circuit was built using the 
hand computation solution. Many special purpose con- 
trol circuits have been built using this solution. They 
have been in operation for six months without failure. 


CONCLUSIONS 
General Conclusions 


With the procedures discussed in this paper, reliable 
resistor-coupled transistor logical circuits can be de- 
signed and built. Furthermore, it has been shown that 
the circuits can be optimized according to certain cri- 
teria. The solution can be programmed on a small digi- 
tal computer (E101). Circuits for the general purpose 
system and intermediate system were designed, taking 
worst case conditions into consideration, and with mini- 
mum required transistor base current (m1) as a circuit 
design criterion. The validity and accuracy of the de- 
sign procedures were proven by building and testing 
typical circuits. 


Circuit Considerations 


The design procedures discussed and the circuits 
tested produced much useful information concerning 
the effects of certain circuit parameters on the circuit 
solution. 

It was found that certain intermediate variables 
(Vow, Vari, and Vgp3) have second-order effects on the 
circuit solution. It was found that the collector load is 
essentially equal for the two worst cases, and therefore 
the transistor collector eurrent remains essentially con- 
stant. It was found that the current through the 
transfer resistors (Rr) connected to the collector of a 
conducting transistor is negligible. These facts were 
used to simplify the circuit design procedures. 

The effects on the circuit design of the variation of 
circuit parameters were studied by observing their in- 
fluence on the minimum required base current ([z1), 
using a computer program. Doubling the value of the 
maximum transistor collector cutoff current (Ico) more 
than doubled Js;. Doubling the value of the minimum 
large signal current gain (B) decreased I by a factor 
of 5. Decreasing the power supply tolerance (a) from 
+5 per cent to +2 per cent decreased Jm by a 
small amount, however, decreasing the resistor tol- 
erance (p) from +5 per cent to +2 per cent, as well as 
a, caused Ip; to decrease by 50 per cent. A decrease of 
the number of inputs (m) or outputs (z) causes J, to 
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decrease, however (m) has the greater effect. Doubling 
the collector power supply voltage (Voc) causes Ip: to 
decrease by about 50 per cent. A change in the base 
power supply voltage (Vzgz) has little effect. Lowering 
the minimum required base voltage (Vgy) lowers Tn. 
Raising the maximum collector voltage under saturated 
conditions (Vcv) raises I:. Those variations which 
cause a decrease in J: will have a favorable effect on 
other circuit considerations, e.g., maximum number of 
circuit inputs (m), maximum number of circuit outputs 


(1). 


Comparison of Computer and Hand Solutions 


A comparison of the computer and hand solutions in- 
dicates that each method has certain advantages. 

It is possible to calculate a solution to about the same 
accuracy using either method. The chances of making a 
computational error on a computer are remote. The 
transistor is approximated by equivalent circuits which 
give simple mathematical expressions for use in the 
computer solutions and will lead to some small inac- 
curacies in the calculated results. In a hand solution the 
transistor characteristics are used directly. 

The time to complete an optimized hand solution for 
one set of circuit conditions, using the methods of this 
paper, was about four hours. The time to calculate the 
same solution on the E101 computer was about one 
hour. However, it took the equivalent of about six full 
working weeks to obtain the computer program. If it is 
desired to obtain a small number of circuit solutions, 
the hand solution methods are satisfactory. If it is de- 
sired to obtain a large number of solutions and to study 
the effects on these solutions of varying circuit parame- 
ters, then the investment of the necessary programming 
time has been worthwhile. 

With a hand solution, the engineer does the compu- 
tation and decision making. He therefore gains a better 
understanding of the circuit operation. The computer 
solution does the decision making as well as the com- 
putation. 

It might be advantageous to place the hand solution 
on the computer, using the computer to do all the cal- 
culations and permitting the engineer to make the de- 
cisions. 


APPENDIX | 


The following definitions and symbols are used in this 
paper. 


Vze=base supply voltage. 
Vec=collector supply voltage. 
1/G:=R,=collector load resistor. 
1/G2,= R2=resistor connecting a base to Vaz. 
1/Gr = Rr =transfer resistor connecting a collector to 
a succeeding base. 
Vau = Vease, up. Base voltage of a cutoff transistor. 
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- which correspond to the worst case for Vgui. In Fig. 12 
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BASE INPUT CHARACTERISTIC EQUIVALENT CIRCUIT 


Fig. 11—Base input characteristic equivalent circuit. 


Vou = Vooltector; Up- Collector voltage of a conducting 
transistor. 
m=number of inputs to a base. 
m=number of outputs from a collector. 
Iz =base current. 
Vep=V Base) Down Base voltage of a conducting 
transistor. 
Vev= Vooltector, Down. Collector voltage of a cutoff 
transistor. 
Ic =collector current. 
Ico=cutoff collector current, emitter open cir- 
cuited. 
B=large signal current gain I¢/Iz. 
p=resistor tolerance. 
o=voltage supply tolerance. 


It will be convenient to represent the forward-biased 
base input characteristics of the grounded emitter con- 
figuration by an equivalent circuit (Fig. 11). 


Vz=base input equivalent voltage. 
Rz=base input equivalent resistance. 


The symbol V stands for the magnitude of a voltage. 
Polarity is indicated by + or —. 


APPENDIX II 
Calculation of Vou. 


Vou is the maximum collector voltage of the stage 
being designed. Jz: has been defined as the transistor 
base current flowing under worst conditions for Vey. 
Veu is the collector voltage and I; is the base current 


the two characteristics for Jg=Ig, and Ip=Ip, are 
shown. The load line is shown as a horizontal line at 
I¢=Ic1=BI,; in the region of saturation because the 
collector currents flowing in the two worst cases are 
almost equal. The actual value of the collector voltage 
under worst conditions for Ver. (Is=Iz1) is Vous (Fig. 
12). The collector characteristic for Iz =I: can be ap- 
proximated in the region of saturation by a straight 
line between the origin and the intercept of the Zz =Jp1 
curve with the Ve=Veu1 line. The approximation 
Diaries) Sit 


Veu1 = Veui = Veu1 — 
blpi Ip 


will give a value of Vou1> Vem which is on the safe 
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Fig. 12—Collector characteristics for calculation of Vey. 


le. Since Jz: is very close to Im, the straight line as- 
mption of the collector characteristic in the region of 
turation gives a good approximation of Veuix. 


APPENDIX III 


In the computer solution discussed, an iteration rou- 
ie is used. A value is assumed for Vep: and a circuit 
ution calculated. A value for Vepi’ is then calculated 
»m the circuit solution. This process is repeated with 
ferent assumed values until Vepi'= Veni. The be- 
vior of this solution is described by a hyperbola (Fig. 
where Vepi' is expressed as a function of the as- 
med value Veni. This equation can be obtained by 
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solving (14), (15), and (18) simultaneously for G, and 
Gr as functions of Von:. The results are substituted in 
(23) to obtain Veni’ (the calculated value) as a function 
of Ven: (the assumed value). This equation takes the 


form of an equilateral hyperbola where the y asymptote 
is defined by 


We (Ie), (34) 
the x asymptote is defined by 
nV (6) us (a) 
k= mVep1 + Vou + Veur) — {= ~ + i} » (35) 
tB1 


and the distance /] from the vertex to the intersec- 
tion of the asymptotes is defined by: 


12 nVeo (Ico. 
ee, 1 
gents a co} B Ge ) 


+ 2m[pVep1 + (20 + ¢)(Veur + Vews)]} . (36) 


ACKNOWLEDGMENT 


The authors wish to acknowledge the guidance, help, 
and encouragement of Dr. Howard E. Tompkins, who 
suggested the basic circuit for this paper and served as 
faculty supervisor for the thesis on which this paper is 
based. 


On the Analysis of Sequential Machines” 


R. G. GILLESPIE} anp D. D. AUFENKAMPt 


.Summary—In this paper we indicate briefly how the methods of 

algebraic solution of Markov chains with constant transition 
babilities apply to the analysis of sequential machines. Mealy’s 
del of a sequential machine is assumed. A stochastic matrix is 
ociated with each such machine to provide a starting point for the 
lysis. ‘‘Closed” sets of states are then characterized, for example, 
the appropriate theorems about Markov chains. A technique is 
mulated for reducing the connection matrix of any sequential 
thine to a canonical form. 


TN HE methods of the algebraic solution of Markov 
chains with constant transition probabilities pro- 


~- vide techniques for the analysis of the matrix 


* Manuscript received by the PGEC, March 21, 1958. 

‘+ Boeing Airplane Co., Seattle, Wash. ; 

tt Lockheed Missile Systems Div., Palo Alto, Calif. Formerly at 
:d College, Portland, Ore. 


representation of any sequential machine. Considering 
the formulation established in two previous papers! ? 
these techniques can be used to reduce the connection 
matrix of a sequential machine to a canonical form. 
We associate with the state diagram of a sequential 
machine a matrix, P=[p,;|, called the probability tran- 
sition matrix, as follows. If there is at least one input 
taking the machine from state z to state 7, we designate 
a real number p;;(0<p,; <1) as the transition probability 


1F, E. Hohn, S. Seshu, and D. D. Aufenkamp, “The theory of 
nets,” IRE TRANS. ON ELECTRONIC COMPUTERS, vol. EC-6, pp. 154— 
161; September, 1957. 

2D. D. Aufenkamp and F. E. Hohn, “Analysis of sequential 
machines,” IRE TRANS. ON ELECTRONIC COMPUTERS, vol. EC-6, 
pp. 276-285; December, 1957. 
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from state z to state 7. If there is no input connecting 
4 to j, then p;;=0. In the event of a terminal state, k, 
we assign py, =1 without essentially altering the struc- 
ture of the machine. 

In this paper we impose the restriction that 


q 
Pa = 1 
73 


(where g is the number of states). Under this restriction 
the probability transition matrix P becomes a stochastic 
matrix.*:* We also define an initial probability distribu- 
tion {a,}, where a, is the probability that the machine 
is initially in state k. The stochastic matrix P, together 
with the initial distribution {a,}, completely defines a 
finite Markov chain with constant transition probabili- 
ties. A sequential machine so described will be called a 
Markov machine. 

If a transition from state 7 to state k in m steps can 
occur via the sequence of states j, fi, jo, °°: , Jn, R, 
then the probability that the machine passes through 
this path, if initially in state 7, is Dj;,Dj1j. °° * Dig—ake 
Hence, the sum of the probabilities over all the distinct 
paths between states z and j is the probability of finding 
the machine in state k after a sequence of 7 inputs if 
the machine is initially in state 7. This probability, 
called a higher transition probability, is denoted by 
pij™. We have, for example, 


q 
Pr? = pny Pir = DY Pirpne™ 
rh 


The higher transition probabilities can be calculated by 
taking powers of the probability matrix P. For instance, 
Prtt=[p,,+) ] = PP». If the initial probability of find- 
the machine in state 7 is a;, then the unconditional 
probability of finding the machine in state k after a 
sequence of ” inputs is given by: 


qd 


a = D5 apn”. 


g=1 


The state k can be reached from state 7 if there exists 


some n such that p;,>0. A set of states S will be . 


called closed if no one step transition is possible from 
any state of S to any state outside of S, that is, if pj, =0 
whenever state j isin S and state kis not. We note that 
if state j is in S and if state k is not, then p;, =0 for 
all n. If state 7 isin S, then the sum of P;, extended over 
all those v for which state v is in Sin unity, 2.e., 


Dd iv = 1. 


3 W. Feller, “Probability Theory and Its Applications,” John 
Wiley and Sons, Inc., New York, N. Y., vol. 1, p. 309; 1950. 
; 4 Jt may be that the nature of the sequential machine is such that 
‘ the. branch=weiglits of the state diagram are already transition 
probabilities. As long as the condition 


q 
> ba = 1 
j=1 


holds, the subsequent development is relevant. 
5 Feller, op. cit., p. 309. 


IRE TRANSACTIONS ON ELECTRONIC COMPUTERS 


June 


Thus, in the jth row of P” the elements corresponding 
to states in S add to unity. Therefore, if in the matrix P- 
all rows and columns corresponding to states outside 
of S are deleted, the remaining submatrix is again a 
stochastic matrix. Hence, closed sets of states may be 
studied independently of the remaining states. A closed 
set of states is called irreducible if it does not containa 
subset of states which in itself forms a closed set. a] 

Next, we apply the method of generating functions 4 | 
to determine the number of irreducible sets of states” 
in a sequential machine.® Define for fixed 7 and k: 


P55). pee 


n=l 


Multiplying by sp,; and summing over 7 obtains 


q j 
SDE PriPix(s) = Puu(s) — por, (1) 

a 
where g is the number of states. For each fixed k there 
is a system of g nonhomogeneous linear equations in the 
g unknowns Py(s), Pox(s),---, Pals). The deter- 
minant of the system, D(s), is thus a polynomial in s of 
degree less or equal to g, while the P,,(s) are rational 
functions of s with the common denominator D(s). One 
of the roots of the equation D(s) =0 is s=1 since 


s(X pu) -1=0. 
j=l 


On the other hand, each closed set of states contributes - 
a factor of (s—1) to the determinant, D(s), of a system 
of equations formed as above. This assertion is evident 
from the following observation. Suppose the system of 
equations in the unknowns Px(s) were rearranged 
such that the coefficient matrix had the form 


Ad 102209 
OR = = Oe 


O¢.0.; Ie0 
alle pues 7 


where closed sets of states are indicated by each of the 
square submatrices, 4, B,---, Y, while Z, although a 
square matrix, does not represent a closed set of states. ; 
The a, 6, --+,A, and 0’s are rectangular submatrices, 
and, in particular, the 0’s represent null matrices. 
The determinant of the above matrix is | A| | BI 
qc | Y| eae Since A represents a closed set of states, 
there is at least one factor of (s—1) in | A | . If A repre- — 
sents more than one closed set of states, then there must 
be more than one factor of (s—1) in | A | : Infact, if A 
represented, say, two irreducible sets of states, then re- _ 
writing A in the form j 


5 Let ao, a1, + + + be a sequence of real numbers. If A(s) =ao+-ais i 
+a2s?-+ + + + converges in some interval —sy<s<so, then the func- 
tion A(s) is called the generating function of the sequence {a;}. 
Feller, op. cit., page 212. 
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Ay 0= 10! 3 
A'=|0 Az 0 
La’ ou i 


nere A, and A» represent the irreducible sets, we see 

at (s—1) must be a double factor of | A bi: We reason 
ilarly for | Bl, Sore | Y|. The determinant of Z, 
wever, can not have a factor (s—1) since Z does not 
resent a set of states, 7.e., 


| (Xp) -1 0 


| 
| 


r s=1, where the summation is over those states 
presented by Z. Hence the theorem: 


‘Theorem I: The multiplicity of the root, s=1, in 
.D(s)=0 indicates the number of irreducible sets of 
states in a sequential machine. 


In the event that the equation D(s)=0 does not 
ive any multiple roots the following procedure can be 
ed to determine the p;,“™ from which the structure of 
e machine can be determined. The technique is out- 
ed only; the reader is referred to Feller for its deriva- 
mn.’ We note that since the P,,(s) are rational it fol- 


ws that there exist coefficients, a4, o,, + + +, on™, 
ch that 

| Oe gi, ©) op 

| P= a (2) 
| iS siy So” Sar 

HCTe Sj, So, °° *, Sq are the roots of D(s).o. Now, con- 


ler the two systems of equations 


: es 2 bite (i = 1, 2,-*-, 9), (3) 
: e 

: q 

: Yn = $ Se Vk Prem (m = iy 2, Cen q)- (4) 
| k=1 


re two systems have a common determinant, and 
éy have nontrivial solutions only for values of s for 
iich the determinant vanishes. Observe that the 
ove determinant equals the determinant, D(s), of 
m1), Let si, So, --*-4 S_ be the roots (there are at 
m= 2q).or J(s)=0. For-each.s,, r=1,,2,- -- 5 q, the 
lutions {94, PRS xq” and fy, > Ghar yg} are 
termined to an artibrary multiplicative constant. The 
) in (2) are obtained from the relation 


on 6” = Crt yy, ™, (5) 


iere the constants c, are determined from® 


q 
1 = c >, ay, (6) 
v=1 


7 Feller, op. cit., ch. 16. , 

8 Although the solutions x,” and y,“ are determined only to a 
merical factor, if we replace x; by ax; and the yz by By., 
n ¢, will be changed into c,/a@8 and consequently the oj,” in (5) 
| remain unchanged. 
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As an example, consider the following probability tran- 
sition matrix: 


ag HOO 
1000750 
0} 0.5.04 A) 
1EUQ Os, 0 


The systems of equations indicated by (3) and (4) are 


SX SXq SV1 
t= ee ; Die eee aioe 
SV1 
X2 = SH, yo = ar 
Gy as Soy y3 = 0 
x4 = SH ys = SVs 


Since the multiplicative constant factor is arbitrary, we 
set x4=1, yo=1, and solve the resulting systems to ob- 
tain: 


1 2 1 D 1 
i ie! at Ce ee 5 ae x3 = S$; eee ag ees 
§ s? Ss St § 
and 
2 2 1 
Via oS ine eee ae 
§ S Ky 


The roots of D(s) =0 are s;=1 and s.= —2. Thus 
{aj} Picea. { xj} = {—-4, 1, —2, 1}, 
{ym} = {2,1,0,0}, and {y,@} = {—1, 1, 0, 0}. 


From (6), the c, are obtained: c;=}3, @=%. Calculating 
the 0; from (5) and substituting these values in (2) 
we see that pis™ = pia™ = po3™ = pou™ =0, while the 
remaining p;,™ are different from zero. Therefore states 
1 and 2 form a closed set of states which is evident from 
the state diagram in this example. 

For machines in which more than one closed set of 
states exist the following method of analysis has been 
developed. Given the probability transition matrix P 
=[p,;| of a Markov machine having g states, consider 
the system of homogeneous equations 


k 
t= Lo Pyty k= 
g=1 
This system will have a nontrivial solution if its deter- 
minant is zero. However, we know that a nontrivial 
solution exists; namely; %7=x;) 7 = 1) 2)~ ~:, ¢-itence, 
the following theorem: 
Theorem II: A necessary and sufficient condition that 
a set of states, { i, Oy, ar is}, be a closed set 1s that in 
the system of homogeneous equations, 


q 
AH aS PrjXj, 


j=l 
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there exists a subset of t equations in the t unknowns, 
Niyy Vigy °° * Vige 


Therefore, we now have a method of determining 
closed sets of states within closed sets of states, obtain- 
ing eventually the irreducible sets of states. In applying 
the above procedure the work is simplified if we group 
those equations whose right-hand member consists of 
two terms, etc. Examination of the system for subsets 
of k equations in k unknowns will yield all the closed 
sets of states. 

As an example, consider the probability transition 
matrix: 


GAN Ome O'eLOMEt eNO uLOn eon ONO 
Del One EAOeRO Dee ONO 10 
Te sOe mC RO) On is 00 
Aen Ome 0.50 he ON 0. 07-05-00 
Slice OO WA. G0 Oso 
Oe Om. er Osh O% Odie nO sO 
7a Ome WOL RO! 0752 0810 20026 
Sal tO) 2 Oo OunO OF 0 
OMe Omg 0 06.0 Oe, <0 
fOmnG te Unga Oot 20 200 
Forming the equations 
a 
Xe = DY past; 
j=l 
and rearranging them in accordance with the above 


suggestion, we have 
xX = 0) 
Ka 04 
t7 = ZX2 aie 5X6 
5 = 34%, + 245 


et 2 
xg = 3%, + 3v5 


1 = aX1 + 2X5 ae zX8 


Xe = $41 + 243 + Baz 


CORRECTION 


Joseph O. Campeau, author of “The Synthesis and 
Analysis of Digital Systems by Boolean Matrices,” 
which appeared on pages 231-241 of the December, 
1957 issue of these TRANSACTIONS, has requested that 
the following corrections be made to his paper. 

On page 238, in paragraphs two and four, the term 
2k" should be 2**—1. 

Immediately after (21), the following statement 


IRE TRANSACTIONS ON ELECTRONIC COMPUTERS 


_ published thesis of R. G. Gillespie presented to the Division of 


Jun 


%3 = au ak aX7 a 3X3 
Xo = e%1 =f: Bu3 = EX6 +e 5X9 
io = 3X1 + $4 + ox6 + Gre. 


States 4, 6, 5, 8 and 1 form a closed set of states. De- 
composing this set we see that states 4 and 6 each form 
an irreducible set of states as does the set consisting of 
states 5, 8, and 1. 

Theorems I and II permit the reduction of the prob- 
ability transition matrix of any sequential machine to 
the form 


OO Om sre ® 
i TG ee NZ, 


where A, B, ---, Y are square submatrices represent- 
ing closed sets of states. Z is a square submatrix not 
necessarily representing a closed set of states. The 
a, 8, +--+, are rectangular arrays of entries which are 
all zeros if and only if Z represents a closed set of 
states. Each of the A, B,--- may perhaps be further 
reduced into matrices of similar form 


Al 20) Be ZON AO 
O03 Ba 6 


Wedd hence ee abe? 
a’ Biep SN WHE 


and so on, until we arrive at the irreducible sets of 
states and the resulting canonical decomposition. 

The preceding development serves to indicate very 
briefly the powerful methods for analysis of sequential 
machines provided by the techniques developed in 
studies of Markov chains.?® 


* A more complete study of these applications appears in the un- 


Mathematics and Natural Science, Reed College, in partial fulfill- 
ment of the requirements for the Degree of Bachelor of Arts. 


should be added: where the 4;;’s are the values a par- 
ticular input 4; can take. 


The ranges on 7 and jin (22) should be 
nm Wee en ey 


pean 
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Summary—The use of digital computers in continuous control 
ems is discussed. A comparison is made of the two major types 
heital machines, namely the GP (general purpose) and the DDA 
ital differential analyzer), and characteristic features of each are 
sidered. Certain advantages and limitations of each type are 
cribed, together with their indicated areas of application. Methods 
rocessing input data from analog and digital sensing elements are 
cribed, as well as means of supplying control signals to output 
TOS. 


ADVANTAGES OF DIGITAL COMPUTERS 


\\ NALOG and digital computers, separately or in 
combination, may be used in automatic control 
systems. Where applicable, a digital computer 

ers the following advantages. 


.) As much precision as may be required in any par- 
ticular application. 

:) Improved reliability because of relative freedom 
from problems encountered in electromechanical 
or electronic analog computers, namely precision 
of components, wear, and stability. 

') Relative insensitivity to mechanical or electrical 
“noise,” which limits the complexity of problems 
soluble on an analog computer. 

-) Facilities for storage of large quantities of data 
that may be required in complex systems. 

)) The ability to perform logical as well as linear or 

nonlinear mathematical operations. This is a com- 

mon requirement in applications calling for a 

master control to supervise an entire system. 

») Flexibility derived from its logical capacity and 

- the fact that its operations are controlled by a 

program rather than specially designed or inter- 

~ connected physical components. 


“ REVIEW OF CHARACTERISTICS OF PRESENT GP 
AND DDA MACHINES 


[he flow of information through a digital computer 
4 control system is shown in Fig. 1. The chief differ- 
ses in the two types of machines are related to the 
t that a DDA has a semifixed program, whereas that 
the GP is not fixed at all. Also, the arithmetic and 
itrol circuitry in a DDA are extremely simple. 

[The GP computer solves problems in a manner simi- 
to that used by an operator and desk calculator. The 
sic Operation is to obtain two pieces of data from 
rage, perform a required operation on them, and 
nsfer the result back into storage, or use it as an 


* Presented at Symposium on Computers in Simulation, Data 
luction, and Control, New York, N. Y., March 21, 1957. In Mr. 
un’s absence the paper was read by Dr. Harold K. Skramstad, 
tional Bureau of Standards. 

+ Genesys Corp. (subsidiary of Chance-Vought Aircraft Co., 
), Los Angeles, Calif. 
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Fig. 1—Information flow through a digital computer 
in a control system. 


output command. Whole numbers are transferred be- 
tween the operational units, and the machine may ac- 
cordingly be referred to as being of the absolute type. 

The number of components external to the main 
storage is usually large, since control of the arithmetic 
unit must provide adequate sequencing and program- 
ming logic. Also, switching circuitry must be provided 
to connect the arithmetic unit to any part of the 
storage. 

In a GP machine, time is consumed for search opera- 
tions, 7.e., obtaining access to data and instructions in 
storage. Also, in many arithmetic operations, a con- 
siderable amount of redundant computing may be done 
if the values of the two items of information change very 
little from one computing cycle to the next. This is par- 
tially offset by the fact that ina GP machine, computing 
may be done in large increments, and therefore, the 
sampling rate can be reduced to the minimum value 
consistent with stability of a control system. 

To solve a specific problem, not only initial values, 
but a detailed sequence of instructions for solving the 
problem, 7.e., a program, must be entered into the main 
storage unit. 

The DDA! operates through an arrangement of 
switching and storage elements whose behavior is anal- 
ogous to that of a system governed by the equations 
being solved. The information transferred between op- 
erational units in a DDA consists of single increments 
of the variables. Therefore, this type of machine is often 
referred to as an “incremental analyzer.” 

The operational steps required in a differential ana- 
lyzer are permanently wired into the logic network. Asa 
result, paired words comprising the “Y” and “R” regis- 
ters of a digital integrator are read from the internal 
store. Because these words are permanently paired, and 
each pair is always operated upon in sequence, no time 
is consumed for “search” operations. 


1E. L. Braun, “Design features of current digital differential ana- 
lyzers,” 1954 IRE NATIONAL CONVENTION RECORD, pt. 4, pp. 87-97. 
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The arithmetic circuitry is simple for two reasons. 
First, the only arithmetic operation required is sum- 
mation. Secondly, the summation is performed serially, 
a bit at a time. After each binary place of an integrator 
is operated upon, the result of that operation is sent to 
the magnetic drum for storage. 

The control circuitry of the DDA has built into it a 
semifixed program. Asa result, the solution of a problem 
requires only the following operations to be performed: 
1) determination of how operational units are to be in- 
terconnected; this information is inserted into the ma- 
chine by a simple code, 2) scaling of the problem, and 
3) insertion of initial conditions into the drum storage. 

It should be emphasized that the simplicity of the 
arithmetic and control circuitry ina DDA is not achieved 
by compromising performance, but is rather the result of 
the nature of the machine. The DDA is designed to 
simulate continuous systems by integration of a set of 
continuous variables starting from a set of prescribed 
initial conditions. The integration process is approxi- 
mated to any desired degree of precision by means of a 
simple addition and storage process. All other mathe- 
matical operations that may be required, e.g., multipli- 
cation, division, function generation, differentiation, 
can be derived for continuous functions, from this 
process. 

For reference, the characteristics of the DDA which 
tend to simplify its control circuitry are listed below. 

1) It has a semifixed program. 

2) The only store to which access is required by each 
operational unit is the so-called dz store, which 
holds the current outputs of all the operational 
units. Its size is usually of the order of only a few 
hundred bits, depending basically on the number 
of operational units, and whether a binary or ter- 
nary representation of positive and negative values 
is used. 

3) Only one bit at a time is read, recorded, or oper- 
ated upon and always ina predetermined manner. 

4) The arithmetic unit is comprised of simple adder 
circuits only. 

The mode of operation of a GP machine requires that 
the control provide access to a main storage unit, hold- 
ing several thousand words, and also that transmission 
links can be established between the main storage unit, 
the arithmetic unit, sequencing control, and input-out- 
put equipment, at the option of the programmer. Be- 
cause of the great diversity of possible programs, a 
number of built-in transfer instructions must be made 
available to facilitate the task of programming. 

If one were to generalize on the relative programming 
flexibility of GP and DDA machines, some features of 
the GP machine that could be cited are: its facility for 
using subroutines, the availability of more logical 
operations, better adaptability for the addition of new 
instructions, and the ease of programming access to any 
part of storage. The importance of programming flexi- 
bility depends, in part, on whether improved programs, 
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and data on the control system’s parameters are ob- 
tained, and, if so, when. 

It is possible to build a GP computer with only one 
built-in arithmetic operation, namely subtraction, and 
to derive all other operations required during a computa- 
tion by means of programming. However, this is not 
very useful, except possibly for purposes of teaching or 
amusement, because of the extremely cumbersome pro- 
gramming procedure that this entails. Also, it means an 
extremely inefficient use of storage capacity and a great 
reduction in speed of computation, compared to a GP 
machine having built-in arithmetic units for addition, — 
multiplication, and division, as well as subtraction. 


AREAS OF APPLICATION OF GP AND DDA MacuINES 


The principal considerations that influence the design 
of a computer for continuous control system applica- 
tions are given below. 

The rate at which solutions have to be obtained. The re- 
sponse time of a digital computer in a control system is a 
function not only of the computer’s speed, but also of 
the complexity of the system in which it is incorporated. 
In other words, since a digital computer functions by 
performing the steps of a computational routine serially 
the more steps there are, the longer it takes to perform | 
a single cycle of computation. 

The specified accuracy of the solutions. The subject of 
error analysis is outside the scope of this paper. How- 
ever, a few comments are in order. Errors may arise 
from two sources (exclusive of equipment failures), 
those inherent in the nature of the computer and those 
peculiar to a specific problem. 

A source of error common to GP, DDA, and other 
machines is that due to round-off. The DDA is also 
subject to a special truncation error due to remainders 
in the R registers. 

Another source of error is that of phase shift. This is’ 
caused by time lags in the transmission of data between 
parts of the system. 

To minimize errors, the forms of the equations chosen’ 
must be such that they adequately represent the sys- 


acteristics of the computer. 

Ina GP machine, the over-all error is governed largely 
by the numerical method and program selected. In the 
DDA, the general numerical method is relatively fixed, 
but there is still considerable choice both in the func- 
tions chosen and in the way they are generated. A char-. 
acteristic of the DDA not shared by the GP machine is 
that an increase in the scale proportionately increases 
the time required forasolution, 7. e., product of the band-- 
pass and reciprocal of the precision is a constant. 

The form(s) in which information is (are) to be trans- 
mitted to various parts of the system. The principal differ- 
ences between the GP and DDA machines in the inter- 
nal transmission of information have been mentioned. 
Transmission from sensing elements and to actuators is 
discussed in the sections on input-output equipment. 
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Fig. 3—Function generation requiring a GP machine. 


The exact areas of application of DDA’s and GP ma- 
ines within the control system spectrum have not 
en established to date. There are a number of reasons 
this. First of all, there is the fact that this is still a 
w field. Secondly, the characteristics of many of the 
ntrol systems in which a digital computer may be in- 
3 have not been adequately specified. Third, there 
-he important point that a decision in respect to which 
ne is better suited for a particular application cannot 
made on the basis of the general nature of each ma- 
‘ne alone. This decision is influenced to a large extent 
the types of components available for fabricating 
sh machine. For example, before the advent of the 
‘gnetic drum store, the utility of the DDA was some- 
at questionable. Now, with new, compact, high- 
sed, reliable components becoming available, it may 
desirable to use these to replace the drum for applica- 
ns calling for considerably higher speeds. However, if 
1er components are used, e.g., if the DDA’s opera- 
nal registers are fabricated from magnetic cores and 
nsistors, its area of best application and the econom- 
of its competitive position with GP or analog ma- 
nes is shifted with respect toa DDA using a magnetic 
4m. Accordingly, the statements made in the follow- 
- discussion concerning the relative advantages of 
JA and GP machines must be tempered by the under- 
ng assumption that a DDA using a magnetic drum 
re is being considered. 
The DDA is the digital equivalent of an analog type 
simulator, and suited for real time control systems 
-hin bounds set by speed requirements. If the registers 
a DDA’s operational units are stored on a magnetic 
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Fig. 4—Establishment of a computed initial condition in a DDA. 


drum, the computing iteration rate is limited by the 
angular velocity of the drum. The pulse repetition rate 
is a function of the angular velocity and diameter of the 
drum, and the pulse packing density. 

The reason why some computations can best be solved 
by a DDA and others by a GP machine can be seen by 
reference to Figs. 2 and 3. Fig. 2 shows successive values 
of a computed function at different sampling times. This 
function is such that the maximum change in its value 
between successive sampling times is always less than 
some value K. The number, K, represents the maximum 
number of increments that can be generated in a given 
time interval (S;.:—S;) in a DDA. Therefore, if the 
value of the function to be computed does not change 
by more than K, in the time interval (S;.1—S;), the use 
of a DDA is indicated because of its relatively greater 
efficiency in the computation of continuous functions. 
In Fig. 3, the computed function is such that the maxi- 
mum change in its value between successive sampling 
times can be greater than K. Therefore, a DDA cannot 
be used. 

Two other areas where a GP machine is indicated are: 
1) where a large table of stored values must be con- 
sulted, and 2) where the interval of time allowed to es- 
tablish correct initial conditions is too short for a DDA 
to do so. 

The situation referred to by 2) is illustrated graphi- 
cally in Fig. 4. At some time, é,, an “alert” time, a signal 
is made available to the computer indicating that a new 
set of initial conditions is to be inserted into a computa- 
tional loop. If the initial values of some of the variables 
must, in turn, be derived from a computational process, 
the following limitation must be considered. The DDA 
cannot generate more than some maximum number, say 
K, of increments during the interval At. Therefore, if 
during the allowed set up interval, At, the computed 
function must change by a number of increments greater 
than K, special means must be used to establish the cor- 
rect conditions at time, fo. 


An Example of “Real-Time” Computation 


An important application of digital computers is to 
provide guidance information automatically aboard 
high speed airborne vehicles. To achieve this purpose, 
the vehicle utilizes measuring and detecting instru- 
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ments which monitor its position in space and the en- 
vironment, and controls which enable it to change its 
position in accordance with current requirements. 

The computer receives its inputs from such sources 
as radars, gyro compasses, air speed indicators, roll and 
pitch indicators, barometric and radar altimeters, etc. 
It transforms this data to produce outputs in the form of 
control signals to the auto pilot, automatic tracking 
signals to the radars, display signals, etc. Derived quan- 
tities include present position in terms of latitude and 
longitude, distance to destination, etc. 

Essentially, the guidance problem consists of resolv- 
ing the aircraft’s motion, as sensed by its instruments, 
along the axes of some chosen coordinate system. The 
computation problem is, therefore, basically that of 
solving trigonometric and algebraic equations, plus the 
integration of certain variables. The mathematical 
operations that will enter into the computation may be 
classified as follows: 

Addition and subtraction. 

Multiplication and division. 

Generation of trigonometric functions (which may be 

restricted to generation of sines, cosines). 

Generation of inverse trigonometric functions. 

Integration. 

Simple function generation: exponentials, absolute 

values, etc. 

Analog-to-digital conversion. 

Digital-to-analog conversion. 


The computation may be performed either by an ab- 
solute or incremental computer. However, the nature 
of the problem particularly lends itself to an incremental 
technique, since the variables are, in general, continu- 
ous and have limited rates of change. 


GENERAL CHARACTERISTICS OF THE INPUT SYSTEM 


A control computer’s function is to supply control 
signals to output servos. These commands are based on 
information received from sensing elements. Ideally, the 
information supplied by the sensors to the computer, 


and by the computer to the servos, should be in a form - 


directly usable without further conversion. Generally 
speaking, this will not be the case. 

An input sampling system for a digital computer in a 
control system may comprise, in general, the following 
elements. 

1) Sensing devices (analog and/or digital). 

2) Data converters (analog-to-analog, 

digital, digital-to-digital). 

3) Buffer storage. 


analog-to- 


The buffer performs the storage and switching func- 
tions required to transmit various encoded inputs to the 
computer in sequence. A timed sequence of orders to 
control these operations is supplied by the control unit 
of the computer. 

These elements perform the following types of opera- 
tions in synchronism with the computer: measuring, 
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Fig. 5—An analog input system. 


encoding, and time-division multiplexing. 

Fig. 5 illustrates a system for sampling a number of — 
analog inputs in sequence. Two important means of — 
equipment simplification are illustrated: time-sharing, 
and analog-to-analog converters. 

To allow time-sharing of encoders, the outputs of © 
instruments that are of the same form are combined by 
the use of sequencing switches controlled by signals 
from the control unit. For a large system, the cost of the. 
sequencing switches and associated control equipments 
is outweighed by the simplification in encoding and 
buffering equipment. 

By the use of analog-to-analog converters, the outputs 
of various instruments may be transformed to forms 
more amenable to encoding, and the time-sharing of 
encoding and buffering equipment increased. A simple 
example of an analog-to-analog converter is a potentiom- 
eter producing a voltage proportional to the position 
of a shaft to which it is attached. 


Analog-to-Analog Converters 


The number of possible representations of data 
which input instruments can produce is usually limited — 
to four: 1) shaft position, 2) voltage, 3) time, and 4) 
frequency. f 

It is relatively simple to convert a time interval to 
binary representation. Therefore, one way of standard- 
izing the input problem would be to convert shaft posi- — 
tion and voltage to a time interval. There are two prin- _ 
cipal methods of voltage-encoding into the time domain. | 
One method derives timing information by means of a 
generated sweep voltage and the other by means of a 
phase shifter. Fig. 6 shows block diagrams of each 
method. ¢ 

In sweep-timing a direct voltage varies linearly from 
zero to a maximum expected value during each sampling 
period. When the unknown voltage equals the sweep 
voltage, the coincidence is detected. The time interval 
between the start of the ramp waveform, and the in- 
stant when it equals the unknown voltage is propor- 
tional to the value of the voltage. 

In phase shift timing, the analog quantity is repre- 
sented by the phase shift of a variable voltage with re- 
spect to a fixed reference voltage. The difference i 
phase is converted to a time interval by applying each 
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Fig. 6—Voltage-to-time encoders. (a) Generated sweep timing. 
(b) Phase shift timing. 


Itage to a comparator. The interval is defined by the 
mes at which the reference and variable voltages, re- 
ectively, pass through 0°. This method is advanta- 
ous if the original analog quantity is easily converted 
a phase shift. It is used mainly to convert rotating 
aft data, since phase shift data can be obtained by 
sans of synchros. 

The nonlinearity of the ramp voltage generator or the 
ase shifter and the uncertainty in the coincidence 
cuits limit the accuracy of the voltage-to-time do- 
Lin encoders. 

Fig. 7 illustrates how a shaft position may be con- 
rted to a time interval by means of a resolver and a 
ase-shift type of voltage-to-time encoder. 


valog-to-D igital Converters 


‘Conversion of a Tvme-Interval to a Number: Under dis- 
ssion of analog-analog converters, means were shown 
- converting a voltage to a time interval. Fig. 8 shows 
w a binary representation is obtained by counting the 
amber of clock pulses that appear in this interval. 
(e required pulse repetition rate of the clock source 
pends on the specified precision and the input sam- 
mg rate. 

Direct Conversion of a Shaft Position to a Binary Code: 
*vices for converting shaft position to a binary code 
Ly be either of the incremental or absolute type. The 
tremental type produces a pulse whenever the shaft 
cates an amount specified as a unit increment. In the 
solute type, the position of the shaft is monitored at 
times. The latter method provides the advantage of 
t being susceptible to an accumulation of error re- 
‘ting from lost pulses or “noise.” 

‘\ widely used shaft position converter utilizes a 
led disk coupled to the shaft to be monitored. If 
tre are restrictions on the size of the disk, two or 
sre disks may be driven in some suitably coupled ar- 
agement, with the less significant binary places being 
\d from one disk and more significant places from the 
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Fig. 8—Conversion of a time interval to a binary number. 
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Fig. 9—A coded-disk type of angle encoder. 


other. These disks may be read either optically or by 
brush contacts. They may be scanned either in a serial 
or parallel manner. To avoid ambiguities in a read-out, 
a reflected binary type of code or a special arrangement 
of brushes may be used. A schematic of a coded disk 
type of angle encoder is shown in Fig. 9. 

There are a number of advantages to the coded-disk 

type of angle encoder: 

1) Conversion of data from analog to digital form 
can be made without intermediate analog conver- 
sion steps. 

2) The accuracy and repeatability of data using this 
system is relatively independent of variations in 
supply voltages, frequency, and temperature. 

3) Maximum time-sharing is possible, and switching 
is performed on digital rather than analog signals. 

4) The conversion time is compatible with the rate 
at which input data is sampled. 


The Forms of Digital Inputs 


Digital data can be obtained directly and/or from 

data converters in the following forms. 

1) A whole number proportional to the value of an 
input variable, either serially, in the form of 
weighted pulse trains, or in parallel. 

2) A number of pulses porportional to a change in 
the input. 

3) A pulse rate proportional to the sampled value of 
the input. 


Absolute measure has the advantage of not being sus- 
ceptible to a permanent offset caused by a lost incre- 
ment. Incremental measure may be satisfactory if it is 
inside a broader closed loop with zero steady-state error. 
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Fig. 10—Conversion of a frequency to a binary number. 


An example of a system where incremental measure 
may be satisfactory is one in which new absolute refer- 
ences are inserted at frequent intervals. 

There is an important difference in the way input 
quantities are processed by a GP ora DDA. Ina DDA, 
transmission of information between operational com- 
puting elements is by means of unitary weighted pulse 
trains. A special register, termed a servoregister, is used 
to generate a unitary weighted pulse train from the ab- 
solute binary input data. In the input problem, as well 
as within the computer, the characteristic limitation of 
the DDA, namely not being able to compute at a rate 
greater than K, presents itself. If the output of a servo- 
register is to be used as an independent variable input 
somewhere within the computer, then increments can 
be accepted from the servoregister only at a rate <K. 


Digital-to-Digital Encoders 


Conversion of a Frequency to a Number: Variable pulse 
rate data can be obtained from sources like pressure 
meters which produce a frequency proportional to pres- 
sure, or Doppler radars which emit a beat frequency 
proportional to the rate of change of range. This data 
can be converted to a binary number as in Fig. 10. 


OUTPUT SIGNALS TO CONTROL ELEMENTS 


Data in absolute form, 7.e., as numbers, can be con- 
verted to,control signals in the form of a voltage input 
to a servoamplifier. This requires a holding register and 
a voltage encoder. There is essentially no difference in 
the way this is done in a GP machine or a DDA. 

Incremental signals, too, can be obtained from either 
a GP machine or a DDA. In the former case, increments 
are obtained simply by computing differences between 
successive control signals. 

A method of converting a number to a shaft position, 
frequently employed formerly, is shown in Fig. 11. In 
this scheme, the positioning error is produced by a 
voltage comparison. 

A system now used extensively is shown in Fig. 12. 
Here, the positioning error is produced by a digital com- 
parison. The signal to the servoamplifier can be of the 
ON-OFF type or made proportional to the error over 
some specified range by use of an error register. The 
error register is set to the semiproportional digital dif- 
ference of the monitored and desired shaft position. The 
conversion from the semiproportional digital error, con- 
tained in the error register, to error voltage is affected 
by using a binary weighted network. The output of this 
network is sent to a holding or desampling device which 
presents the output data to the controlled element as 
continuous, 7.e., smoothed data. The simplest method is 
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Frc. 12—Conversion of a number to a shaft position by a 
digital comparison process. 


to clamp the last sample until a new one is computed. 

A relatively simple type of servodrive used exten- 
sively in aircraft control systems is the de pulsing, or | 
notched, motor. It produces an increment of rotation 
for each pulse input. 

A voltage proportional to a pulse rate can be ob- 
tained simply by causing each pulse to add a fixed 
amount of voltage to a capacitative storage circuit. 


CONCLUDING REMARKS 


It is difficult to compare different classes of machines, 
especially when so many variations of design exist and so 
many more are possible within each class. Many items 
must be considered, e.g., the type and complexity of the 
application, precision and speed requirements, choice of 
components. Many items which enter into an integrated 
system design are too detailed and variable to be treated 
here. However, it is not unusual for these details to viti- 


-ate arguments, based on general considerations alone, in 


favor of the superiority of one type of machine. For 
example, consider the problem of reliability. Often, — 
reliability must be obtained at the cost of other factors. 
There are many devices that may be used in an attempt | 
to improve reliability; proper selection and rating of 
components, duplication of equipment, duplication of 
programs, programmed checks, built-in checks, etc. 
Each of these affects the design of a particular machine, — 
and may affect one design to a greater degree than an- 
other. 

In short, conclusions about the aptness of a computer 
for a specific application cannot be reached without ref- 
erence to specific computer embodiments, and the whole 
system of which the computer is a part. Computer “A” 
is better than Computer “B” only if the system includ- 
ing “A” more nearly satisfies specified requirements. 
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Computers in Process Industry Control" 
WILLIAM F. GUNNINGT 


Summary—Successful inclusion of a digital computer in the real- 
ne control loop of an industrial fluid process plant is dependent on 
practical steps: 1) the formulation of adequate mathematical 
odels, 2) the development of satisfactory data acquisition tech- 
ques, and 3) the proof of sufficient reliability. Examples are taken 
9m the petroleum industry. 


INTRODUCTION 


| HIS paper considers some past, present and future 
) applications of computers in simulation, data re- 
; duction, and control in the process industry. The 
rocess industry, by definition, is largely made up of 
ympanies whose raw materials, intermediates, and 
roducts are primarily fluids. While most of the exam- 
les used in this paper are taken from the petroleum in- 
ustry, the basic problems are found in plants producing 
ich diversified items as paint, catsup, plastics, and 
ottled gas. 

Over the past 20 to 30 years, the process industry has 
2en a leader in the development and application of 
-rvomechanisms for what might be called “minor 
op” industrial control. Many of the automatic regu- 
‘tors now in use contain ingenious, reliable, and inex- 
ensive analog devices for addition, differentiation, and 
itegration. Although they are seldom glamorized by 
ye term “computer,” it is difficult to distinguish be- 
veen the information-processing they perform and that 
ene by components of present day analog computers. 
The ease of handling afforded by fluids, in combina- 
on with this simplified version of automatic control, 
as resulted in an industry that needs comparatively 
‘w workers per dollar of product produced. At first 
ance, it would hardly seem justifiable, therefore, to 
Hd a hundred thousand dollar digital computer to 
ach an installation. It most certainly cannot be justified 
n the basis of saving the money now spent on manual 
yntrol operations, as is the case in less automized ap- 
lications such as accounting and engineering design 
alculation. 

Minor-loop process control has its limitations, how- 
yer, especially in relation to process efficiency, called 
yield,” and in quality control. For example, present 
ay electronic regulators can hold process temperature 
ariations to within a few hundredths of a degree, yet, 
yen with such tight control on such a primary variable 
5 temperature, it is often not possible to produce con- 
stently the exact product required. The complex fac- 
ors leading to truly effective control generally cannot 
= measured directly, but must be determined by com- 
ining information from several primary measuring 


'(* Presented at Symposium on Computers in Simulation, Data 
eduction, and Control, New York, N. Y., March 21, 1957. 
+ Beckman Instruments Co., Fullerton, Calif. 


instruments. And this, then, is the justification for con- 
sidering the use of major computing elements in process 
control. 

The stored program digital computer was proposed 
as the logical candidate for this role at least five years 
ago. While it is true that both analog and digital com- 
puters are available today with speeds capable of real- 
time control, and while several significant steps have 
been made toward their application in such control, 
much additional work remains. The problems may be 
grouped roughly as follows: 


1) The formulation and verification of satisfactory 
mathematical models. 


2) Development of adequate data acquisition tech- 
niques and equipments. 

3) The accumulation of convincing proof that the ad- 
vanced computer and control equipment required 
is sufficiently reliable for process control. 


THE MATHEMATICAL MODEL IN PLANT 
DESIGN AND OPERATION 


The chemical and petro-chemical industries are pres- 
ently in the early stages of using both analog and digital 
computers as aids in 1) the design of new plants, and 2) 
in the operation of existing plants. 

The basic equations of mass transfer, fluid dynamics, 
and chemical reaction rates are considered to be too 
complex to be used directly in the design—by classical 
methods—of new plants or processes. Equations take 
the form of nonlinear differential equations and partial 
differential equations that also include arbitrary func- 
tional relationships based on empirical data. As a re- 
sult, the design of new plants has been based largely on 
extensions of existing working designs. Because of this 
empirical approach, designers have included substantial 
safety factors. Each new design represents only a small 
deviation from its predecessors. 

On the other hand, the airframe industry has been 
virtually forced to pioneer in the use of computers in 
design. Since a contrivance built with the generous 
safety factors used in other industries would never get 
off the ground, these people have had to use comput- 
ers to aid design, in order to have product at all. The 
process industry does not have the same compelling 
impetus to use computer technique but, nonetheless, 
has some important economic reasons for being inter- 
ested. For example, if company A produces a plant 
which is more efficient in design, or produces a better 
product, then company B, addicted to the status quo 
and avoiding the risks of trying new techniques, may 
have a product, but no business. 
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Automatic control in the process industry necessar- 
ily differs from other forms of “automation.” ‘The proc- 
ess industry has an advantage over “Detroit type 
automation” in that its materials are much easier to 
handle and control. On the other hand, it is much more 
difficult to define, in exact mathematical terms, the 
operations which are to be performed. Defining the 
problem of metal cutting is a matter of specifying where 
to place the work with respect to the cutting tool, and 
the rate at which the process will be carried out. The 
inter-actions in the production process in a chemical or 
petroleum plant are very complex, involve a consider- 
able amount of cross-coupling and nonlinearity, and 
can be a mathematical trap for the casual promoter of 
automation. 

Before we can hope to mechanize the decisions in- 
volved in controlling a refinery, we must be able to 
formulate, in mathematical-logical terms, the exact 
response that is to be made for each operating contin- 
gency. What is needed, then, is a reasonable mathe- 
matical model of the process to be automatized. 

The level of the mathematical model can be quite 
varied. The following 5 classes that apply in the pe- 
troleum industry are listed in increasing order of gen- 
erality or scope. 


1) A model can be made of an zudividual component 
in a production process such as a distillation tower or a 
heat exchanger. An example will be given later. 

2) A model can be formulated of a presently operat- 
ing process, including its operating personnel. This 
model may avoid the formulation of detailed models of 
the first class. Such models are already being used to 
schedule efficient production operations in existing re- 
fineries. The model takes into account such factors as 
the storage facilities available, raw materials available, 
and various other economic factors, as well as descrip- 
tions of the plant yield characteristics. Manipulation of 
the model results in efficient operating instructions for 
use by the plant operators in setting automatic produc- 
tion control devices. 

3) A model can include a complete refinery involving 
several plants or processes. Since the products of one 
process can be used as an input or as a fuel at another 
part of a refinery, the scheduling can be a very compli- 
cated affair. Texaco has used an IBM 705 to manipu- 
late such a model to find the most economical and effi- 
cient balance in assignment and selection of various 
crudes and end products. 

4) A more general mathematical model can include 
company-wide operations. This model would include 
such factors as: a) the scheduling of crude production, 
including a model of the company’s oil reserves, b) 
crude transportation and shipping schedules, c) the 
refinery program as it relates to the division of assign- 
ments among several refineries, d) the assignment of 
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storage and distribution facilities for the finished prod- 


ucts, and e) market information based on past sales. — 


The Standard Oil Company of Ohio is reported to be~ 


simultaneously simulating three refinery operations 
on an IBM 705 Computer.! 

5) An ultimate extension of the scope of the model 
includes an entire industry. In addition to the factors 
mentioned before, the techniques of the theory of games 


scelinaaiemn ce 


will no doubt be used to explore the effects of various — 
marketing and production strategies on an industry- 


wide basis, if, indeed, this is not already being done. 
(It might be strategically important to be quiet about 
such activities.) 


As an example of a model of the first class, let us con- 
sider one of the basic tools used in chemical engineering, 
the multiple tray distillation tower. This device is very 
much like a many layered coffee percolator. The input 
material is fed into the center usually as a hot gas. The 
light components bubble up through the trays and 
eventually come out the top asa gas, while the heavier, 
less volatile components condense and drip down to- 


ie 


ward the bottom. Intermediate products may be ex- 


tracted from intermediate trays. There are many other 
variations. 

Many people have reported the results of using auto- 
matic digital computers to aid in the design of distilla- 
tion towers. The first efforts appear to have used an 


IBM 602. A recent paper describes a typical calcula-— 


tion using a 650 and involving a 40 plate tower handling 
15 fluid components.? In the latter, a set of tempera- 
tures and concentrations was determined in 10 to 15 
minutes of computing time. 

A mathematical model useful in designing a tower as 
such is not, however, necessarily best for investigating 


control operations. As has been true of missile guidance | 


and control problems, analog simulation has been more 
helpful than digital simulation. Williams’ reports a case 
which required the fairly large array of 84 amplifiers and 


16 servos to handle a simplified 5-section, lumped model — 


of such a tower. The liquid-vapor balance equations in 


-each tray are nonlinear. The equations used in this sim- 


ulation were simplified (although still nonlinear) in order 
to fit on the analog equipment available. Further, it was 
only possible to consider a two-component stream, 
whereas much more complex streams are typical. The 
analog computer approach, however, allows thousands 
of experiments to be made to explore the system behavior 
with assumed characteristics of various control schemes, 
thus approaching the techniques developed by the air- 
craft and guided missile control people. It is clear that 
analog simulation methods can operate in real-time and 


1 Business Week, April 7, 1956. 


? Bonner, “Solution of Multicomponent Distillation Problems on 
a Stored Program Computer.” 
3 Indus. and Eng. Chem., vol. 48, p. 1008. 
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yw the actual physical controllers proposed for a 
tem to be connected into the simulated model to 
ify the part of the mathematical model used to de- 
ibe their action. 

Jne should always worry about the validity of a 
thematical model, particularly if it has been neces- 
y to make simplifying assumptions. Various methods 
ve been developed to take experimental data from 
its in production. The usual technique involves su- 
imposing a small sinusoidal variation as a “dis- 
bance” to a working piece of equipment and measur- 
“its frequency response. The range of frequencies in 
ypical experiment of this type is from about 1/10 to 
cycles per minute. Time constants representing cor- 
‘s of the frequency response are in the range of 1 to 4 
nutes. It is interesting to note that the instrumenta- 
n used in such a test is sensitive enough to get good 
n and phase data for temperature variation ampli- 
les down to +0.015°F.4 With this sensitivity, it is 
ssible to conduct a frequency response experiment 
th such minute variations that the normal ink trace 
a production recorder-controller chart will not move 
-ceptibly. 

Measurements of gain and frequency response for 
‘ification of mathematical models involving flow 
‘ulation are plagued by noise or interference due to 
w turbulence and other acoustic disturbances. Satis- 
tory reduction of such data has, however, been 
uieved through the use of auto- and cross-correlation 
hhniques.® 

4 large amount of work still has to be done in the es- 
vlishment of suitable mathematical models. Several 
portant devices used in chemical engineering have 
parently not been adequately described by models. 
eked towers loaded with china brick, coated with 
‘ious substances, is one example. A boiling heat ex- 
inger involving both fluid and vapor phase remains to 
‘adequately described. 

Che collection of information by automatic data log- 
:s (about which more is said in the following section), 
i the conversion of this information into a form ac- 
stable to digital computers will facilitate the forma- 
n and verification of accurate mathematical models 
-arger sections of the process by exploiting statistical 
hhniques. A recent use of regression coefficient anal- 
5 on an established process (not in petroleum) has 
closed that some 60 separate independent variables 
-st be measured in order to predict the quality of a 
zle product with reasonable confidence. Cross-cou- 
ia and nonlinear relationships need not be avoided in 
mulating mathematical models, as these may be 
nipulated by a computer to allow the behavior to 
Sompared with actual results measured with the data 
ging equipment. 


A. R. Aikman, J. Instr. Soc. Amer., vol. 3, no. 10, p. 412. 
Reswich, Control Eng., p. 50; June, 1955. 
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Data ACQUISITION 


Most process industry plants are operated in a way 
that requires the logging of hundreds of variables. In 
many instances, the logs are prepared manually. Thus, 
log taking ordinarily represents an analog-to-digital 
conversion on the part of the man who writes down the 
numbers on the log sheet. The results are usually con- 
verted to units more meaningful to operating engineers 
and accountants. There is a growing trend in the proc- 
ess industry to use computer techniques, both analog 
and digital, to mechanize this operation. These auto- 
matic data loggers, as they are called, are usually inde- 
pendent of the recorders and controllers concerned with 
regulating the process. 

Automatic data logging is made difficult by the va- 
riety of measuring instruments used, and the nature of 
the data that must be collected. To be economically 
feasible, data logging equipment must be multiplexed, 
or time shared, by a large number of different kinds of 
information sources. Since these sources vary widely in 
their sensitivity, or electrical response to a change in 
the measured phenomenon, a dynamic range of 100,000 
to 1 is often required in the multiplexed equipment. For 
example, thermocouple temperature signals are an im- 
portant class of process data, requiring high sensitivity. 
Ten microvolts per digit and sometimes even 1 micro- 
volt per digit are typical values. The analog-to-digital 
conversion device must be flexible enough to handle the 
entire range. 

Another problem is that signals collected from vari- 
ous parts of a plant that may be as large as a square 
mile in area cannot be counted upon to be referenced to 
a common ground potential. It is typical to find a volt 
or more of difference in “ground potential” between two 
different remote signal sources in such a situation. This 
means that, for consistent accuracy, the input signals 
must be measured by balanced or floating instruments 
independent of the “ground.” Neither of the two input 
terminals for a channel may be common to any of the 
input terminals on any other channel. 

Finally, it is important that data loggers have the 
ability to do some computing. The analog-to-digital 
conversion is of little value unless the digital printed 
log entries presented are in familiar and convenient 
units. This requires at least the following sorts of mathe- 
matical operation, all now available on present-day 
equipment: 


1) Zero offset and scaling. This is the operation 
Y=A+BxX. 

2) Nonlinear function generation. Operations 1 and 2 
suffice for converting thermocouple voltage measure- 
ments, for example, to temperature in degrees Fahren- 
heidt, Centigrade, or Rankin. 

3) Square Root. This operation is necessary when 
flow is being measured as the square root of the pres- 
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sure drop measured across an orifice plate in a pipeline. 

4) Multiplication and Division. To obtain consistent 
measurements of gaseous flow, it is necessary to “cor- 
rect” for variations in temperature, pressure, and den- 
sity. 

5) Integration. This operation is essential if flow rate 
measurements are to be used in making a mass balance 
on flows in and out of a process or to calculate the 
amount of material moving into storage. 


Integration is the most involved of these operations. 
Typically, many different integrals or totals must be 
maintained and regularly reported. Ordinarily, one 
hour totals as well as 24 hour totals must be calculated. 
Electrical analog methods have been used in data log- 
gers. A pneumatic totalizer is available which combines 
the operations of integration and square root. Coinci- 
dent current, magnetic-core storage is also used to 
keep running totals stored without drift. 

An important development in process industry con- 
trol is the concept of “stream analysis control,” that is, 
the measurement of the actual chemical composition of 
the stream as it flows through the process. (The term 
“end point” analysis has also been used.) 

In the past, process stream analysis has been con- 
ducted primarily on a sample basis by laboratory tech- 
niques, using wet chemistry or various types of spectro- 
metric equipment. Because of the time lag, ranging from 
minutes to hours, stream analysis as such has had vir- 
tually no part to play in direct process control. This in 
turn has forced plant operators to run their processes by 
such “inferential” variables as temperature, pressure, 
or pH. (The inference is that if these are within bounds, 
the product stream is, too.) And this has held up the de- 
velopment of automatic control since the relationship 
between these inferential variables and the property 
one really wishes to control are complex and often ob- 
scure. 

Certain laboratory-type analytical techniques are 
inherently fast enough to effect on-stream analysis, if 
computing capacity is available to analyze the record 
quickly. Such a device is the infrared spectrophotometer. 

A typical record from an infrared spectrophotometer 
represents the “frequency response” of a mixture of ma- 
terials; that is, the relative amount of light absorbed is 
plotted vs wavelength. One might guess from the com- 
plexity of detail of such a record that there is enough 
information to identify a large number of complicated 
components and, indeed, this is the case. However, the 
extraction of the desired information from such a record 
involves a considerable amount of computation. The 
Spectro-Sadic introduced about 5 years ago is an im- 
portant specialized data acquisition device. It converts 
an essential part of the information, represented by a 
(mass) spectrograph chart, into digital language that is 
acceptable to a computer. 

When spectrographic methods are used in stream 
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analysis, the identities of the different stream compo- 
nents are usually known, but their relative amounts are 
unknown. In this case, the computational job is merely 
that of multiplying a vector by a matrix. Twenty to 
thirty components in the vector is typical. Such a prob- 
lem, of course, can be handled and checked in a few sec- 
onds by most intermediate sized electronic digital com- 
puters. A matrix inversion is required for qualitative 
analysis, and this, of course, takes a little longer. How- 
ever, in order for the information to be most useful, the 
computer must be available whenever it is required. 
Machines with the “program interrupt” feature can 
efficiently handle such spectrographic analysis and also 
be available for other computations involved in the con- 
trol of a process or accounting work. 

Spectrographic instruments are a powerful analytic 
tool, but their output represents a cross-coupled and 
inter-twined presentation of data in what is basically 
an inconvenient form. The use of a computer to make 
sense out of such data and introduce it into the real- 
time process control loop has served as a prime example 
of the way digital computers can be used. On the other 
hand, simpler, more effective continuous stream analysis 
methods (for example, nondispersive infrared and ultra- 
violet analysis or gas chromatography) are now becom- 
ing available. These eliminate the need for computation 
to determine the quantitative chemical composition. 
Freed of the data reduction task involved in spectro- 
graphic analysis, digital computers may bring their con- 
siderable talents to bear on the important job of com- 
puting control actions based on stream analyzer data, 
combined with other pressure, temperature, and flow 
data. 

Present data logging speeds of 1 to 10 data samples 
per second are adequate for the alarm limit scanning 
and data logging functions for which the equipment has 
been designed. It seems likely that higher speed sam- 
pling will be required to allow the most to be made of 
sophisticated, computer-aided,  statistically-oriented, 
investigations aimed at extracting the control relation- 
ships from test data taken on existing pilot or produc- 
tion operating units. 

Real gains in plant operating economy are possible if 
production is freed from the constraint of steady-state 
operation. Variations in feed stream content, or a 
change from one product to another, will no doubt re- 
quire higher information processing rates. Technically 
feasible, but costly means are available for increasing 
sampling speeds by several orders of magnitude. 


SYSTEM RELIABILITY 


One modern refinery in Southern California is a high- 
production operation spread over many acres. As is 
usual, the process is regulated and the operating data 
are presented by means of instruments in a centralized 
control room. One of these instruments recently failed 
in a very subtle manner. This “upset” took 5 days to 


1958 


Gunning: Computers in 


straighten out, and cost the company on the order of 
3 million dollars in lost production. This illustrates not 
only the need for more sophisticated process control in 
one of our most advanced industries, but also underlines 
the need for absolute reliability, day after day, in every 
element in the control system. 

The mean-free-time between failures of the analog 
control equipment presently used is of the order of 6 
months. The mean-free-time of large present-day, high- 
speed, general-purpose, electronic, digital, computing 
machinery using vacuum tubes, is of the order of one 
day. Industrial management has been justifiably reti- 
cent about conversion to the more sophistocated, but 
less reliable hardware, and of assigning to it the in- 
creased system responsibilities required by more com- 
plete automatic control. 

Much is being done through the use of solid-state de- 
vices and other improved engineering techniques to ex- 
tend the reliability of basic digital computing equip- 
ment. Much more needs to be done before such equip- 
ment can be considered as a serious contender for 
real-time control. However, the inertial properties of 
most controlled systems open several avenues to an 
“effective” system reliability, wherein the failure of a 
part of the system may be permitted so long as it has 
negligible consequence to the process. 

Sophisticated error detection and correction tech- 
niques can be included since the basic information rates 
in the process being controlled are low when compared 
with the information processing rates available in elec- 
tronic computers. Marginal checking techniques have 
proved to be helpful in predicting system failures caused 
by components that drift or change gradually. In order 
that such techniques may be used, it is necessary either 
to have two computers or to be able to allow the system 
to “coast” over marginal checking periods of a single 
computer. “Breathing spells” in the schedule of a digital 
computer can be obtained by using it primarily to corre- 
late information and to produce instructions not for the 
process, but for analog servoregulators similar to those 
now used in the process industry. The process itself 
would continue to be regulated for minor variations and 
disturbances by the servoregulators, although between 
corrections by the digital system it would not necessar- 
ily be operating at optimum. During such intervals, the 
digital equipment could be marginally checked or minor 
repairs undertaken. 

Such an analog-digital computing and controlling 
team will afford new opportunities for reliable operation 
and error detection. The analog equipment can easily be 
arranged to reject set-point instructions from the digital 
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computer that are “unreasonable.” The digital equip- 
ment would, of course, be signalled concerning any such 
rejection. Furthermore, failing or drifting instruments or 
controllers can be detected by proper programming and 
mathematical model formulation. The short-time re- 
sponse of a thermocouple to temperature variations is 
perfectly adequate to assure the stability of tempera- 
ture control. The long-time drift of the voltage-to-tem- 
perature conversion relation of a thermocouple or the 
drift in the characteristics of the associated temperature 
regulator, will allow the controlled temperatures to shift 
slowly. Such a subtle shift will usually produce a change 
in product quality. Since stream-analysis information 
will be a part of the over-all control feedback loop, ad- 
justment in the regulator control point can compensate 
for such thermocouple or regulating device drifts. 

It appears that it is possible, therefore, to consider 
using a single digital computer in a control loop if 1) 
error detection methods, 2) error and failure prediction 
methods, and 3) modular design techniques are incor- 
porated. This last feature allows rapid maintenance to be 
performed by inexperienced personnel while the process 
coasts on the last set of instructions generated by the 
computer. It is conceivable that the computing system 
could itself schedule the optimum time for maintenance 
by predicting periods of smooth operation when no 
change in raw materials or other disturbances are called 
for by the over-all program. 


CONCLUSIONS 


The process industry is in the early stages of develop- 
ing and applying accurate mathematical models and 
computer techniques to aid in the design of new plants. 

Real-time control of plants using digital computer 
techniques is being approached from two directions. 
Digital computer mechanized operations research 
methods are scheduling the operations of large systems. 
Data acquisition equipment is introducing digital tech- 
niques at the detailed minute-by-minute operating 
level. The middle ground of combining primary instru- 
ment data with economic data in real-time control is 
still largely unexplored. The pieces of the large picture 
are being fitted together like a jig-saw puzzle. Many of 
the economically important pieces are already in place. 

A combined computer control system that 1) uses the 
best of today’s analog and digital building blocks, and 2) 
takes advantage of available redundancies, can achieve 
an operating reliability greater than that possible using 
either technique alone. Such a system will allow a prac- 
tical realization of complete, real-time computer control 
of an important part of our productive economy. 
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Aspects of Real-Time Simulation* 
WALTER F. BAUERT 


Summary—Present day digital computers are too slow for the 
comprehensive real-time simulation of complex electronic systems 
such as those involving guided missiles. Computers appearing 2-3 
years hence will be sufficiently fast for these applications. There is 
and will continue to be fertile application fields for analog-digital 
computer combinations. An existing large-scale (Univac Scientific 
1103A-Epsco converter—Electronic Associates) analog-digital system 
is described. System design and specifications are discussed which 
relate to real-time speeds, necessary accuracies, digital computer 
programming, and synchronizing the digital computer with the re- 
maining system components. Techniques for exploiting each com- 
puter’s advantages are mentioned. Applications in the guided missile 
field and in other fields are discussed. 


the computing groups at the Ramo-Wooldridge 

Corporation made the decision to establish a 
combination analog-digital computer system for the real- 
time simulation of ballistic missile systems. The plans 
and specifications subsequently made will come to frui- 
tion next month with the delivery of a high-speed, high- 
capacity analog-digital converter which will effect the 
tie-in of a large-scale digital computer and a large-scale 
analog computer. In the following I would like to de- 
scribe this system to you and comment on three addi- 
tional items related to this large-scale system: tech- 
niques for using such a large-scale analog-digital sys- 
tem; advances in digital computer systems which will 
make them competitive with analog computers and ana- 
log-digital combinations for simulation; and the future 
of analog-digital combinations. 

The comprehensive simulation of a ballistic missile 
system implies the precise handling of many low-fre- 
quency quantities, as well as the handling and integra- 
tion of many quantities varying at high frequency. Fur- 
thermore, certain aspects of the system are inherently 
of a sampled data nature and involve discrete data, 
while other parts of the system are inherently of a con- 
tinuous nature. Also, the requirement exists that all 
computations proceed at extremely rapid rates in view of 
the need to make tractable the simulation of a large 
number of “flights.” The requirement is, of course, a 
direct consequence of the fact that a number of meas- 
ured quantities are degraded by noise and a sampling 
must be made to determine the statistics of the random 
variables of the flight under various conditions. 

The Univac Scientific 1103A computer which has been 
performing much of the simulation is too slow to perform 
the entire set of computations of the extended simula- 
tion model at the desired rate, that is, a real-time rate. 
The large-scale Electronic Associates analog installa- 


i PPROXIMATELY two and one-half years ago 


* Presented at Symposium on Computers in Simulation, Data 
Reduction, and Control, New York, N. Y., March 21, 1957. 
+ Ramo-Wooldridge Corp., Los Angeles, Calif. 


tion, on the other hand, is too inaccurate to provide the 
precision necessary. Thus, despite the fact that these 
computers represent nearly the ultimate in commer- 
cially available computers in their respective classes, a 
large-scale analog-digital converter was ordered from 
the Epsco Company to provide a hybrid system with the 
desirable speed and precision. 

The characteristics of the Univac Scientific Model 
1103A computer and the 300-amplifier Electronic Asso- 
ciates computer installation are not emphasized here, 
for their characteristics can be found elsewhere. The 
converter equipment has 15 channels of analog-to- 
digital conversion and 15 channels of digital-to-analog 
conversion capacity. Ninety microseconds are required 
for conversion of all 15 channels from analog to digital 
quantities. The accuracy of the equipment will be 10 
bits on a dynamic range of 18 bits which implies a 0.1 
per cent or 1 millivolt accuracy, whichever is larger. 
Data to and from the channels can be presented or re- 
ceived at the 1103A computer under digital computer 
control with flexibility as to the channels chosen. Con- 
versions can be made upon command of the digital 
computer or upon a signal received from an external 
clock source. Specifications for the conversion equip- 
ment were written to allow the data transfer between 
the 1103A and the converter to be done quickly and 
with a minimum of digital computer programming 
inconvenience. The characteristics of the control and 
some of the techniques for programming, as well as a 
description of the system, have been recorded in a 
paper.’ 

One of the design features of the conversion equip- 
ment will exploit the “program-interrupt” feature of the 
1103A computer. This feature allows a signal from an 
external source to interrupt the computations. In usual 
operation, an external clock signals the conversion of all 
analog-digital channels and all digital-analog channels, 
and simultaneously interrupts the 1103A computer to 
“inform” it that new data is now available as a result of 
the conversions. With this technique the programmer, 
in laying out the computations, need pay attention only 
to the logic of the program and not to the timing of 
operations; in other words, the computations in the 
digital computer can proceed approximately asynchro- 
nous to the conversion frequency and the tedious count- 
ing and time estimation on the part of the digital com- 
puter programmer is unnecessary. 

As to the techniques of use of this system, certain 
loops of the simulation which do not have a profound 
effect on those variables, which need to be measured 


*W. F, Bauer and G, P. West, “A system for general purpose 
analog-digital computation,” J. Assoc. Computing Machina vol. 4, 
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precisely, are handled in the analog computer. The loops 
where precision is important are handled in the digital 
computer. In some cases the integration of certain vari- 
ables is actually accomplished by using both the digital 
and analog computers. Recently? some of these tech- 
niques have been presented. The technique of using 
the analog and digital computers together consists es- 
sentially in separating the signal into low- and high- 
frequency components and handling the low-frequency 
componentsdigitally, and the high-frequency components 
analog-wise. The ideas here are new and have neither 
been explored nor exploited to any significant extent, 
although the general technique of separating a signal 
into high- and low-frequency components is not new. 

A number of recent advances in digital computers en- 
hance the digital computer as an instrument for real- 
time simulation, control, and data reduction. One of 
these is the program-interrupt feature mentioned above 
which, to my knowledge, was first placed on the 1103 
computer installed at the N.A.C.A. installation in 
Cleveland described by Turner and Rawlings.’ This 
feature will undoubtedly be seen more and more in 
digital computers for all of these real-time purposes, for 
so frequently input to the computer occurs at an un- 
predictable time or in asynchronism to the computer 
operation. Another improvement beginning to be seen 
in digital computers is the time clock which can be in- 
terrogated under program control. This time clock can 
be a digital clock inside the digital computer, or an 
“analog” clock which must be sampled by the digital 
computer by means of an analog-digital converter. 
The first digital clock installed on a computer was that 
installed on the Whirlwind computer at the Massachu- 
setts Institute of Technology. A third major area of im- 
provement is the increased input-output speed allowed 
by the buffering systems, or the simultaneous compute- 
while-transfer operation which will appear first in a 
commercial model in the IBM 709 and Univac Scientific 
1105 computers next year. The increased high-speed 
storage capacity is likewise important, for in the digital 
computer storage can always be traded for time. 

The speed of the computer is perhaps the biggest 
item for real-time control. I recently heard of a large 
analog installation performing a real-time simulation of 
a complex system in approximately one minute per run. 
This simulation is probably the most complex ever 
attempted in this country. Because a digital check was 
necessary, the computations were completely pro- 
grammed for a drum computer. As I recall, it required 
20-odd hours to perform one simulation run. Consider- 
ing a factor of 50 as the difference in speed between this 
drum computer and a computer of the 1103A or 704 


2 J. L. Greenstein, “Application of AD-DA Verter System in Com- 
bined Analog-Digital Computer Operation,” presented at Pacific 
General Meeting AIEE, June, 1956. 

3. R. Turner and J. H. Rawlings, “Realization of randomly 
timed computer input and output by means of an interrupt feature,” 
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variety, these latter computers could perform the com- 
putations in about 25 minutes. The Univac-LARC com- 
puter, which will appear in late 1958, being about ten 
times as fast as the 1103A or 704, would perform the 
computations in about two and one-half minutes. The 
IBM STRETCH computer, which will probably appear 
in 1960, will be able to perform the computations for 
one simulation run in approximately one-half minute, 
that is, twice as fast as real time as with the analog com- 
puter and probably with considerably more accuracy. 

Although indicative of things to come, these facts do 
not necessarily carry with them the implicit forecast 
that the digital computer will make the analog computer 
completely obsolete early in the next decade. For one 
thing, the cost of these ultra high-speed digital com- 
puters will be prohibitive for many uses. 

It is likewise interesting to contemplate the future of 
analog-digital systems. Before the ultra high-speed 
digital computers appear in considerable numbers, it 
appears that the only possibility for the speed required 
for comprehensive real-time simulation of complex 
guided missile systems will come from these analog- 
digital combinations. It is entirely possible, also, that 
these systems will be very much competitive economi- 
cally with the large-scale digital systems for many, many 
years to come. The analog-digital system carries with it 
the promise of flexibility for, obviously, the analog and 
digital computers can each be used separately on a wide 
variety of problems. 

Many other possibilities of use appear. As examples, 
the digital computer can provide precision multiplica- 
tion, function generating, or precision time-delay func- 
tions on a time-shared basis while it is performing com- 
putations independent of the analog computer. The 
availability of the conversion equipment will, of course, 
allow for playing FM-FM telemeter tapes into the 
digital computer for computations required in data re- 
duction. 

I would like to pose this question which I hope will be 
provocative: is it possible that the $200,000 analog- 
digital installation consisting of $50,000 of analog equip- 
ment, $50,000 of conversion equipment, and $100,000 
digital computer will have more over-all computing 
potential than a $200,000 analog installation or a 
$200,000 digital installation? When one regards the 
digital computer as an equivalent number of precision 
multipliers, function generators, or precision time de- 
lays, the prospect may not seem too exciting; however, 
when one considers that the digital computer can be 
used for all of these things at different times as well as 
supply certain logical operations not possible on the 
analog computer, prospects seem considerably brighter. 

It is platitudinous to remark on the importance of 
simulation to the complex electronic system to a tech- 
nical audience. The expenditure of large sums of money 
for computing equipment to investigate guided missile 
characteristics is trifling compared with the money nec- 
essary to perform the equivalent number of firings of 


136 


the actual missiles. I believe it can be stated unequivo- 
cally that the modern guided missile would not be pos- 
sible without the modern simulation device. The actions 
of the complex system are only clearly and reliably un- 
derstood when the complete mathematical formulation 
of the entire missile system has been accomplished and 
when the resulting equations have been solved a large 
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number of times to determine the correct parameters 
and the correct detailed methods of missile guidance. 
It is tempting to say that the peaceful future of the 
free world depends on modern weapon technology and, 
therefore, depends directly on the electronic computer. 
As one with a heavy bias toward computers, I shall re- 
sist that temptation. 


Digital Information Processing for 
Machine-Tool Control’ 


ALFRED K. SUSSKINDt 


Summary—Filexibility, accuracy, elimination of manual skill, and 
greater work potential result from the addition of digital data process- 
ing to machine-tool control. Several special-purpose digital com- 
puters for machine-tool control have now been developed. These 
perform only simple computational tasks. The more sophisticated 
tasks require the addition of general-purpose large-scale computers. 
This combination, while only partially explored, appears very promis- 
ing. It should lead not only to improvements in translating a part 
design into the finished piece, but should ultimately assist in the 
design procedure. 


INTRODUCTION 


OR several decades now automatic machine tools 
is which the tool is controlled by analog methods 

have been successfully used. Here the instructions 
take the form of a master which is to be duplicated, or 
templates, cams, limit stops, etc. In the use of machines 
so controlled, considerable effort, both in labor and 
time, is incurred in tooling and setup. The first term 
refers to preparing the instructions (e.g., making the 
masters) and the second refers to inserting them into 
the machine (e.g., properly aligning several templates). 
Where production in substantial quantities is involved, 
the difficulties inherent in the analog control techniques 
have been of no serious consequence and for quantities 
such as are involved in automobile manufacture, spe- 
cial-purpose analog-controlled machines have amply 
demonstrated their high merit. But where production 
on a modest scale is involved, such as in the aircraft in- 
dustry where perhaps a few tens of units are produced 
before design changes occur, there has arisen a need for 
new control techniques which are more flexible and 
more adaptable to rapid change. At the same time, de- 
mands for accuracy have become increasingly stringent. 
These needs can now be met through the use of tech- 
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niques which control machine tools on the basis of digi- 
tal information. 

Adding digital information processing to machine- 
tool control results in several desirable features. First, 
the control instructions are in a form such as punched 
paper tape, cards, etc., all of which can be quickly 
changed when a different part is to be manufactured. 
Second, there is no limit to the accuracy with which the 
instructions can be specified or processed, and the limi- 
tations on the accuracy of the machined part are due 
entirely to the characteristics of the basic machine tool. 
Third, there is no manual skill involved in preparing the 
instructions such as there is in making masters, tem- 
plates, etc. Fourth, there is no longer any limit on the 
machine-tool functions (number and type of tool mo- 
tions, number of cutting tools, number and type of 
supervisory functions such as coolant ON-OFF, etc.) 
which can be brought under automatic control. Fifth, 
the machining can be performed at all times under de- 
sired optimum conditions (e.g., any prescribed feed rate 
within the capability of the machine tool may be ob- 
tained). Sixth, the potential exists for including in the 
information processing more than the faithful transla- 
tion of previously determined dimensions. While we do 
not know yet how to exploit this potential, it seems 
reasonable to expect that means might be found for 
letting the digital information-processing equipment 
contribute to certain parts of the design process. Thus, 
it would not appear out of the question to expect that, 
given the dimensions of mating surfaces and mounting 
requirements, the computer could design the part for 
minimum weight in light of prescribed stress require- 
ments, as well as determine appropriate machine-tool in- 
structions. 

While this last point may not deserve the prominence 
given to it here, it serves to emphasize that the tech- 
nology of digital machine-tool control is still far from 
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Fig. 2—Parts machined on experimental system. 


completed. This is certainly not surprising, for the in- 
clusion of digital information-processing equipment ina 
machine-tool system is a recent development. It started 
in 1949, when the Servomechanisms Laboratory of the 
Massachusetts Institute of Technology undertook the 
development of a numerically controlled milling ma- 
chine. 


Tue M.1.T. EXPERIMENTAL SYSTEM 


An experimental system was completed in 1952 and 
is shown in the photograph of Fig. 1. The machine tool 
is a standard model with three motions, x, y, and 2, 
which form a Cartesian coordinate system. The informa- 
tion-processing system contained in the vertical panels 
accepts at its input a series of instructions, each of 
which contains four pieces of data, punched on paper 
tape: the dimensional increments Ax;, Ay;, and Az; 
through which the center of the cutting tool is to move 
from where the previous instruction has brought it, and 
the time interval At; during which that motion is to be 
completed. The computing circuits then translate the 
instructions into a set of three pulse trains (one for each 
of the three motions), with as many pulses in each train 
as there are half-thousandths of an inch (the sensitivity 
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Fig. 3—Workpiece and corresponding tool path. 


arbitrarily chosen here) in the corresponding increment; 
and with all three pulse trains nearly uniformly spread 
in time over the prescribed time interval. Finally, each 
pulse train is converted into an analog signal (here a 
shaft rotation), suitable as input to the corresponding 
drive mechanism on the machine tool. Since each in- 
struction results in pulse trains and hence in shaft 
rotations and machine motions at a nearly constant 
rate, it follows that for each instruction the tool moves 
in a straight line. As a result, curves in this system are 
approximated by straight lines (e.g., circles by poly- 
gons). Each straight line may be arbitrarily small so 
that the approximation can be made as good as desired. 

A large variety of parts, most of them designed for 
aircraft, have been made on this experimental system. 
Some of these parts are shown in the photograph of 
Fig. 2. 

It should be emphasized that because the controlled 
quantities are the machine-tool drives, the proper mo- 
tions of the machine slides which will result in the de- 
sired cuts must be prescribed rather than the cuts them- 
selves. Because of the nonzero dimensions of a milling 
tool, the slide motions are not equal to the desired cuts. 
This is illustrated in Fig. 3 for a planar series of cuts, 
where KR is the radius of the cutting tool. To cut the path 
AB, the center of the tool must follow the path A’B’, 
aud the specification of the latter is prescribed on the 
punched paper tape in the form of its components. Two 
sets of specifications are required for this example. If it 
is assumed that the path A’B’ lies in the plane of the 
machine-tool table, the instructions consist of Ax,, Az, 
and Af, (the time during which the first motion is to be 
carried out), and then Ax, Aze, and Af. 

Since the machine-tool motions are the controlled 
quantities, the control system can assure only that the 
slide motions are carried out as specified. But this does 
not guarantee that the finished workpiece will have the 
desired dimensions; indeed, because of tool deflection, 
part warpage, etc., the workpiece will differ somewhat 
from that which would result from ideally rigid com- 
ponents. These deviations could be overcome only by 
monitoring the workpiece itself. However, conditions at 


Fig. 4—Director system. 


the cutting surface due to the presence of chips and 
coolant make feedback from the workpiece impossible. 
Not only does this restriction limit the accuracy with 
which parts can be made, but it also introduces an addi- 
tional step in the information processing, for it requires 
that the dimensions of the desired part be translated 
into the proper tool path. Tool-center correction, as this 
translation is called, is thus an inherent requirement in 
information processing for milling-machine control by 
digital techniques. 

In the experimental system, the information-process- 
ing equipment performs two functions: it carries out a 
linear interpolation between given points on the tool- 
center path in digital form, and it translates the result 
into an analog signal. There is no feedback between the 
machine tool and the director, as the information- 
processing equipment is called. It follows that separa- 
tion of the director from the machine tool does not 
affect operation. If the output of the director were re- 
corded on some medium such as magnetic tape and the 
magnetic tape were used to control the machine tool, 
the director would be free to prepare magnetic tapes for 
other parts or other machine tools, and it could operate 
in a convenient time scale with the magnetic tape play- 
back speed appropriately adjusted to correspond to the 
difference in time scale between director and machine 
tool. This is the approach taken in a second director 
developed by the M.I.T. Servomechanisms Laboratory 
for the Giddings and Lewis Machine Tool Company and 
completed in 1955. Fig. 4 shows a picture of the equip- 
ment. 


Tue GIDDINGS AND LEwis DIRECTOR 


While basically similar to the earlier experimental 
system, this second director contains several advances. 
A block diagram of this machine is shown in Fig. 5. The 
output of the director is a multichannel tape, with a 
group of six tracks containing the analog signals for 
controlling up to five machine-tool motions and the re- 
maining seven tracks containing supervisory signals 
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Fig. 5—Director system block diagram. 


such as coolant ON-OFF, etc. The motion information 
is in analog form, expressed as phase shift of a 200-cycle 
carrier. Parameters are chosen so that for each 0.1000 
inch of travel of a given machine-tool motion the cor- 
responding signal is shifted 360 degrees in phase with 
respect to a reference 200-cycle signal which occupies 
the sixth track of the first group. 

The reference and phase-shifted signals are generated 
electronically by a group of five decoders which receive 
at their inputs the corresponding pulse trains. In this 
machine too, each pulse train contains as many pulses 
as there are half-thousandths of an inch of travel in the 
desired machine motion and the pulse train is also 
nearly uniformly spread over the specified time interval. 

Next in the block diagram is the linear interpolator. 
Its function is to generate the five pulse trains in re- 
sponse to the specified Ax;, Ay;, Az;, Ag;, AO;, each spread 
over the same time interval At;. The interpolator here is 
similar in principle to that used in the experimental di- 
rector and is described by McDonough and Susskind [1]. 

The input to the director takes the form of punched- 
paper tape. Here the standard Flexowriter code, with 
a minor modification, is used, so that every digit of the 
instructions occupies one line of tape. Reading is ac- 
complished by means of a standard Flexowriter reader, 
which makes information available serially, line by line. 
But the interpolator requires the entire instruction in 
parallel and, in order to avoid a discontinuity in the in- 
formation processing, a new instruction must be im- 
mediately available in complete form when the previous 
set of pulse trains is completed. It is the function of the 
storage unit to retain an instruction until no longer 
needed and to assemble the next instruction by per- 
forming a serial-to-parallel conversion while the previ- 
ous instruction is being interpolated. In the Giddings 
and Lewis director the storage unit is made up of 
square-loop magnetic material cores, assembled as 
stepping registers for the serial-to-parallel conversion. 

The block diagram of Fig. 5 also shows a box labelled 
Translator. This unit carries out a translation between 
Flexowriter code and the internal language of the di- 
rector. While both are in binary coded decimal form, 
the use of a different internal code (5-2-1-1) leads to 
greater computational ease inside the director. The 
unit marked Switch separates instructions designating 
the motions (Ax;, Ay;, Az;, Ag;, Ay;, At;) from supervisory 
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instructions. The latter are directly translated by the 
Control-Signal Generator into tape-recording excitation. 

Finally, there is a unit labelled Visual Display, which 
monitors the operation of the decoders. It puts on dis- 
play, in decimal form, the distance in inches through 
which each of the five machine-tool motions has been 
commanded to move from the start of the program. This 
feature allows human observation of the information 
processing and is an aid in program debugging and di- 
rector trouble shooting. Decimal representation was 
chosen here and in the input (paper tape) for conven- 
ient man-machine communication. 

The director also contains a group of control circuits 
which perform such functions as power-supply control, 
initial clearing of the system, time-scale settings, and 
automatic slow-down of machine-tool commands when 
successive instructions would otherwise result in ex- 
cessive accelerations. Means for manually modifying all 
programmed feed rates and inverting sign conventions 
are also provided. The latter permits processing of right- 
and left-handed parts from the same punched-paper 
tape. 


OTHER DIRECTORS 


Both of the directors described above perform but 
one major computational task: linear interpolation be- 
tween points in space through which the center of the 
cutting tool moves. The question may well arise 
whether or not so limited an ability is adequate for 
milling surfaces. The answer is affirmative if the execu- 
tion of useful work at reasonable cost is accepted as the 
criterion. But this is not to say that directors with higher 
computational ability might not be more successful. 

For example, it is not difficult to imagine a director 
system which accepts at its input a series of points on 
the workpiece itself rather than on the required tool path, 
and which computes the required tool path that will 
result in connecting these points by a second degree 
curve such as a parabola or circle. This approach is de- 
scribed in the literature [2]. One could extend the com- 
putational ability of the director even further, so that 
the required input instructions are reduced to the 
point where either the workpiece is described only in 
terms of the equations of the surfaces which define it or, 
in the absence of such equations, by a series of key 
points such as are now given at successive stations of an 
airplane wing. In the latter case, the director would 
derive appropriate equations of any degree connecting 
them. This director should also be capable of choosing 
proper feed rates in view of machining requirements 
such as available spindle-drive horsepower and _ per- 
missible cutting tool loads, and in view of required work 
tolerances and dynamic characteristics of the machine- 
tool drive servomechanisms. Clearly, an information- 
processing problem of considerable magnitude results. 
It appears that its solution requires the use of machines 
comparable to the large-scale computers (1103, 704, 
etc.) and that the flexibility of these machines is 
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needed. One might then conclude that because these 
computers are already in existence and becoming more 
easily available there is no need for developing sophis- 
ticated machines aimed specifically toward the control 
of machine tools. It may well be that the standard com- 
puters, coupled with simple directors which perform 
the final data expansion (linear interpolation in the ex- 
amples described above), form the optimum informa- 
tion-processing system for machine-tool control. The 
director is then akin to a computer terminal device 
which takes over the time-consuming, but fairly simple, 
last steps of the processing. 


DIRECTORS AS COMPUTER TERMINAL DEVICES 


Studies have started at the M.I.T. Servomechanisms 
Laboratory which are directed toward making use of 
large-scale computers in the information processing for 
machine-tool control. The first phase of this work, car- 
ried out by Siegel [3], has considered only the case of 
planar configurations consisting of circles and straight 
lines. It has resulted in the preparation of programs for 
the Whirlwind computer which enable persons with no 
previous knowledge of digital computers to use Whirl- 
wind to prepare complete instructions for the experi- 
mental director. Given the specifications of the circles 
and straight lines which make up the desired part, the 
size of the cutting tool, the tolerance required, and the 
desired machine-tool feed rates, all stated in language 
closely approximating ordinary English, the computer 
has been programmed to prepare the control tape for 
the director. In its operation, the computer first trans- 
lates the input language, chosen from the point of view 
of human convenience, into its own language and de- 
duces the nature of the computations needed. Next, the 
computer determines the complete contour of the part 
to be made. In the third step, the computer finds the 
straight lines which approximate the contour within the 
accuracy specified. Then Whirlwind computes the cut- 
ter-center path corresponding to the straight lines. 
Next, the tool-center path, together with feed-rate 
specifications, is coded into the form required by the 
director. Finally, the complete commands are printed 
out. on punched paper tape, ready for insertion into the 
director. 

At the present time, the Laboratory is extending this 
work to the general three-dimensional case, without re- 
striction on the types of functions which specify the 
desired surfaces. Computer programs will ultimately be 
written for generally available machines, e.g., 704, 
rather than the Whirlwind computer. Another extension 
will consist of including the capacity to control machine 
tools with more than three motions. Again, the aim is 
to construct programs in such a way that they can be 
used by personnel with no previous computer experi- 
ence. It is expected that these programs will be of maxi- 
mum benefit if used by design or production personnel 
without familiarity with computer technology. 

Studies in the area of utilization of general purpose 
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digital computers for determining director instructions 
are now being carried on at several places, all primarily 
concerned with aircraft production problems. This work 
is in the early stages, and it will be some time before 
significant results can be expected. 


OTHER Tyres OF MACHINE TOOLS 


So far, the discussion has centered on milling ma- 
chines. For these tools, information processing require- 
ments are greatest. The path of the tool must be con- 
trolled at nearly all times since the tool is in contact 
with the work and is removing metal almost continu- 
ously. There is a large group of other machine tools in 
which the tool is in contact with the work only intermit- 
tently. 

Consider, for example, a drill press making a series of 
holes in a plate. The tool is first placed over the desired 
location of the hole. Next, the tool travels into the work 
to the desired depth and is then withdrawn. In the next 
step, the tool is put into position for a second drilling 
operation, and the process continues. Note that there 
is no longer a requirement for path control, for after 
the tool withdraws from one hole and until it is ready 
to plunge in for the next hole, no metal is being cut. 
Thus the control system need only assure that before 
plunging starts, the tool is at the desired coordinate 
and the path taken after the previous withdrawal is of 
no consequence except that it affects the time spent 
between successive drilling operations. 

It follows that for this drill press one need only specify 
the coordinates of successive holes, the depth of each 
hole, and the rate at which drilling is to occur. It is quite 
simple to design a control system which accepts these 
specifications in digital form, and even though the data 
processing required here is negligible, a digital control 
system may be advantageous because of the ease with 
which the input specification can be prepared, edited, 
and automatically and accurately read. Examples of 
such systems are given by other authors [4]-[6]. In 
these systems, the position instructions are read and 
then compared, either in digital form or in analog form 
after they have been decoded, with the feedback signal 
indicating the position of the controlled motion. The 
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result of the comparison is used to bring the controlled 
motion into alignment with the specified position. 


FUTURE DEVELOPMENTS 


So far, all applications of digital information process- 
ing to machine-tool control have been concerned with 
the production of parts which were designed for conven- 
tional production techniques. As the new technology 
becomes more fully exploited, it will lead to modifica- 
tions in established design procedures. For example, 
much effort is now being expended on producing draw- 
ings which are graphical representations of some of the 
part specifications. However, for digital information 
processing, a symbolic representation is needed, such as 
equations or tables of points defining the required part. 

It seems reasonable to expect that the benefits of 
digital information processing can be extended beyond 
the area now explored which covers only a portion of 
the over-all process of translating the need for a part 
into the finished product. A more comprehensive ap- 
proach should become the subject of future research. 
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Realization of Randomly Timed Computer Input 
and Output by Means of an Interrupt Feature’ 


L. R. TURNER} anv J. H. RAWLINGS} 


Summary—The Lewis Flight Propulsion Laboratory of the Na- 
tional Advisory Committee for Aeronautics is assigned the primary 
task of research into the basic problems of the design and operation 
of aircraft power plants. Since World War II this has primarily meant 
research in the field of jet and turbine engines. 

Although these engine types are simple in concept, it was quickly 
realized that they could not be developed adequately by research on 
small-scale models. It has been necessary to build large facilities in 
which research could be conducted under conditions approximating 
those expected of actual flight. 

With new facilities under construction or being planned, it be- 
came clear in the late 1940’s that only some form of mechanization 
could cope with the growing problem of data-processing. The result- 
ing studies led to a plan in early 1950 for the construction of a data- 
recording system aimed toward digital data-processing techniques. 

In January, 1952, the first model of an automatic digital pressure 
recorder with mechanized, manually-controlled computing was put 
into operation. 

From this point the system gradually evolved in concept and in 
working equipment to its present form in which a central editing and 
computing station serves a number of remote research facilities. The 
operation of data collecting and input editing is described in con- 
siderable detail by Ryskamp.! It is sufficient at this time to give a 
brief description of the system. 


SYSTEM DESCRIPTION 


CENTRAL Automatic Digital Data Encoder 

(CADDE) is connected by telephone cables to 

nine remote experimental stations. The CADDE 
system is designed to record up to 500 steady-state 
quantities in any record and, when it is available, auto- 
matically connects itself to the requesting facility and 
controls the sequence of operations required to record 
the desired voltage, pressure, and speed measurements. 
These are recorded together with identification and mis- 
cellaneous control information on magnetic tape, which 
then serves as a legal record of the experiment and also 
as an input store for data processing. The magnetic tape 
may also be used to produce high-speed plots of the un- 
processed data or to type selected information on an 
electric typewriter. 

The computing unit of the Central Data Processing 
System isa Remington Rand Univac Scientific Com- 
puter, Model 1103, which has been modified to adapt 
to the task of communicating with the data recorder 
and with remote devices for reporting of results, in order 
to achieve on-line data processing. 
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Mechanically, the most important problems to achieve 
on-line data processing arise from the wide difference 
between relatively high speeds required for editing and 
computing, and slow speeds of inexpensive devices for 
data transmission. Another factor is the awkward varia- 
tion in the rate at which results may be produced. This 
rate may vary from perhaps one hundred per second to 
as low as one every three or four seconds. 

Fortunately, the quantity of information which is to 
be displayed remotely for immediate inspection (per- 
formance coefficients or plots) is sufficiently small for 
adequate speeds to be obtained by multiplexing a few 
efficiently used slow-speed devices, such as the punches 
of Flexowriters, while the bulk output is directed to- 
ward higher speed devices for deferred printing. 

This degree of multiplexing is, however, not the com- 
plete answer, because the sum of computing time and in- 
put-output time is still greater than desired. It is neces- 
sary to interlace computing and information transfers 
to achieve required speeds. It is well known that almost 
all of the cycle time of input-output devices is available 
for computing if coding which is clever enough can be 
prepared. In practice, the difficulty of such coding is so 
great that little use is usually made of this theoretically 
available time. 

A satisfactory solution to these problems which re- 
quires no external buffers and generates no important 
coding difficulties was achieved by the development of 
the Lewis Laboratory’s version of Interrupt, which was 
first put into operation in February, 1956. 

In this version of Interrupt, programs in the high- 
speed store of the computer drive the input-output de- 
vices and provide for the necessary information trans- 
fers. These programs are in the form of subroutines 
which are referenced automatically at the time required 
by the device itself by means of the interrupt feature. 

Interrupt is permitted only when the interrupting de- 
vice is ready for an immediate information transfer. 
As each device for which interrupt has been enabled 
requires service, a signal to the computer causes the 
normal sequence of command procurement from the 
store to be interrupted and, at the same time, an arbi- 
trary command is planted in the input register for imme- 
diate execution. The command counter is not advanced. 
In nearly all cases the injected command will be the 
Return Jump command of the 1103 with the addresses 
of the exit and entry lines of the appropriate service 
subroutine. A properly coded subroutine will provide 
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for the next information transfer and will return com- 
puting to the point of interruption. 

While experience has shown that the interrupt mode 
achieves high efficiency of all components without add- 
ing significantly to the difficulties of coding, there are, 
however, a few important ground rules which must be 
strictly followed, both in hardware design and in coding, 
if an interrupt system is to work at all. The importance 
of adherence to the standards was forcefully given 
point during the early operation in the interrupt mode. 
It was learned that a problem, or nearly all of a group of 
similar problems, could be solved correctly with incor- 
rect service subroutines and with the actual omission of 
logically required hardware. 

As an example of a coding error which might go un- 
detected, it is clear that the failure to set an interlock 
against overwriting an output store before it has been 
translated and used would be unimportant if the output 
were always complete before the interlock was tested. 
Similarly, the complete failure of a hardware interlock 
intended to prevent a difficulty of rare occurrence might 
go undetected for a very long time. 

The heart of these difficulties is a subtle difference 
between interrupt and noninterrupt modes of operation. 
In the conventional mode, not only does the result of 
calculation depend only on the initially stored code and 
on the data, but so does every intermediate state of the 
computer. With interrupt, because of the random tim- 
ing of interruption, this is no longer true; in fact, even 
with a correct program, only the answers are expected 
to be predictable. One of the startling effects of this 
variability is that a programming or coding error, by 
producing a variable error or sometimes none at all, 
may be indistinguishable from a random computer mal- 
function. Many hours were spent in discovering this 
peculiarity. 

The elimination of these hazards requires well-coor- 
dinated care both in the hardware design and in coding. 
The problem can be most easily stated from the coding 
standpoint by considering the successive tasks of an 
interrupt invoked subroutine. It is assumed that at the 
time of interrupt the computer is doing useful work. 

The subroutine must first form a record of the con- 
tent of the arithmetic registers (A and Q in the 1103) if 
they are to be used for housekeeping calculations. It 
must then prepare for the next step of information 
transfer and correctly modify any interlocks required to 
keep the input-output and the computational routines 
in step. Finally, it is necessary to restore the status of 
the arithmetic registers and the program counter. A 
further interrupt by the same device may not be per- 
mitted until these tasks are completed. If an earlier 
interrupt were permitted, the program reference stored 
by the return jump would be to some point in the serv- 
ice routine and an unintended tight loop would result. 

From the hardware standpoint it is necessary to 
maintain an electronic connection between the input or 
output device and the computer until the transfer is 
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complete. It is also necessary to prevent an interrupt by 
any device that would refer to the subroutine in current 
use. In the 1103 these conditions were easily met. The 
information transfer command is used to enable a fur- 
ther interrupt and, at the completion of its execution, 
cause the current input-output device to be discon- 
nected. The information transfer command is the next- 
to-last command in the service routine. The last com- 
mand is a jump which restores the program counter to 
its content at the time of interrupt. 

The hardware contains a delay in enabling the next 
interrupt for one command after the information trans- 
fer command, thereby assuring that the problem status 
can be completely restored. The command structure of 
the 1103 is compatible with these simple controls be- 
cause the arithmetic registers need not be used for input 
and output. 

The electronic hardware executes several other im- 
portant functions. With several interrupting devices, it 
is necessary to control through coding the ability of each 
device to initiate interrupt. This is done in the 1103 by 
means of the External Function (EF) command. One 
EF command is executed to set an enabling flip-flop 
when a task is generated. A different EF command lo- 
cated in the service routine resets this flip-flop when the 
task is finished. 

Other functions of the hardware include the selection 
of a single device for interrupt when several may be 
competing for this service, and the usual safeguards 
against premature or improperly repeated information 
transfers. 

At the present time the interrupt feature is used by 
the control and reading mechanism of the tape handlers 
on the CADDE recording system and for two Western 
Electric high-speed tape punches. Equipment is near- 
ing completion to add four slow-speed punches, which 
will be used to transfer information through paper tape 
buffering to remote typewriters and digitally operated 
point plotters. Further details are given in the following 
sections. 


INTERRUPT 


When computation demands input or output, devices 
external to the computer are activated for interrupt. 
Any device, selected under program control, may inter- 
rupt the computer if it is ready for an immediate trans- 
fer of information. By interrupting the computer the 
normal sequence of command procurement is suspended 
and, instead, the interrupting device causes a command 
to be injected for immediate execution. Injected com- 
mands are emitted from a pluggable COMMAND 
EMITTER having a position available for each device. 

The command selected for injection is a return jump 
command because its execution results in recording the 
contents of the Program Address Counter on the oc- 
currence of interrupt. The injected commands provide 
entries to the service subroutines for the interrupting 
devices. 
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Devices which are ready for immediate information 
transfer compete for computer time. Conflicts arising 
from the fact that several devices may be ready at the 
same time are automatically resolved by an ANTI- 
COINCIDENCE NET. 


OPERATION 


The system is illustrated in a block diagram [Fig. 1(a) 
and 1(b) |, on the next page. Only two devices are shown 
because of space limitations. 

Consider now that the computer is in the normal 
mode of operation, executing commands from consecu- 
tive address locations in storage. An EXTERNAL 
FUNCTION COMMAND is executed, turning on the 
ability of the OUTPUT DEVICE to interrupt. If there 
has been a previous use of this device, a GATED 
TRANSFER COMPLETE will have cancelled the 
READY SIGNAL generated in the previous cycle. 

The OUTPUT DEVICE is not yet READY. When 
the OUTPUT DEVICE can accept new information 
immediately (its INPUT THYRATRON REGISTER 
cleared of previous information, and plate voltage re- 
applied), it sends a DEVICE RESUME signal to the 
READY CONTROLS. Having now received a signal 
turning INTERRUPT on, a GATED TRANSFER 
COMPLETE, and a DEVICE RESUME, the ready 
controls can generate a READY signal. If the ANTI- 
COINCIDENCE NET is empty (no established pri- 
ority), the READY signal is generated on the next 
NET ENTRY PROBE. If a second device has mean- 
while interrupted the computer, the ANTICOINCI- 
DENCE NET may not be empty. In this case the 
READY signal will not be generated until a GATED 
TRANSFER COMPLETE to the second device’s 
READY CONTROL cancels the second device’s 
READY signal and releases the ANTICOINCIDENCE 
NET. 

NET ENTRY PROBES are not generated while 
there is an established priority in the ANTICOINCI- 
DENCE NET. This insures that once a priority is es- 
tablished it is not disturbed until the device which has 
the priority is served with its share of computer time. 
After the transfer is complete, all devices which are 
READY compete for computer time. 

The READY signal from the OUTPUT DEVICE 
READY CONTROLS is sent to the ANTICOINCI- 
DENCE NET. If there are no other READY signals 
in the net which are higher than the READY signal 
from the OUTPUT DEVICE, then the OUTPUT DE- 
VICE will receive PRIORITY. 

The PRIORITY signal is sent to SEQUENCE CON- 
TROL. An enable is generated in SEQUENCE CON- 
TROL associated with the OUTPUT DEVICE and 
sent to the COMMAND EMITTER, which in turn 
enables the INPUT MIXER in preparation for the in- 
jection of a command into the INPUT-OUTPUT 
REGISTER. The COMMAND EMITTER is plug- 
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gable, and is plugged in the position reserved for the 
OUTPUT DEVICE with a return jump command. 
The v-address portion of the plugged command is the 
location in storage of the subroutine which will serve 
the OUTPUT DEVICE. The u-address portion is 
plugged with the address of the last command in the 
service subroutine. 

An INTERRUPT request is generated and sent to 
MAIN CONTROL of the computer where it will re- 
quest an alteration of the command procurement por- 
tion of the computer cycle. The request may be denied 
temporarily if the INPUT-OUTPUT REGISTER is in 
use or if the computer is in a repeat sequence. When the 
request is granted, the INTERRUPT signal transfers 
the contents of the INPUT MIXER (the injected return 
jump command) to the INPUT-OUTPUT REGISTER. 
The injected command is then sent from INPUT-OUT- 
PUT REGISTER to ARITHMETIC EXCHANGE 
REGISTER and then into the main control where it is 
executed. 

The INTERRUPT signal is also sent into SE- 
QUENCE CONTROL where it initiates an INTER- 
RUPT lockout. The lockout is not released until SE- 
QUENCE CONTROL receives a TRANSFER COM- 
PLETE (WRITE) from main control of the computer. 
This lockout prevents the OUTPUT DEVICE from 
interrupting its own service subroutine. 

When the injected return jump command is executed, 
the computer records the address of the command which 
would normally have been the next to be fetched for 
execution in the interrupted problem. It then jumps to 
the address given by the v-address portion of the in- 
jected return jump command and finds the subroutine 
which services the OUTPUT DEVICE. 

The next-to-last command in the subroutine is an in- 
formation transfer instruction. A WRITE ENABLE 
is generated as a result of the execution of this transfer 
instruction and sent to SEQUENCE CONTROL. SE- 
QUENCE CONTROL directs the WRITE ENA- 
BLE to enable the OUTPUT DEVICE THYRA- 
TRON REGISTER. The TRANSFER COMPLETE 
(WRITE) is also a by-product of the execution of the 
transfer instruction and is sent to SEQUENCE CON- 
TROL. SEQUENCE CONTROL directs this signal to 
the READY CONTROLS of the OUTPUT DEVICE 
as a GATED TRANSFER COMPLETE, cancelling 
the READY signal. The timing of INTERRUPT and 
TRANSFER COMPLETE signals prevents a further 
interrupt before the execution of the last command in 
the service subroutine. This guarantees that the inter- 
rupted problem is always resumed if the last command 
of the sburoutine is a jump to the re-entry address of 
the interrupted problem. 

Any device which is READY and HIGH in the 
ANTICOINCIDENCE NET may now interrupt the 
computer. Interrupt for the OUTPUT DEVICE is still 
turned on and it will interrupt again when it is READY. 
To turn off its ability to interrupt would require execu- 
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Fig, 1—Interrupt communication of INPUT and OUTPUT DEVICE with Remington Rand 1103 computer block diagram. 
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Fig. 2—READY CONTROL. 


tion of another EXTERNAL FUNCTION command. 
This command is normally located in the subroutine 
servicing the OUTPUT DEVICE. 


READY CONTROLS 


READY CONTROLS are necessary to insure that 
information is transferred only when immediate ac- 
ceptance is possible and that only one transfer is made 
per device cycle. 

READY, for a particular device, means that it can 
accept or supply new information immediately. In 
operation all devices in the system are tested periodi- 
cally to determine if they are ready. The testing is done 
by the NET ENTRY PROBE. The circuitry for this 
testing procedure is illustrated in Fig. 2. 

When a device interrupts the computer and is serv- 
iced by the execution of a subroutine, the execution of 
the information transfer command within the subrou- 
tine results in a GATED TRANSFER COMPLETE 
to the READY CONTROLS of the interrupting device. 

The GATED TRANSFER COMPLETE turns the 
READY signal off after its use so that the device can- 
not interrupt more than once per device cycle. When a 
device has completed its mechanical work, or when its 
buffer register is clear and ready for more information, a 
DEVICE RESUME is generated and sent to the 
READY CONTROLS of the device. Requiring a DE- 
VICE RESUME insures that the READY signal can- 
not be turned on again until the device has accom- 
plished its assigned task. 
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Fig. 3—ANTICOINCIDENCE NET. 


READY 
OUTPUT DEVICE 


Each device in the system must have itsown READY 
CONTROLS 


ANTICOINCIDENCE NET 


The ANTICOINCIDENCE NET is necessary be- 
cause more than one device can be READY at the same 
time. 

Each device in the system is assigned a position in 
the ANTICOINCIDENCE NET, which has priority 
assigning properties. 

READY signals from device READY CONTROLS 
are sent to the ANTICOINCIDENCE NET illustrated 
in Fig. 3. An output from any position in the net can- 
cels any other outputs from the net lower than itself. 
With multiple inputs, only the output which is highest 
in the net survives. 
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SEQUENCE CONTROL 


PRIORITY signals are sent to SEQUENCE CON- 
TROL, which is essentially an information traffic di- 
rector. It directs INFORMATION TRANSFER EN- 
ABLES and TRANSFER COMPLETE signals to the 
proper place and at the proper times. 

For command injection, it directs a COMMAND 
EMITTER ENABLE to the pluggable position of the 
COMMAND EMITTER reserved for the interrupting 
device at the time appropriate for command injection. 
At the appropriate time for information transfer, SE- 
QUENCE CONTROL directs a TRANSFER EN- 
ABLE to the input or output register of the device re- 
ceiving or delivering information. When a transfer has 
been completed, SEQUENCE CONTROL directs a 
complete signal to the READY CONTROL of the de- 
vice which has been served with computer time. 

The circuitry is illustrated in Fig. 4. The flip-flop 
shown is initiated by the occurrence of INTERRUPT 
and reset by any TRANSFER COMPLETE. Its func- 
tion is to sequence information transfer enables in the 
order in which they are to be used (for multiple trans- 
fers, a multistage flip-flop distributor would be required) 

INTERRUPT REQUEST is generated in SE- 
QUENCE CONTROL. It is locked out anytime FF I 
is set to “1,” preventing any device from interrupting 
its own service subroutine. 

NET ENTRY PROBES are generated in SE- 


QUENCE CONTROL and sent to all READY CON- 
TROLS. NET ENTRY PROBES are cut off by the 
presence of a priority signal. 

Expansion of SEQUENCE CONTROL to include 
more devices requires additional AND circuits to sup- 
ply one COMMAND EMITTER ENABLE for each 
device and one INFORMATION TRANSFER EN- 
ABLE for each device. One TRANSFER COMPLETE 
GATE is required for each device, and the OR circuit 
receiving device priorities must have an entry for each 
device. 


ALTERATION OF COMPUTER CYCLE 


Alteration of part 1 of the computer cycle results in 
fetching the next command to be executed from the 
INPUT-OUTPUT REGISTER rather than from the 
next consecutive address location in storage. 

The normal functions of command procurement occur 
on MAIN PULSE 6 and MAIN PULSE 7. When IN- 
TERRUPT occurs, the functions of MAIN PULSE 6 
are altered, but MAIN PULSE 7 retains its normal 
function. 

Alteration of the computer cycle is illustrated in Fig. 
S(a) and 5(b). If the INPUT-OUTPUT REGISTER 
contains information from a previous use or if a com- 
mand is being executed which will make immediate use 
of the INPUT-OUTPUT REGISTER, then a request 
for INTERRUPT is denied. 


—— 


—— 


ee 


| 1958 


Turner and Rawlings: Randomly Timed Input and Output 


147 


INTERRUPT 
TO SEQUENCE CONTROL 


AND 


VO 
ADVANCE MAIN PULSE . G 
DISTRIBUTER , 


INTERRUPT DEMAND 


TIMED INTERRUPT REQUEST 


FF I 


ie) | 


VO | 


INTERRUPT PROBE 


vo 4 
a INHIBIT 
O | 


INPUT MIXER 


STORAGE CLASS 
CONTROL RESUME 


TRANSMIT INPUT-OUTPUT 
REGISTER TO EXCHANGE 
REGISTER 


TRANSMIT 

PROGRAM ADDRESS CTR 
TO STORAGE ADDRESS 
REGtSTER 

ADVANCE 

PROGRAM ADDRESS CTR 


STORAGE CLASS CONTROL 
STORAGE REFERENCE 


STORAGE CLASS CONTROL 
FAULT DETECTOR 


CLEAR PROGRAM 
AND REPEAT SEQUENCE 


For | CONTROL 
poo} E= 3 e) O O O e) 
MPO iNTERRUPT INPUT-OUTPUT CP FOLLOWING MAIN PULSE MAIN PULSE INPUT -OUTPUT REPEAT 
REQUEST REGISTER NOT MP6 DIST. FF I-"O" DIST. FF O-"0" REGISTER IN INSTRUCTION 
IN USE (+) USE (+) (+) 
(WAIT WRITE) (MCT) 
(a) 
INTERRUPT NORMAL 
MP Housekeeping MP Housekeeping 
6 Transmit INPUT MIXER to INPUT-OUTPUT REGIS- 6 Transmit PROGRAM ADDRESS COUNTER to STOR- 
TER AGE ADDRESS REGISTER 
Set FF I, Fig. 4=1 Advance PROGRAM ADDRESS COUNTER 
Set FF II, Fig. 5(a) =1 Initiate READ 
Initiate READ Clear X 
Set WAIT INTERNAL WAIT INTERNAL 
Clear X 
CP Transmit INPUT-OUTPUT REGISTER to EXCHANGE | CP Clear PROGRAM CONTROL REGISTER 
REGISTER Clear REPEAT SEQUENCE CONTROL 
Clear FF I, Fig. 5(a) 
Clear FF II, Fig. 5(a) 7 Clear STORAGE ADDRESS REGISTER 
Clear PROGRAM CONTROL REGISTER Transmit EXCHANGE REGISTER TO PROGRAM 
Clear REPEAT SEQUENCE CONTROL CONTROL REGISTER 
STORAGE CLASS CONTROL resume Initiate 2-ysec delay 
if Same as NORMAL 


(b) 


Fig. 5—Alter part 1 of computer cycle. 


A request for INTERRUPT is denied when the com- 
puter is in a repeat sequence because the PROGRAM 
ADDRESS COUNTER does not contain the address 
of the next instruction to be fetched. 

The CLEAR for the PROGRAM CONTROL 
REGISTER and REPEAT SEQUENCE CONTROL 


has been shifted unconditionally from MAIN PULSE 6 
to the following clock pulse so that denials to INTER- 
RUPT REQUEST, based on the contents of the PRO- 
GRAM CONTROL REGISTER or REPEAT SE- 
QUENCE CONTROL, are effective throughout the 
time when INTERRUPT can occur. The shift in timing 
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does not affect operation of the computer in the NON- 
INTERRUPT mode. The purpose of FF I is to time 
the INTERRUPT REQUEST. In case the request is 
denied, it must be removed so that the normal sequence 
is not upset. FF I is reset unconditionally on the clock 
pulse following MAIN PULSE 6. 

INTERRUPT REQUEST and INPUT-OUTPUT 
REGISTER IN USE (+) are ANDED to enable 
GATE VO 1, so that the timed request does not occur 
if there is a WAIT WRITE. This is necessary because 
of the random timing of devices that clear WAIT 
WRITE and are not included in the INTERRUPT 
system. The entries to the OR circuit from the MAIN 
PULSE DISTRIBUTOR FLIP-FLOPS are permissive 
and specify the time at which INTERRUPT can occur 
(MAIN PULSE 6 and the following clock pulse). 

Generation of INTERRUPT DEMAND results in 
an INTERRUPT. 

INTERRUPT DEMAND inhibits the normal house- 
keeping function of sending the contents of the PRO- 
GRAM ADDRESS COUNTER to the STORAGE AD- 
DRESS REGISTER and advancing the COUNTER. 
The COUNTER is left intact on MAIN PULSE 6 when 
INTERRUP® “oceurs. On MAIN’ PULSE “6 (AD- 
VANCE MAIN PULSE DISTRIBUTOR) the con- 
tents of the INPUT MIXER are sent to the INPUT- 
OUTPUT REGISTER, READ is initiated (but not 
allowed to go to completion), and FF II is set to 1. The 
purpose of FF II is to inhibit the reference to STOR- 
AGE initiated by the READ and to inhibit the detec- 
tion of a storage class control fault. The fault would 
occur as a result of inhibiting the reference to storage, 
since the inhibition results in no storage class selection. 

On the clock pulse following MAIN PULSE 6, the 
contents of the INPUT-OUTPUT REGISTER are 
sent to the ARITHMETIC EXCHANGE REGISTER 
(X), and a storage class control resume is generated 
(WAIT INTERNAL is set on INITIATE READ). 
Flip-flop II is reset. 

On MAIN PULSE 7, the injected command is sent 
from the ARITHMETIC EXCHANGE REGISTER 
(X) to the PROGRAM CONTROL REGISTER for 
execution. This is a normal function for the 1103. 


GROUND RULES 


The following rules characterize the system and also 
indicate some of the design considerations. 

1) A device may not interrupt its own service sub- 
routine. INTERRUPT REQUEST (SEQUENCE CON- 
TROL) is locked out when INTERRUPT occurs and 
this lockout is not released until a TRANSFER COM- 
PLETE is given. The TRANSFER COMPLETE is a 
by-product of executing a transfer command. There- 
fore, the location of the transfer command, within the 
service subroutine, is important. It must immediately 
precede the exit from the subroutine, unless it can be 
shown that the device in question cannot become 
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READY again before the exit is made. 

2) A device may not interrupt more than once per 
device cycle. 

3) A device may not use time-shared equipment for 
longer than is necessary for acceptance of the interrupt 
request and execution of the service subroutine. 

Cancellation of the READY signal by GATED 
TRANSFER COMPLETE and release of the INTER- 
RUPT lockout (SEQUENCE CONTROL) by TRANS- 
FER COMPLETE insure rules 2 and 3 above. 

4) A device may not interrupt if there is an incom- 
plete previous assignment. This means completion of 
the actual mechanical assignment. DEVICE RESUME 
is the signal indicating completion and is required to 
produce a READY signal. 

5) An established priority may not be disturbed until 
it is cancelled after use. NET ENTRY PROBES are 
not generated while there is an established priority. 

6) If INTERRUPT is to be turned off for a par- 
ticular device the EXTERNAL FUNCTION command 
for this purpose must be located within the service sub- 
routine, unless it can be shown that the EXTERNAL 
FUNCTION command will be executed before the de- 
vice can become READY again. 


TIMING AND INTERLOCKING 


Time allowed for switching transients in the ANTI- 
COINCIDENCE NET and SEQUENCE CONTROL 
is determined by the timing of the following pulses: 


1) NET ENTRY PROBE—MAIN PULSE 0 and 
ANTICOINCIDENCE NET empty. 

2) TRANSFER COMPLETE (READ)—MAIN 
PULSE 2 and EXTERNAL READAMCESiG): 

3) TRANSFER COMPLETE (WRITE)—first 
MAIN PULSE 0 to occur after write transfer in- 
struction is executed. 


4) INTERRUPT probe—MAIN PULSE 0. 


The minimum time allowed is given in the following sit- 
uation: 


1) The ANTICOINCIDENCE NET has an entry 
(READY signal) from a device competing for 
computer time. 

2) An input device with higher priority is interrupt- 
ing the computer. 


The TRANSFER COMPLETE which releases the 
ANTICOINCIDENCE NET is given on MAIN 
PULSE 2 of the READ transfer. The READY signal 
is waiting to establish a new priority. The next interrupt 
probe to occur will determine the earliest possible time 
that an attempt will be made to make use of the newly 
established priority. The next interrupt probe will oc- 
cur on MAIN PULSE 0 of the command following the 
transfer instruction. The minimum time will be 22 usec. 

When the information transfer instruction is executed 
in a service subroutine, the timing must guarantee that 
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Fig. 6—Timing. 


at least one more command can be executed before a 
further interrupt. 

The SEQUENCE CONTROL lockout of INTER- 
RUPT request is relased on MAIN PULSE 2 of the 
READ transfer command. The preceding MAIN 
PULSE 0 probed INTERRUPT REQUEST and found 
none because of the lockout. INTERRUPT REQUEST 
can be generated after MAIN PULSE 2 of the READ 
transfer instruction, but will not be probed to set up an 
INTERRUPT DEMAND until the following MAIN 
PULSE 0. Normal command procurement is therefore 
guaranteed on the MAIN PULSE 6 following execution 
of the READ transfer command. The timing is illus- 
trated in Fig. 6. 


PROGRAMMING 


Program interlocks must be set up to guard against 
several situations which will arise. For output, storage 
space is reserved to stack the results of computation. As 
soon as the results are stacked in this file, an output de- 
vice is assigned to empty it. This is done under program 
control by executing an EXTERNAL FUNCTION 
command turning INTERRUPT on for the output de- 
vice. When the file becomes empty, INTERRUPT 
must be turned off. Turning INTERRUPT on too soon 
or off too late will result in some nonsense in the output. 
An interlock must also be provided which prevents over- 
writing the file with new results before it has been emp- 
tied of previous results. If computation is finished be- 
fore a file has been emptied of results, the computer 
must be kept running until the output job is complete. 
Provisions for saving the contents of the ARITHMET- 
IC REGISTERS and remaking them before reenter- 
ing the interrupted problem must be included in the 
service subroutines. 


The nature of INTERRUPT with randomly-timed 
devices is such that the thin line which has separated 
program difficulties from hardware difficulties is now 
invisible. A problem with a program error may be suc- 
cessfully run many times before the peculiar circum- 
stances which produce a failure show up. When they do 
show up there is little or no evidence to be analyzed. 


CONCLUSION 


Communication of external input and output devices 
with the Remington Rand 1103 computer has been ac- 
complished on an INTERRUPT basis. The devices 
compete for computer time whenever they are ready. 
Service subroutines stored in the memory of the com- 
puter are referenced by the devices themselves when 
INTERRUPT occurs. 

A device may not solicit computer time unless it is 
ready for an immediate transfer of information. The 
computer is never required to stop and wait for some 
mechanical work to be accomplished. The devices in the 
system are required to do all preparatory work “on 
their own time.” 

Development of a system of communications based 
on INTERRUPT, with program controlled use of in- 
ternal storage as a buffer, has eliminated the undesir- 
able effects of entirely sequential programming of in- 
put, computing, and output. Multiple device INTER- 
RUPT permits a mixture of input, output, and com- 
puting, at optimum speeds for all. 

The experience thus far with interrupt has been very 
favorable. In the first use of the feature, which was on 
an engineering problem with a large output, effective 
computer speed was very nearly doubled. A somewhat 
better improvement is expected in data processing prob- 
lems when the full ration of equipment is installed. 
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Questions and Discussion* 


Standards, 


The Symposium on Computers in Simulation, Data Reduction, and Control held at the 1957 
National Convention of the Institute of Radio Engineers, was sponsored by the Professional 
Group on Electronic Computers and presented in the Starlight Roof of the Waldorf-Astoria on 
March 21, 1957. Five invited papers were followed by questions and discussion between the 
audience and a panel that included the authors, except E. L. Braun who could not be present, 
Harold K. Skramstad of the National Bureau of Standards, Morris Rubinoff of the University of 
Pennsylvania, Eric Wolf of Lincoln Laboratory, Walter F. Wakefield of the General Electric 
Company, and the session organizer and chairman, R. D. Elbourn of the National Bureau of 


Chairman Elbourn: In choosing the subject of this 
symposium we took our cue from a remark of Jay For- 
rester, who pointed out that between 1945 and 1950 we 
were learning how to build the first computers; from 
1950 to 1955 we were primarily concerned with solving 
mathematical and engineering problems; and in the 
period that we are in now—1955 to 1960—the principal 
interest is in business applications and data processing. 
He predicted that from 1960 to 1965 our principal in- 
terest would be in automatic control applications of 
computing machinery. 

We had to use a long and awkward title for this ses- 
sion because there did not seem to be any one word 
which described this group of applications, although 
they have very recognizable characteristics in common. 
Most notable, perhaps, is the fact that the computer is 
no longer the prima donna at whose beck and call every- 
thing else starts and stops. Instead, we now have the 
computer dancing to the tune called by some other part 
of the system. It may have to be running in real time; 
it may be interrupted; it may have to respond to an 
unexpected need of the process control in a plant. Usu- 
ally we find both analog and digital information in the 
same system, and thus the necessity of converting from 
analog to digital and from digital to analog; and often 
we find the need for communication between men and 
machines. 

Question from the Audience: On what experience do 
you base the statement of one failure per hour for pres- 
ent-day digital computers? 

Chairman Elbourn: If I recall, the figure was one 
failure per day, and an average of about 23 hours and 
some minutes has been achieved by the Whirlwind. 
Perhaps the author had some others in mind. 

Mr. Gunning: Whirlwind is probably the best exam- 
ple. I tried to get some examples from people who are 
using commercial machines, and I got the general im- 
pression that a few hours is more appropriate for the 
regularly available machines of the large class. I am not 
talking about the smaller machines. 

Question: What do you call a small machine? 


* Manuscript received by the PGEC, February 14, 1958, 


Gunning: I will say what I mean by a large machine 
and exclude all the others. I mean machines of the IBM 
700 series or the Whirlwind—machines of that type. 

Dr. Rubinoff: I think the comment made by one of 
the speakers, Mr. Bauer, was an interesting one. He 
raises what I consider an important question, namely, 
the optimum interconnection of economical units of 
different types of computers, where I guess by “eco- 
nomical units” I mean the small computers; I will ex- 
clude the large ones. In particular, I was thinking about 
one speaker, who favors the use of a DDA and perhaps 
even the use of a digital computer and DDA intercon- 
nection. Another speaker favors the digital and the ana- 
log interconnected through a converter. 

The obvious question is, what would we win if we 
tied together a digital computer (a relatively inexpen- 
sive one), a DDA, and an analog computer, 7.e., all three 
relatively inexpensive computers? Would there be any 
real advantage from the standpoint of the amount of 
work per dollar? Would there be any inherent disad- 
vantages? 

Dr. Bauer: With the present state of the digital dif- 
ferential analyzer, I don’t feel that it is a good com- 
petitor with the more conventional analog computer in 
an analog-digital hybrid system. However, I under- 
stand that at least one manufacturer is coming out with 
a digital differential analyzer that has a much higher 
integrating frequency, by an order of magnitude, as I 
recall. In this case I think the prospect of using conven- 
tional general-purpose electronic computer techniques 
with differential computer analyzer techniques seems 
very excellent. 

I don’t know why anyone would ever want to put the 
conventional analog computer in with the DDA and the 
GP computers, provided that, as I said, the digital 
differential analyzer were of the more modern, ad- 
vanced type which, unfortunately, is not on the market 
right now. 

Mr. Braun:! I would like to return to a point men- 
tioned in my paper. It was stated there that absolute 


‘Since Mr. Braun was not present, his contributions to the dis- 
cussion were added during editing of the transcript. 
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computation was appropriate in some areas and incre- 
mental computation was better in others. There are 
many ways of combining absolute and incremental com- 
puting elements into a computing system. As yet, no 
proof has been offered of what is generally the best com- 
bination, or even that there is a “best” combination. The 
problem hinges on two important considerations: the 
division between relatively high- and low-speed comput- 
ing requirements in a system, and the components avail- 
able for mechanizing the computer. Several absolute- 
incremental hybrid types of computers have been pro- 
posed and/or are being built. However, not enough 
experience has been gained to allow one to state with 
certainty the nature of the “best” computer even for a 
given class of applications. If one could determine this 
for elements currently available, the answer would be 
subject to change as new elements appeared. The follow- 
ing situation could occur. The appearance of certain new 
elements might dictate that one abandon a logical de- 
sign of type “A” in favor of one of type “B.” Later, the 
appearance of other elements could cause one to go toa 
design of type “C” or possibly to revert to a design of 
type “A.” 

It may be expected that such oscillations in design 
philosophy may continue for some time to come. 

I agree with Dr. Bauer’s questioning of the advisabil- 
ity of combining a DDA and a GP digital computer with 
an analog computer. In all probability, DDA’s can be 
developed that can cause the obsolescence of present 
analog computers for such applications as the simula- 
tion of high-speed aircraft and missile performance. 

Rubinoff: There was a question in the middle of Dr. 
Bauer’s answer, and I would like to say one word on 
that. He said that he didn’t know why an analog com- 
puter might want to be tied to a digital computer if it 
has a DDA. In the design of a control system, you may 
have to include analog components such as _ servo- 
mechanisms, in which case you might just as well in- 
clude some analog computation. 

We ought to be a little clearer as to what we mean by 
an analog computer. It seems to me that a complete 
control system has to have such things as motors and 
sensory equipment. The sensory equipment might well 
be digital, although I suspect that it would be analog. | 
guess you can make digital motors, but I don’t know 
of any effort to develop them at this time. I feel, there- 
fore, that all three do tie in together. 

I would also like to comment on the other speaker who 
mentioned the DDA. I got the impression from speaking 
to someone recently that there is a real advantage in 
including a DDA even in a digital computer of high 
speed. The advantage of the DDA is that once the 
necessary derivatives are available, the actual integra- 
tion process of a slow DDA is comparable in speed with 
the integration time of a high-speed digital computer, 
because of the way the information is processed in the 
DDA to obtain optimum efficiency of the computing 
units. 
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Braun: Dr. Rubinoff has just raised two questions 
and has also presented speculations on the answers. 
With the first speculation I disagree. With the second I 
concur, and would like to add some substantiating com- 
ments. 

First, it does not follow that because you happen to 
have analog sensors and effectors you should also have 
some analog computation. Analog computation would 
have to be justified by taking other criteria into con- 
sideration, e.g., such things as reliability, cost, required 
speed of computation, precision, etc. 

Second, it is true that there are advantages to adding 
a DDA toa GP computer. One of the points that I have 
tried to stress is that the DDA is a very efficient means 
of computing continuous functions. Its principal limita- 
tion in existing embodiments is its speed. However, it 
is possible at the present time to increase this speed by 
at least two orders of magnitude (from approximately 
100 to 10,000 iterations per second). This will extend 
the DDA’s range of usefulness to the point where it will 
not only be superior to the GP machine in more areas, 
but will also be fast enough to replace analog computers 
in many high-speed applications. 

Mr. Glassman (Burroughs): The integration per- 
formed in a DDA is not real-time in the sense that it is 
in an analog computer, is it? 

Rubinoff: This wasn’t implied at all. It was a question 
of having the information arrive at the right point at 
the right time so that there is no wasted time in trying 
to arrange to get the information when you want it, as 
you do in a conventional digital computer. There are no 
instructions, for example, no housekeeping, no adjust- 
ing, none of that. Everything goes in a smooth, rapid 
manner. 

Bauer: Real-time only means the effective speed at 
which computations are performed, not the method. 
Serial computation can be real-time, just as well as 
parallel computation. What is not true is that in the 
serial operation every variable in the machine, or its 
analog, has exactly the value it should have at every 
instant. 

Braun: Integration in a DDA is performed in a dif- 
ferent manner than in an analog machine, and it is not 
in “real time” in the sense that the time taken for the 
process is determined by the “real” response time of an 
electronic circuit. However, as generally used in comput- 
ing circles, the term “real time” refers not to the manner 
but rather to the rate of a computation. What is usually 
meant when it is said that a computer operates in “real 
time” is that it operates sufficiently fast to keep up with 
a physical process being followed or controlled. Since 
different physical processes take place at different rates, 
the term “real time” is relative. Its use always implies 
reference to the rate at which a specific physical process 
occurs. 

Glassman: I have one more question with regard to 
the DDA. When taking variables in from the outside 
world, you would like to have them, as I understand, 
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come into the DDA in an incremental manner. How do 
you accomplish this? 

Braun: There was insufficient time to discuss the 
input-output problem in the oral presentation, though 
it is discussed in the written presentation. Therefore, I 
will only comment briefly here on how input variables 
are treated in a DDA. 

It is not necessary to have external inputs enter a 
DDA in an incremental manner. In fact, inputs are 
usually entered very much the same as in a GP machine. 
Special registers, termed servo registers, are available 
in a DDA, and they can be used to generate a unitary 
weighted pulse train from absolute binary input data. 
The output of the servo register is the signum function 
of its contents. This output is also fed back to the input 
of the register with a sign reversal. Thus, whether the 
contents of the register is a positive or negative number, 
the register tends to drive its contents to zero. This 
occurs at the rate of a single increment per sampling 
cycle, and in the process a unitary weighted pulse 
stream is generated which more or less follows, on the 
average, the variations of the input variable. 

Dr. Skramstad: One way to get variables into the 
DDA from an outside analog source would be to convert 
the analog input to digital form and transfer directly 
to a DDA register. Another way of doing this would be 
to produce an analog voltage proportional to the con- 
tents of a register, compare this with the analog voltage 
to be inserted, and use the difference to cause incre- 
ments to be added or subtracted from the register until 
balance is obtained. It could be arranged so that when- 
ever the rate of change of the variable exceeds the 
counting rate that the DDA can follow, the machine 
counts by 2’s or 4’s or at whatever speed is necessary in 
order to follow. 

Mr. B. K. Smith (Beckman Instruments): The usual 
way of getting an analog input into a DDA is to per- 
form a digital-to-analog conversion on a number repre- 
senting that input, and then “update” that number ac- 
cording to the polarity of the output of an analog com- 
parator. 

I would like to ask whether anyone in the audience or 
panel knows of an incremental computer wherein the 
sampling rate is made a function of the rate of change; 
that is, where an input is observed to be frequently 
changing, but at much longer intervals in the normal 
scan sequence. It seems to me that a speed increase of, 
perhaps, a hundred-fold could be obtained in systems 
where only a few inputs vary at any given time. 

Chairman Elbourn: Apparently no one knows of such 
a device. 

Bauer: Some of us on the panel have referred to the 
so-called interrupt feature, and we mentioned that we 
are going to see a great deal more of this. As you know, 
it is a standard feature on the 1103 computer, and is 
also available on the 704 computer. On neither of these, 
nor on any computer that I know of, are multiple inter- 
rupts available. 
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Mr. Turner whispers in my ear that at N.A.C.A. 
they can distinguish interrupts, and that is the point I 
wish to make: you must not only have the interrupts, — 
but you will need computer distinguishable interrupts, 
too, so that if three or four or more outside equipments 
happen to be interrupting, the computer will be able 
to distinguish among them. 

The other point about the interrupt is this: we are 
seeing for the first time in digital computers the com- 
puter working on two problems simultaneously—on a 
time-sharing basis, of course. For instance, one idea we 
have in our shop is to supply, let us say, 100 cells for the 
analog computer to use just as it wishes. Then the ana- 
log computer, by means of the interrupt and the analog- 
digital converting equipment can steal time without the 
digital computer’s knowing about it, so to speak. That 
is, the digital computer can be working on an entirely 
different problem, be interrupted at times by the analog 
computer, and therefore essentially do two problems 
simultaneously. 

Mr. Turner: [ would like to extend that comment a 
little. There is a point in connection with interrupt 
which must not be forgotten in making a computer solve 
two problems at one time: the base problem must be 
guaranteed to be free of errors. Actually, our second at- 
tempt at interrupt was to completely override a standby 
problem with provisions for returning to this problem 
later. 

It gave rise to a rather amusing situation. One of our 
problems involved a computer-output typewriter. The 
man who had the control over the interrupt had a push- 
button. Whenever he pushed it, it would automatically 
interrupt the data processing problem. He managed one 
time to interrupt one line of typing four times. 

Bauer: I don’t know if it is proper for panel members: 
to ask questions of other panel members, but I have a 
question to ask of Professor Susskind. He mentioned the 
fact that it was impossible, or at least very difficult, to 
put feedback resulting from measuring the effect of the 
cutting device back into the computer and thereby re- 
adjust the position of the cutting device. I wonder 
whether the loop can’t be extended to make an ap- 
proximate cut, go back over the cut, and measure and 
put the resulting measurement back into the computer. 
Following this there would be some calculations, and 
then another cut, and so forth. Thus one would close 
the loop in a longer fashion rather than directly as 
would first be conceived. This would involve more ex- 
tensive computer operation, I am sure, but I wonder if 
it is feasible. 

Mr. Susskind: I don’t see any reason why this could 
not be done. I suspect the only limitation here is not 
technology but economics. 

Mr. Wakefield: May I ask a question? In a system 
which includes general-purpose equipment and DDA’s 
and analog equipment, how do you check on the output 
when it is required to be continuous? In a general- 
purpose equipment, you can stop the computation to 
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check the results. How do you do this in DDA or analog 
equipment? 

Bauer: I am not quite sure whether I understand the 
question. I guess you check out the whole system in the 
same way as you would check out an analog computer 
-where you have the same difficulty of stopping in the 
middle of the problem; that is, run it all the way 
through and see if the answers are reasonable. Of course, 
there are many static tests that can be placed on the 
equipment before it is tied together and run as a system. 

Susskind: May I put in a comment here which 
slightly contradicts one of the other speakers? In ma- 
chine-tool control one usually speaks of a data processing 
system which is small enough to be made at least as 
reliable as the basic machine tool. I am not at all worried 
about putting any amount of electronics on the ma- 
chine-tool provided I can justify it economically, and I 
am not worried about downtime on the machine-tool or 
spoiling parts. We have probably cut in our lab, strictly 
on an experimental basis, 400 or 500 parts, and we 
spoiled only two of them. One was spoiled by the ma- 
chine, the other by the operator. 

Glassman: | would like to ask the speaker from Beck- 
man a question. | think you indicated that the digital 
computer in a control loop could be checked by the ana- 
log servo controller. If the digital computer made an 
error, there would be some method to check by the 
analog control. I don’t see how you can do that. Is that 
what you said or did I misunderstand you? 

Gunning: That is what I said, and by that I meant 
to imply something along the lines of incremental data 
transmission—the sort of thing that you and the panel 
have been discussing. For example, assume that the 
set point of the servoregulator is allowed to move only a 
small increment each time the digital computer gives it 
an instruction. If the digital computer were to request a 
change that is very much larger than the maximum 
value that you know would be “reasonable” for the proc- 
ess, then the analog regulator would take the small 
step that is allowed to take and signal back to the digital 
computer that an unreasonable request has been made. 

Mr. Weiss (Waldorf Instrument Co.): I have a ques- 
tion directed toward the gentleman who presented the 
paper in reference to the 15-channel A and D converter. 
Was that converter a separate unit from the 1103, or 
was it strictly fed into the 1103? 

The second part of my question is, was it strictly a 
15-channel input which is electro-input, or was it a com- 
bination electro-mechanical? 

Bauer: The converter equipment is a tailor-made 
piece of equipment by the Epsco Company, of Cam- 
bridge, Mass. It is of very high speed, it is very accurate, 
and it is very expensive. It is completely electronic with 
speeds of 90 microseconds per analog-to-digital conver- 
sion, and it is made up of somewhat standard com- 
ponents and techniques, however, with some extra 
equipment, such as an oscillator for clock pulses and 
control mechanisms, which give it flexibility in use. It 
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has fifteen channels of analog-to-digital conversion and 
fifteen channels of digital-to-analog conversion, giving 
it a 30-channel capacity. I am not sure I have answered 
all of your questions. 

Weiss: The other half of my question was this. 
Among the input channels, do you have a system by 
which you feed in mechanical data also, as well as elec- 
trical data? 

Bauer: The converter will, of course, read into the 
digital computer anything that can be presented on the 
plus or minus 100 volts scale of the analog computer. 
This means, for example, that we can take FM tape and 
play it through the proper equipment into the converter 
and, in turn, into the digital computer. 

Glassman: In connection with that question, I don’t 
understand why you are using 1-millivolt accuracy 
when 100 volts is what you get out of the converter. 
Why are you using such a converter as that? 

Bauer: | believe that the 1-millivolt accuracy comes 
from the desire to measure and convert precisely those 
quantities which are small, as well as the quantities 
which have a large voltage value. 

Question: Is it 1-millivolt accuracy or resolution? 

Bauer: We refer to it as accuracy. I am not sure that 
I know the exact difference between them; I am sure 
there are various interpretations of those two words. 

Rubinoff: I would like to say a little more about the 
original question that Dr. Bauer raised. I hesitate to 
start on it because it is a topic I have done quite a bit of 
thinking on and there is a danger of my overstepping 
available time. But I would like to give you some indi- 
cations of the direction of the thinking. 

We have been talking today about digital computers 
that are capable of doing more than one problem be- 
cause they have an interrupt feature, about digital 
computers which are tied in with analog computers, and 
about digital computers which are tied in with DDA’s. 
The question I would like to raise is, what is the proper 
size of module that we ought to be talking about? We 
seemed to be talking about the module at the level of a 
complete computer. And here we are talking again from 
the standpoint of a complete converter tied to a com- 
plete computer tied to another complete computer, one 
computer being digital, the other analog. 

What little bit of thinking has been done so far indi- 
cates that maybe the size of module is wrong. Maybe 
because everybody said ten years ago that it is very 
difficult to build a digital computer, the feeling has 
grown that a complete computer is a basic module. 
Perhaps we have forgotten that that was something 
that was important ten years ago but not necessarily 
important today. 

As I said, I don’t want to go into a long oration, but I 
do want to put across this last thought, that perhaps a 
module ought to be just a storage register or a group of 
storage registers and another module just a counter, an- 
other module just an accumulator. Perhaps the whole 
should be set up in such a way that all of these are inter- 
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connected in such a manner as to make optimum use of 
the equipment. 

This again ties in to the difference between a DDA 
and a digital computer. A DDA works much more effi- 
ciently on the job for which it was built because it is 
set up to do that job without a lot of housekeeping and 
getting prepared and shifting around of data. Even 
though it is inherently a slower, less costly machine, it 
competes favorably in many cases with a high-speed, 
expensive digital computer. 

The question is: Can we take smaller modules and 
put them together in such a way that we keep these 
modules active more nearly 100 per cent of the time, 
rather than called into play perhaps 10 or 15 per cent of 
the time as is the current situation? 

Mr. Flesch (Federal Telecommunication Laborato- 
ries): I would just like to remark on this last comment, 
that it seems that this would replace the program by a 
logical design for each problem. 

Rubinoff: Not at all. As a matter of fact, the plan that 
we are working on is to eliminate altogether the man 
who schedules the problem for the computer. Why not 
just put in the priority information at the outset and 
let the whole system put itself together in any way it 
pleases to handle the problems as they come, in accord- 
ance with the priorities that are put on them? At pres- 
ent this sounds like a real dream, as digital computers 
may have sounded fifteen years ago, but, upon first 
examination at any rate, it looks as though this is quite 
feasible. There are, in fact, some really attractive fea- 
tures, and it is not nearly as expensive as one might at 
first think, because so much is saved on the digital com- 
puter equipment, the actual computation equipment, 
that you can afford to put it into the priority determina- 
tion, of if you like, into the master control. I don’t think 
it really makes life any more difficult for the program- 
mer, particularly now that we have automatic pro- 
gramming. 

Chairman Elbourn: Dr. Rubinoff, your remarks re- 
mind me a little bit of microprogramming, building the 
digital computer in such a way that the more elemen- 
tary operations of the machine are available, such as a 
simple transfer from Register A to Register B, so that 
either by plugging or programming one can synthesize 
complex computer operations and possibly make a DDA 
out of a general-purpose computer if desired. 

Rubinoff: Correct. Microprogramming is what I have 
in mind, but microprogramming in such a way that the 
system is always using all the equipment, and operating 
on several problems at one time. This apparently is not 
too difficult to mechanize. The problem is still an eco- 
nomic one, and it may be years before such a system 
can be achieved dollar-wise, but conceptually there is 
real hope in this area. 

Braun: When I spoke of the relative advantages of 
absolute and incremental computation, I did not mean 
to imply that one should necessarily interconnect a con- 
ventional GP computer and a DDA. For some applica- 
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tions this might be adequate. However, in general, one 
would wish to look at the combination of absolute and 
incremental computation at the level of small modules 
in order to see logical relations that might allow simpli- 
fication of equipment requirements. Storage registers, 
counters, and accumulators, as Dr. Rubinoff suggests, 
are very reasonable choices for modules. Such devices 
are basic for all common types of digital computation, 
that is, for absolute arithmetic units, for digital inter- 
grators, and for binary rate multipliers, as well as for 
control circuitry in general. However, the lower the 
level of integration, the greater the probability that one 
will have a computer that is relatively more specialized 
in function. 

In regard to Dr. Rubinoff’s remark about the speed 
of a DDA, I would like to suggest that the DDA is not 
inherently slower than a GP machine in the solution of 
differential equations. Current DDA’s, operating in a 
completely serial mode with a magnetic drum store, may 
be slower than GP machines having considerable paral- 
lel circuitry and high-speed random access storage, but 
this is hardly a valid comparison. However, as Dr. Ru- 
binoff points out, even these DDA’s compete favorably 
with GP machines in many instances. 

A few general comments may be in order on the sub- 
ject of GP computers. There has been, in the past, a 
great emphasis on this type of machine. In the early 
stages of a new development, the general purpose ap- 
proach is justifiable on the following grounds. First, the 
requirements of potential users usually are not clearly 
defined, and a GP machine provides a customer with 
more assurance that, even if his requirements change, he 
will not be left with a piece of equipment for which he 
will have relatively little use. Second, development costs 
can be spread out over a wider area. However, one 
does sacrifice efficiency for general utility. It is not un- 
reasonable to expect that in the future, as large classes 
of special applications become more clearly defined, 
it will be economical to develop a number of types of 
special purpose computers, each of which is most efficient 
in the solution of a particular class of problems. 

~ Mr. Wolf: I just want to comment on Mr. Gunning’s 

remark on the reliability of a digital computer in a real- 
time industrial control system. While it may be un- 
economical for that purpose, we, at Sage, use a very 
large number of both special-purpose and general-pur- 
pose interconnected digital computers, and, by means 
of redundancy—on a very large scale to be sure—we 
hope to have an over-all system of reliability which is 
characterized more by the six months than by one day. 
While this probably would not be economical for an in- 
dustrial process right now, it does show that it is at 
least technically feasible to have a such control system 
even for systems that are very likely to be more com- 
plicated than any industrial system. 

Chairman Elbourn: With reference to the reliability 
required in industrial control systems, I was a bit shaken 
one day by some visitors to the Bureau of Standards 
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from Imperial Chemical Industries, Ltd., of Great 
Britain. They wanted to know if any work is being done 
on digital computer circuitry that could be used for 
chemical industry process control. And they explained 
that they would not like to have a major breakdown any 
oftener than once in five years. Until then I thought we 
had already achieved rather good reliability; the ques- 
tion showed me how far we have yet to go. To do that 
kind of thing, we shall have to introduce the kind of 
system that Bill Gunning mentioned, namely, a hier- 
archy of control loops with perhaps unsophisticated but 
very reliable things in the lowest control loops, and the 
more sophisticated computing equipment superimposed 
in higher loops, in such a way that if something goes 
wrong there, it only affects the system about as much as 
one martini affects a human being. 

Solid-state components certainly should let us make 
some of these lower control loops as reliable as the 


Logical Machine Design: 


Summary—This paper presents a selection of current biblio- 
graphic references in the field of the logical design of machines. 
Selections were made with the objective of being useful and stimula- 
tory to the wide variety of scientific persons who have some interest 
in this subject. An essential feature of the report is the extensive 
index of significant title words. Each indexed word is given in the 
complete context of the title in which it occurs. 


INTRODUCTION 


N THE PREFACE to “Elements of Symbolic 
| Logic,” Professor Hans Reichenbach says, “. . . the 
logistic approach to philosophy is not bound to a 
certain type of mind, or of milieu or educational system, 
but represents a most successful clarification of ideas, 
by which all forms of scientific pursuit will profit.” The 
design of finite automata is a scientific pursuit in which 
this observation appears to be sharply appropriate. 
The widely increasing emphasis now being given to re- 
search in mathematical logic as applied to computing 
machine networks and components is convincing evi- 
dence. The originality and usefulness of results already 
obtained should provide ample encouragement for con- 
tinued endeavor along these lines. 
The title of this report can be read in two ways, and 
it should be. The area of interest includes not only the 
use of mathematical logic in machine design (as con- 


* Manuscript received by the PGEC, November 8, 1957. 
+ Electronic Components Lab., Wright Air Dev. Center, Wright- 
Patterson AF Base, Dayton, Ohio. 
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pneumatic equipment that is now being used. 

Mr. Lipple (Signal Corps Engineering Laboratories) : 
I was just wondering whether, as a result of some of 
these last comments, it wouldn’t be well to point out 
that where at one time you had general-purpose ma- 
chines, nowadays you tend to get different types of gen- 
eral-purpose machines; maybe the concept of a general- 
purpose machine, or the name, has outlived its useful- 
ness. Maybe we should talk about multipurpose ma- 
chines and realize that for different kinds of jobs we may 
have need in the future for different kinds of multipur- 
pose machines. 

A multipurpose machine for simulation, or for some 
of the real-time problems, might include many of the 
concepts that were discussed this morning. 

Chairman Elbourn: We have come to the end of the 
time for this session. We certainly appreciate your par- 
ticipation. 


A Selected Bibliography” 


DOUGLAS B. NETHERWOOD{ 


trasted with empirical, cut-and-try techniques) but also 
the use of logic machines to help solve difficult or 
lengthy problems in logic and even to generate new con- 
cepts and formulations. See, for example, [24].! 


SCOPE 


The reader may be surprised to find that Boole, 
despite the profound importance of his work (in mid- 
19th century), does not appear in this bibliography. 
Nor do Jevons or De Morgan. The investigations of such 
early writers have been thoroughly revised and assimi- 
lated into more modern works of wider application. This 
report does not offer a historical compilation but lists 1) 
informative papers on what has been done and is being 
done, and 2) stimulatory papers on what needs to be 
done and what might be possible. 

Because the discoveries and techniques of logic have 
proved uniquely effective in the design and analysis of 
discrete-state machines, as opposed to analog or con- 
tinuous-variable types, nearly all titles listed are per- 
tinent to that subject. However, the area of interest has 
been somewhat broadened beyond the limits of the 
central subject alone (the direct use of logic in digital 
and finite-state problems). An attempt has been made 
to assemble a spectrum of material in which almost any 
investigator, regardless of specialized approach, can find 


1 Numbers in brackets throughout this report refer to biblio- 
graphic listings. 
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topics of interest. A particular effort has been made to 
anticipate and accommodate such individuals as, for 
example, an engineer who might be scanning the titles 
in this field and ask: “How have some of these abstract 
ideas been reduced to practice? What are some new 
phenomena and devices that offer further realizations?” 
On the other hand, a mathematician might say: “Given 
these system requirements and unsolved problems, what 
are some theoretical developments and studies of pos- 
sible pertinence to their solution?” 

The range of technical difficulty of the papers listed 
is from rigorous mathematical expositions to popular 
articles and news items, although these latter have been 
confined to a few notices of unusual promise. A few 
papers have been included from related areas such as 
machine philosophy, information theory, game theory, 
neurological analogies, and the mutual influence be- 
tween programming and coding and computer logic and 
design. 

Accessibility and readability of papers have been 
given some consideration, but the more significant titles 
have been listed regardless of probable availability or 
language of publication. An effort has been made to 
minimize the effect of personal bias on the part of the 
compiler. However, his opinion that truth-function 
theory and solid-state physics devices are especially 
promising is undoubtedly reflected in the selections. 


ORGANIZATION OF THE BIBLIOGRAPHY AND INDEX 


Some thought was given to separating the biblio- 
graphic listings into categories. At least half of the titles 
could be so arranged, but the remainder would so over- 
lap the boundaries that it would be necessary to make 
either multiple entries or arbitrary assignments. There- 
fore, a different procedure was selected. 

Authors are given in alphabetical order in the Bibliog- 
raphy; but the Index is the key to this entire work. All 
significant title words appear as entries in the Index, 
together with the remainder of the title in each case. 
Where the significant word is not first in the title, a slant 
line precedes it. By this means, it is hoped that desired 
papers and related ones can be located with ease and 
accuracy. The reasoning for this arrangement is that 
most authors have selected their titles with care, and 
that no better method could be chosen than to give the 
salient words they have settled upon. The fairly common 
practice of listing single words in an index, followed by 
all relevant reference numbers, was rejected because it 
does not show the reader the context and may force him 
to search through dozens of unwanted references to find 
the one desired. 

The use of the Index is straightforward. One has only 
to locate known or probable title words and note the 
associated reference numbers. In searching for titles 
involving phrases such as “Boolean algebra,” “system 
design,” or “logical matrix,” it is recommended that 
both words be checked. 
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READING SUGGESTIONS 


For readers who may be unfamiliar with the litera- 
ture in this field and who would like to look at some in- 
troductory material, it is suggested that some selections 
be taken from among the following references: [30], 
[51], [62], 1107], 1145]; [195], 123A. [243 hizo 
[301], [364], [384], [456]. 

Periodicals in which pertinent articles are frequently 
to be found, arranged in approximate sequence of tech- 
nical difficulty, are as follows: 


Electronic Design 

Automatic Control 

Computers and Automation 

Instruments and Automation 

Control Engineering 

IRE TRANSACTIONS ON ELECTRONIC COMPUTERS 
Journal of the Association for Computing Machinery 
Journal of Computing Systems 

Journal of Symbolic Logic. 


Other Bibliographies 


There are a number of excellent bibliographies in re- 
lated fields. Some of the more important of these are: 


1) Aiken, H. H., and Staff of the Computation Laboratory, Har- 
vard University. “Bibliography of numerical analysis” in “A Manual 
of Operation for the Automatic Sequence Controlled Calculator,” 
Harvard University Press, Cambridge, Mass., pp. 338-396; 1946. 

2) Army Library. “Data Processing Operations and Computers 
for Management,” Special Bibliography No. 10, Dept. of the Army; 
October 4, 1956. 

3) Beach, A. F., ef al. “Bibliography on the Use of IBM Machines 
in Science, Statistics, and Education,” IBM Corp., New York, N. Y.; 
1954. 

4) Bureau of Labor Statistics. “Automatic Technology and its 
Implications; a Selected Annotated Bibliography,” Dept. of Labor 
Bulletin No. 1198; August, 1956. 

5) Chapin, N. “Publications for business on automatic com- 
puters: a basic listing,” and “ . . . asupplemental listing,” Computers 
and Automation, vol. 4, September, 1955; vol. 5, February, 1956. 

6) Chessin, P. L. “A bibliography on noise,” IRE TRAns. ON 
INFORMATION THEORY, vol. IT-1, pp. 15-31; September, 1955. 

7) Church, A. “A bibliography of symbolic logic” (parts I and 
II), J. Symbolic Logic, vol. 1 (no. 4, entire) and vol. 3, pp. 178-212 
(available in separate reprint); December, 1936, and December, 1938. 

8) Clark, F. A. “Computers; Including Fire Control Computers, 
Flight Calculators, Mathematical Computers, Navigation Com- 
puters: A Report Bibliography,” ASTIA Reference Center, Library 
of Congress, Report No. ARC-854. 

9) Comstock, M. “Bibliography on Machine Computation 1945— 
1954,” Brookhaven National Lab., ASTIA No. AD 62539; May, 
1955. 

10) Culbertson, J. T. “Robots and automata: a short history. 
bibliography,” Computers and Automation, vol. 6, pp. 20-24, 38; 
April, 1957. 

11) Hollander, G. L. “Bibliography on Data Storage and Record- 
ing,” Servomechanisms Lab., M.I.T., Cambridge, Mass., ASTIA 
No. AD 15789; 1953. Also in Trans. AIEE, vol. 73, pp. 49-58; 
March, 1954. 

12) McCann, G. D. “Bibliography of Articles on Computing 
Machines and their Applications,” California Inst. Tech., Pasadena, 
Calif.; October, 1949. 

13) Miles, J. G. “Bibliography of magnetic amplifier devices and 
the saturable reactor art,” Trans. ATEE, vol. 70, pt. 2, pp. 2104— 
ZA OSE 

14) Nottingham, W. B. “Bibliography on Physical Electronics,” 
Res. Lab. of Electronics, M.I.T., Cambridge, Mass.; 1954. 
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15) Parke, N. G., III, “Guide to the Literature of Mathematics 
and Physics Including Related Works on Engineering Science,” 
McGraw-Hill Book Co., Inc., New York, N. Y.; 1947. 

16) Serrell, R., and Rupprecht, K. J. “Bibliography on Electronic 
Computing,” (3rd ed. rev.), RCA Labs. Div., Radio Corp. of Amer- 
ica, Princeton, N. J.; July, 1953. 

17) Stumpers, F. L. “A Bibliography on Information Theory,” 
Report, Res. Lab. of Electronics, M.I.T., Cambridge, Mass.; Feb- 
ruary, 1953. Reprinted’ by IRE-PGIT, November, 1953. “... Sup- 
plement,” IRE Trans. ON INFORMATION THEORY, vol. IT-1, pp. 31- 
47; September, 1955. 

18) Wade, W. F., and Kemler, E. N. “Automatic Control 
Bibliography,” Summary Reports, Spring Park, Minn.; 1955. 
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Unit-Distance Error-Checking 
~Codes* 


An interesting set of binary codes can be 
generated which combines the properties of 
error checking, such as single-error detection, 
and unit distance, which means that suc- 
cessive code elements (-bit numbers) of the 
code sequence differ in only one bit. For 
example, the code sequence for n=4 


0000 
0001 
0011 
0111 
1111 
1110 
1100 
1000 

(0000) 


is of this type, in the sense that each single 
binary error in an element of the sequence 
results in either 1) the next element, 2) the 
preceding element, or 3) an element which 
does not appear in the sequence at all, and 
which may therefore indicate a “detected” 
error. 

In this note some possible applications 
and a few properties of these codes are dis- 
cussed, and several as yet unresolved prob- 
lems and conjectures are offered for their 
possible interest. 

These codes are potentially useful in two 
areas. First, they provide a type of error- 
detection code for digitalized analog data, 
in the sense that the most likely binary 
errors are either detected or result in neg- 
ligible analog errors, and the analog-to- 
digital conversion may be accomplished with 
no commutation error. Second, in asyn- 
chronous digital networks without race con- 
ditions, they provide a code for counting 
which enables most simple circuit faults to 
be detected. Thus, if a failure in the circuitry 
causes a flip-flop to fail to change state 
when it should, or to change state when it 
shouldn’t (or a relay to open or close in- 
appropriately), then either the count will be 
in error by one, or, what is more likely for 
large m, the counter will advance to a for- 
bidden state and the error arising from the 
fault may be detected. The price paid for 
this feature is an increased amount of logical 
circuitry in the counter. Such an arrange- 
ment is a type of “single-fault-detecting” 
sequential net, and is significant from a re- 
liability standpoint. 

The usual model of an error-checking 
code represents each of the code elements 
(%nXn—1* * * X2X1) as a vertex (node) of the 
n-dimensional unit cube! with coordinate 
directions x1, %2 °° + Xn1, %n. A set of code 
elements with the single-error-detecting 
property is then represented by a set of 
vertices, no two of which are adjacent, 1.e., 


* Received by the PGEC, February 14, 1958. 
1R. W. Hamming, “Error-detecting and error- 
correcting codes,” Bell Sys. Tech. J., vol. 29, pp. 147— 
160; April, 1950. 
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joined by a single edge of the cube. A code 
sequence with the unit-distance property is 
represented by a chain of adjacent vertices.? 
It is of interest in the design of codes having 
properties such as these to determine the 
maximum number of elements (nodes) 
which may be incorporated into the codes 
(cube), and to derive if possible a systematic 
pattern for generating such maximal-length 
codes for arbitrary m. The present problem 
is therefore one of packing the longest pos- 
sible chain into the ” cube with the restric- 
tion that no two nonsuccessive nodes on the 
chain be adjacent. (The interpretation of 
this chain as a unit-radius tube led to the 
term “snake in the box,” or SIB, code.) 
One may consider open chains (segments) or 
closed chains (cycles); however, only closed 
chains are discussed here. 

Let K(n) be the maximum possible 
length of an SIB chain in the » cube. The 
following values have been derived, prin- 
cipally by trial and error: 


n 2S. a oe 16 7 
K(n) 46578 14. 24 > 46. 
The law of growth is not known, but the 
above sample, and a familiarity with the 
process of deriving the code in »+1 dimen- 


sions from the code in dimensions, suggests 
the rule 


Wie ae 


K(n +1) = es = eed 


2K(n) — 4 for n odd 
This recurrence relation leads to 

K(n) = 3[2* — (-1)"] +3 n2>2. 
The following bounds are easily established: 
K(n) > 2n 


Ka) <—— 241 
n—1 


K(n) > A2”/? for some constant A independent 
of n; A > 4. 


The first of these bounds follows from the 
pattern of the above example for »=4. The 
second results from a count of the total 
number of edges in the cube, both on and off 
of the SIB path. The third results from the 
realization that an SIB path of length 
2K(n) can be formed in +2 dimensions 
from two SIB paths each of length K(n), in 
n dimensions,? 7.¢., 


K(n + 2) > 2K(n). 


These bounds are generally rather poor, 
although the third one confirms the ex- 
ponential growth displayed in the above 
conjectured formula for K(m). 

An alternate description of the structure 
of a unit-distance code sequence is provided 
by an ordered list of the coordinate changes 


2H. E. Tompkins, “Unit-distance binary codes for 
two-track commutation,” IRE Trans. oN ELEC- 
TRONIC COMPUTERS, vol. EC-5, p. 139; September, 
1956. 

3 J. Brenner of Stanford Res. Inst. reports a slightly 
better bound: K(n) >B(2.25)"/2. 
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Fig. 1—The principal (solid line) and alternate Soman 
line) natural SIB chains for 7» =5. K(5) = 


required in passing from each element (node) 
to the next. For example, for the above 
n=4 SIB code with elements (x4, x3, x2, x1), 
we have, listing only subscripts: 


n = 4: 12341234. 


The values of K() listed above are achiev- 
able with the following codes: 


SPAR I) Se 
n = 3: 123 123* 
n = 4: 1234 1234* 
1234 1243 
1234 2143 
nm = 5: 1234524 123524* 
1234524 1234254 
1234523 1243253 
n = 6: 123456423453 123456423453* 


= 7: 12345674534214561241534 
21451675234563216346154 


For n<5 the codes listed exhaust the sym- 
metry types, 7.e., all other maximal-length 
SIB codes for 2<5 may be obtained from 
those listed by a reversal of direction along 
the path, a permutation of variables, choice 
of a different starting point, or some com- 
bination of these transformations. The num- 
ber of symmetry types for 7>5 is not known. 
The starred (*) codes may be referred to as 
“natural” codes, in that they have the fol- 
lowing unique special properties: 1) The 
path is antisymmetric; that is, the first 
and second halves of the path, from any 
starting point, follow the same coordinate- 
change sequence. 2) The reverse path can 
be derived from the forward path by merely 
inverting the order of the variables. 3) K(m) 
of the 2"—K(n) non-SIB-path nodes on the 
n cube may be selected (for N>3) for an 
alternate SIB path of the same symmetry 
type (see Fig. 1). 4) Less precisely, the logic 
circuitry required in conjunction with 
flip-flops or relays arranged to count with 
this code sequence seems to be simpler, by 
almost any reasonable criterion, than for the 
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other symmetry types. The equations for 
the code of Fig. 1 are listed in Table I; here 
Z is the next value of z, Y the next value of 
y, etc. We conjecture that at least one such 
“natural” SIB code exists for arbitrary n. 


TABLE I 


EQUATIONS FOR » =5 SIB COUNTER WITH 
PRINCIPAL AND ALTERNATE PATHS 


Z =2lw +2 +4) + wes 

¥Y =9[w + (« © »)) + we © 2) 

X =s[9 +2+(w@ 9)] + v2(w © 9) 
W =wly + (« @2)] +9 @ ») 

Va Bee +9) + wxy 
a@®b=ab+ab 


The SIB codes discussed above may be 
directly generalized to provide for multiple 
error checking, and to p-nary rather than 
binary codes. 

The author is indebted to Bernard 
Elspas of Stanford Research Institute for his 
interest and participation in the study of 
these codes, and also benefited from contact 
with E. N. Gilbert of Bell Telephone Lab- 
oratories and Joel Brenner of Stanford Re- 
search Institute. 

W. H. Kautz 
Stanford Res. Inst. 
Menlo Park, Calif. 


An Optimum Character Recognition 
System Using Decision Functions* 


Issue is taken with this statement from 
the above paper by Chow:! “The above re- 
sults also indicate that the recognition sys- 
tem described in Section II (correlation 
coefficient) is not optimum for the particular 
signal and noise structure as given in ex- 
amples 1 or 2.” 

This statement may be true generally 
but it is not true for a most important case. 
That is, for the case to be cited, correlation 
is an optimum recognition system. 


Consider the following weights and 
probabilities: 
Wi =9 Wi = 0 for all z 
Wi=W alli¥z 
fi=p  allz. 


This is the case when all characters are 
equally likely. It is equally costly to mis- 
identify them, but costs nothing to iden- 
tify them, and rejection is not considered. 

The detector proposed examines the out- 
put 


X;(V) = 2 WijpiF (V/ai), 


* Received by the PGEC, February 24, 1958. 

1C, K. Chow, 1957 IRE WESCON CoNvVENTION 
REcoRD, pt. 4, pp. 121-129. Also, IRE TRANS. ON 
ELECTRONIC COMPUTERS, vol. EC-6, pp. 247-254; 
December, 1957. 
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but 
Wiipi = 9 fori =7 
=Wp fori7j 
X(V) = Wp X, F(V/as) 
ers | 


I 


we( XPV /a) - FV/a)). 


To find the minimum X; simply means 
to find max F(V/a;). : 
Examine F(V/a;) as given in (29)! 


(27) */"Mo4; 


F(V/a;) = 


The hypothesis that the noise is inde- 
pendent means that the covariance matrix, 
\|o<;2||, is diagonal. Then o;;=0 for 77. Call 
oj=0;;. Also, call K =[(27)*/I1o,]-}. 

Now define A; as follows, where A; is the 
likelihood ratio discussed in reference 5 of 
Chow’s paper: 


A; = F(V/a;)/F(V) 


where now 


Vi — aii)? 
F(V/a:) = K exp (-x oo) 
and 
Aqaba eo (- > Vit/2oi). 


Maximizing A; corresponds to maximizing 
F(V/a;) as does maximizing log A;. There- 
fore, examine 

log A; = log F(V/a;) — log F(V) 
V;- aij 2 
Se = Ss i i) 


5 2 
7 20; 


+ D0 V}?/20;?, 
9) 


—log K 


V; — ai; 2 V,2 
tog n= Se ee 
oj 


oe 


2 2a Ve + dase 


pe : 
a IE 

=) V?/o7 a 220 2:V;/o7 

Dy Cerca) Vaars, 
log A; = 2X (aij Vj/ox?) — 5 2 (a;;?/0;?) . 


2 


The first term on the right side is the scaled 
correlation coefficient; the second term is a 
bias term correction which corresponds to 
the power content of the reference wave- 
form. 

Thus, the correlation coefficient is an 
optimum system if probabilities and weights 
of the system are equal, as might be encoun- 
tered in actual practice, and if the bias term 
for waveform power is considered. 

IvAN FLORES 
Remington Rand Univac 
South Norwalk, Conn. 
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A Note on the Reduction of Switch- 
ing Functions* 


The purpose of this note is to show how 
to reduce switching functions of high num- 
bers of variables while taking advantage of 
redundancies. 

Any switching function of 1 variables 
can be expressed in canonical form as 


M 
+5 am) = Do gibi (M = 2"—1) (1) 


i=0 


f(m, %2, * 


where the p,;’s are the complete products 
(minterms) and the g;’s are binary coefficients 
peculiar to a given switching function. Eq. 
(1) can also be written in terms of complete 
products g; of any group of four of the n 
variables and complete products v; of the 
remaining 2 —4 variables in the form 


15 


R 
+, Xn) =D Dy rested; 


7=0 k=0 


(Ri= Det — 1) ee) 
= a 6;(k) qi 


(heise 0, bY 2,= RR): (3) 


In the language of the hypercube geometri- 
cal model! a switching function can be re- 
duced if and only if there are at least two 
adjacent vertices on the associated hyper- 
cube. If there are adjacent vertices on the 
hypercube for variables, these must ap- 
pear either as adjacent vertices on the 
hypercube for the four variables, or on the 
hypercube for the remaining variables, but 
not on both of these hypercubes. A cell 
which is not a basic cell of the four-dimen- 
sional hypercube is not used in forming 
basic cells of the n-dimensional hypercube. 
Cells which are basic cells for the four- 
dimensional hypercube may separately or 
collectively form portions of basic cells for 
the n-dimensional hypercube.2 When basic 
cells of the four-dimensional hypercube can 
be combined in basic cells of the -dimen- 
sional hypercube, such combinations appear 
as redundancies on the Veitch-Karnaugh 
map in the second step of the reduction 
process. This argument gives the essential 
basis for the reduction scheme. 

In brief, a higher order switching func- 
tion may be reduced by working succes- 
sively with variables in groups of four. The 
reason for choosing four, rather than any 
other equally valid integer, is that sets of 
four variables are easily reduced by means 
of a Veitch-Karnaugh map. The reason for 
applying the map method of reduction is 
that it places in evidence the role of re- 
dundancies (don’t-care conditions). 

The method used suggested is now ap- 
plied, as an example, to the determination 


f (Bi; a2, *- 


* Received by the PGEC, February 27, 1958. 

1 R. H. Urbano and R. K. Mueller, “A topological 
method for the determination of the minimal form of a 
Boolean function,” IRE Trans. oN ELECTRONIC Com- 
PUTERS, vol. EC-5, pp. 126-132; September, 1956. 

B. Harris, “An algorithm for determining minimal 
representations of a logic function,” IRE TRANS. ON 
ELECTRONIC COMPUTERS, vol. EC-6, pp. 103-108; 
June, 1957, 

_ » There are several sets of basic cells from the four- 
dimensional hypercube, and cells from different sets 
may be combined, as has been done in obtaining (5) 
from the map in Fig. 2. 
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Correspondence 
TABLE I 
Carry Dicit For BINARY CoDED-DECIMAL ADDER 
A 
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 iS 
10) 0 0 0 0 0 0 0 0 0 0 d d d da d d 
1 0 0 0 0 0 0 0 0 0 1 d ad d d d d 
2 0 0 0 0 10) 0 0 10) 1 1 d d d d da d 
3 0 0 0 0 (0) 0 0 1 1 1 d d d d d d 
4 0 0 0 0 0 0 1 1 1 1 d d d d d d 
5 0 0 0 0 0 1 1 1 1 i d d d d d d 
B 6 0 0 0 0 1 1 1 i 1 1 d d d d d d 
7 0 0 0 1 1 1 1 1 1 1 d d d d d d 
8 0 0 1 1 1 1 1 1 1 1 d d d da ad d 
9 0 1 1 1 1 1 1 1 1 1 d d d d d d 
AS d d d d d d d d d d d d d d d d 
12 
13 all d 
14 
15 
Ao 
Z s 
0 0 0 d 0 0 10) d 0 10) 1 d 10) 0 1 d 
0 0 0 da 0 0 1 d 0 0 1 d 0 0 1 d 
0 0 d d 0 0 d d 0 0 d d 0 0 d d 
0 0 d d 0 0 d d 0 1 d d (0) 0 d d 
0 0 1 d Ose 0 is d 0 1 1 d (0) 1 1 d aS 
0 0 1 d 0 1 1 d 0 1 1 d 0 il 1 d 
0 1 d d 0 1 d d 0 1 d d 0 1 d d 
0 1 d d 0 1 d d 1 1 d d 0 1 d d AG 
ya 
d d d d 
As We 
0 1 ag d 0 1 1 d 
0 1 1 d 1 1 1 d 
1 1 d d 1 1 d d d d 
Neel 1 d d 1 1 d d 
WW Yo 
Ai 
Fig. 1—First reduction with A variables. 
Be 
Z Ss 
AsAi 8 Ag 
AsAoAi 
IS 
As Ags Ag Ag 
AsAe AsA2 Ag Ag 
AsA1 AiA2 Bs 
vA 
D D D D 
Bs ve 
As As 
Ag A4 
Ao Az D D 
SS Aj 
S VA 
Bi 


Fig. 2—Second reduction with B variables. 


of a switching function for the carry digit in 
a binary coded decimal adder with inputs As, 
Ax, Ao, Ai, Bs, Bs, Bz, Bi. Table I shows the 
addition table for producing a carry. 
The first reduction will be performed on 
the set of A variables, with the correspond- 
‘ing binary coefficients of the complete prod- 
ucts of the B variables shown in groups of 
four on the map in Fig. 1. The reduction is 
carried out by dealing first with the initial 
element in all boxes of the map, then with 
the second element, and so on. 
The resulting switching function is 


C=1As41+A s+03(Ast+AsA2Ai) 
+04(A s+A4A2) +05(As+Asdo+AsA1) 
+16(As+A4)+07(As+As+A1A2) 
+13(As+As+A2) +09(As+Ast+A2+A1) (4) 


where the v; are the complete products of the 
B variables. 

Next, reduction is carried out with re- 
spect to the B variables as shown in the 
map in Fig. 2. The resulting switching func- 
tion for the carry digit is 


C=A5(Bs+B2+Bi)+As(BiB2+Bs) 
+A2(Bs)+A1(BsBi) +A sAi(Bi) 
+ A4A2(Bs) +A241(BsBoBs) +AsAi(B,B:) 
+A4AoA1(BoBi). (5) 


In obtaining (4) and (5) the redundancies 
have been used. 


Joun N. WARFIELD 
Purdue University 
Lafayette, Ind. 
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Logical Design Using the Stroke 
Function* 


INTRODUCTION 


Recently there has appeared an article! 
on a transistor circuit similar to the type 
shown in Fig. 1. The operation of the cir- 
cuit is such that if any of the inputs are 
positive the transistor is turned off and the 
output is negative. If all of the inputs are 
negative the transistor is turned on and the 
output is positive. Positive and negative in 
this instance refer to potential levels of 
ground and —4&;, respectively. A “polarity 
truth table” for such a circuit having two 
inputs is 


OUTPUT 
B 


INPUTS 


Fig. 1—Transistor gate circuit. 


Using a positive logic (+ is identified with 
true and ONE) the circuit represents the 
binary NOR function, 41+A2=B. A nega- 
tive logic (— is identified with true and 
ONE) results in the binary NAND func- 


tion, A 1 * A 2> Vase 


SHEFFER STROKE 


The Sheffer stroke is a function which 
has not been applied to any great extent to 
the design of logical circuits. The function 
is stated as: A1| Ag is true if either Ai or A>» 
or both are false. An examination of the 
“polarity truth table” shows that the circuit 
in Fig. 1 can be represented by this func- 
tion. Using the negative logic convention 
one has 


0= + = fal AT OniOrt 
nie eed Bian Ae 
1 = — = true 42 OOS 
iy al © 
The positive logic can also be used, giving 
0=-—=false 411010 
B= A,| Ao. 
1 = + = true Az1100 
1B OW O at 


* Received by the PGEC, March i4, 1958. This 
work was supported by the Office of Naval Research 
under Contract No. N7onr41906 and the Dept. of the 
Navy under Contract No. N600(167)45550. 

1 W. D. Rowe, “The transistor NOR circuit,” 1957 
IRE WESCON ConvENTION RECORD, pt. 4, pp. 231- 
245. 
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Cc 


Fig. 2—Half adder based on the stroke function using 
circuits of the type shown in Fig. 1. 


For convenience we choose the negative logic 
representation. 

The relations between the stroke func- 
tion and the more familiar Boolean func- 
tions are 


A|A= 4, (1) 
A;| A» = Ay:Ao, (2) 
A,| Az = Ai + An. (3) 


Two other relations, well suited for computer 
circuits, are 


(Ai| Az) | (Ai | A2) = (41 = 42), (4) 
(A; | As) | (A1| 42) = (Ar = 42). (5) 


Using the stroke functions in (1)-(3) one 
can design circuits consisting only of units 
of the type shown in Fig. 1, which will per- 
form all of the logical operations necessary 
in computers. Relations (4) and (5) are in- 
troduced since they appear in registers em- 
ploying full and half adders. 


DESIGN OF CircuITs USING THE 
STROKE FUNCTION 


Using half and full adders as examples of 
the design of logical circuits involving the 
stroke notation and blocks consisting of the 
circuit shown in Fig. 1, one proceeds as fol- 
lows. 

The truth table for a half adder is shown 
below, where X and Y are the binary digits 
to be added, and Z and C are the sum and 
carry, respectively. 


It is readily apparent that the sum, Z, can 
be formed by 


Z= (x Y)\ (XY); 


and the carry is just 


C= XY. 
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Fig. 3—Full adder based on the modified stroke function using circuits of the type shown in Fig. 1. 


The half adder shown in Fig. 2 is con- 
structed of blocks, each of which represents 
a circuit of the type shown in Fig. 1. The 
function of each block is inscribed in the 
block. 

A full adder can be constructed in the 
manner of using two half adders and an 
OR gate but an economy of circuitry can 
be effected by working directly with the re- 
sults of a truth table. 


26, Wil Oil Mla 
War Oa) ih ake Oy. @ a al 
Cag, OO OW i alata 
Zi Oe ORe I OROR 
Ce OO ie Waitt al 


In the above table X; and Y; are the 7th 
digits to be added, C;_1 is the carry from the 
previous stage, Z; is the sum, and C; is the 
carry. No great effort is required to show 
that 


Z,= (Gia= &% = ¥,)| 


|X; = (Ex = 


Vi) | 
IY; = (Ca = Xi)] 


and 
C; = [C4| (Xs| ¥s)] | [C.| (X;| ¥.)]. 


The design of the full adder can now pro- 
ceed on the basis of forming the Z; and C; 
functions from (1) through (5) and then 
following a functional design of the type ex- 
hibited in Fig. 2. 


MOopDIFIED STROKE FUNCTION 


Restricting the circuit in Fig. 1 to only 
two inputs is a waste of the circuit’s capabil- 
ity. True, there is a practical upper limit to 
the number of inputs determined by the 
finite back resistance of the diodes but this 
limit is certainly greater than two. To de- 
scribe the action of a circuit having ” inputs 
we introduce a modified stroke function and 
define it by stating that: A1|A2| +--+ |Anis 


true if A;, Ao,---, Ano, or An is false 
(using the inclusive OR). Thus in negative 
logic a truth table for three variables would 
be 


A,0101 0101 
4.0011 0011 B=A,|A2| As. 
Az O00: 0 stated 
Foe to en a a a 


Returning to the example of the full 
adder we see that we desire to construct Z; 
and C; from combinations of X;, Y;, and 
C;-1, taking all three at a time. Examina- 
tion of the truth table above shows that 
A,|A2|A3 is £0” only when Ay, Aa, and As 
are all “1’s,” therefore it is possible to obtain 
a “O” at any position by use of the NOT 
functions of A1, As, and A3. A list of these 
relations is given below. 


Ai | As Ay 1 el od 
Ail ds|s4e 1 od 
Aylas| Acted: Oda 
AylAy| AY Ya V1 Ot 
Ay Malta led t, Oleieleiet 
Ay Ay | dents OAc el ed 
Ail As Ay Oe et 
Ags) AgileAigee Ou toilet ede apelin! 


Referring to the above list and the truth 
table for the full adder it is readily apparent 
that 


Zi = (X;| Yi] Cu) | (X, | P;| C4) | 
-(X;| Y;| Cea) (x; | Y; | Cat 
and 
Cy = el Vel Cea) OR Paes 
“(Xs | ¥s| Cua) | | Ye] Cu). 
A block diagram of a full adder based on 
circuits of Fig. 1 and the above relations for 
Z; and C; is shown in Fig, 3. 
Various methods of notation can be used 
to designate complex stroke functions other 


than the usual set of brackets, parentheses,. 
etc. In fact, it is possible to designate alb 
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possible functions using only the stroke no- 
tation since the NOT function has an 
equivalent stroke notation. This is done by 
replacing the brackets with multiple strokes 
and the NOT function with its equivalent 
stroke function. As an example, consider the 
‘sum and carry functions mentioned above. 
They could be written 


Z;=X;|| Vil]. ||| X<| Xi] ¥e| PellC I] Xe 
» Xil| Vela | Cea lla] Pe] Vil|Cua| Cr, 

Ci=Xil| Val|Ci ||] Xs Xe Pell Ce [||] ¥6 | 
- Yil|Ca||| Xi] ¥il|Ca] Ga. 


Obviously this is a matter of one’s own taste 
and the authors have found that their 
preference is to use brackets, the NOT func- 
tion, and to omit the stroke in the lowest 
level of the logic so that the sum and carry 
would appear as 


Contributors 


D. D. Aufenkamp, for a biography see 
page 205 of the September, 1957, issue of 
these TRANSACTIONS. 
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Walter F. Bauer (A’53) was born in 
Wyandotte, Mich., on March 21, 1924. He 
received the B.S., M.S., and Ph.D. degrees 
from the University of Michigan, Ann Ar- 
bor. He was an Air Force weather officer 
during World War II, and research engineer 
at the Willow Run Research Center, Uni- 
versity of Michigan from 1951 to 1954. 

He has been with Ramo-Wooldridge Cor- 
poration since 1954, and currently directs 
the Computation and Data Reduction Cen- 
ter, Space Technology Laboratories. His main 
interests are system simulation, computer 
systems, numerical analysis, data reduction, 
and programming. He has been directly in- 
volved with missile simulation on digital 
computers since 1951 with the SEAC com- 
puter at NBS and with the MIDAC com- 
puter at the University of Michigan. He has 
written several papers on computer systems 
and simulation. 

Dr. Bauer is a member of the Associa- 
tion for Computing Machinery and the 
American Mathematical Society. 
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J. V. Blankenbaker (A’54) was born on 
December 24, 1929, in Apache, Okla. Fol- 
lowing two years as an electronic technician 
in the U. S. Navy, he attended Oregon State 
College, Corvallis, Ore., where he received 
the B.S. degree in 1952 in physics and ma- 
thematics. In 1954, he received the M.S. 
degree in physics from the University of 
California at Los Angeles. 

From 1952 to 1956, he was engaged in 
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Z=XiViCi1 | XiViCia | XV Ces | XViCe, 
C= XV iC | XiViCia| XiViCia| XiViCu. 


Higher order functions would of course re- 
quire the use of brackets, parentheses, etc. 
No loss in meaning or understanding results 
from this “simplified” notation but the 
complex notation using multiple strokes does 
provide a convenient means for determining 
the number of transistors in the circuit. For 
example, in the Z; relation there are three 
different single stroke relations, four differ- 
ent double stroke relations, and one triple 
stroke relation making a total of eight 
transistors required to generate Z;. Only 
twelve transistors are required to generate 
both Z; and C; because four of the inter- 
mediate functions are common to both. This 
also is evident when the two functions are 
compared with each other. 


the design of data processors and special 
purpose computers at Hughes Aircraft 
Company. During the summer of 1957, he 
did computer design studies at the Ramo- 
Wooldridge Corporation, Los Angeles, Calif. 
Since 1956, he has been working toward the 
Electrical Engineer degree at the Massa- 
chusetts Institute of Technology. 

Mr. Blankenbaker is a member of Phi 
Kappa Phi and the ACM. 


\7 
“eo 


Edward L. Braun (SM’55) was born on 
February 18, 1925 in Pittsburgh, Pa. He re- 
ceived the B.A. degree in physics in 1947 and 
the M.S. degree in engineering in 1954 from 
the University of California at Los Angeles. 
While in the U. S. Navy, he received exten- 
sive education in electronics, principally in 
radio and radar. 

His professional experience and first 
interest in digital computers began at 
Northrop Aircraft, Inc., where he stayed 
from 1947 to 1952. From 1952 to 1955 he was 
with Librascope, Inc., and the following year 
he was with Lockheed’s Missile Systems Di- 
vision. He was subsequently employed by 
Hughes Aircraft Company and Litton In- 
dustries, Beverly Hills, Calif. Now he is 
director of engineering of Genesys Corpora- 
tion (a subsidiary of Chance-Vought Air- 
craft Company, Inc.), Los Angeles, Calif. 

Mr. Braun has organized and headed a 
number of digital computer development 
programs. He has designed digital computers 
for military computation and control sys- 
tems, and for engineering and scientific com- 
putation. While at the Lockheed Missile 
Systems Division, he also served for a time 
as head of the Digital Computing Center. 

Mr. Braun is a member of the ACM, and 
served a term as membership chairman of 
the Southern California section of the ACM. 
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CONCLUSIONS 


A symbolic notation has been developed 
which fits the transistor gate circuit shown 
in Fig. 1 to the Sheffer stroke function. An 
extension of the stroke function has been 
developed to aid in the application of the 
logic to complex circuits. The two examples 
shown are relatively simple circuits but the 
extension to more complex circuits is a 
simple matter. It is felt that a considerable 
simplification in the design of computer cir- 
cuitry can be achieved if the elemental 
circuit blocks are all of the same form, 
allowing one to express all logical operations 
in the form of the stroke function. 

N. T. GRISAMORE 

L. S. RoTOLO 

G. U. UYEHARA 

George Washington Univ. 
Washington, D. C. 


William Comley was born in Indianapo 
lis, Ind., on September 14, 1921. After com 
pleting the electronics curriculum at Los 
Angeles City College, Calif., he served as a 
communications officer with the U. S. Air 
Force for three years and was for the greater 
part of this time in charge of a number 
of high-powered ground transmitting sta- 
tions. 

In 1947, Mr. Comley joined the engineer- 
ing staff of Lear, Inc. of California as a de- 
sign engineer in aircraft radio. He later 
served as project engineer for that com- 
pany’s line of phonograph pickups and wire 
recorder heads. 

In 1950, Mr. Comley joined Douglas 
Aircraft Company, where he is currently en- 
gaged in the design of analog computing de- 
vices and auxiliary equipment for existing 
facilities. 
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Robert G. Gillespie was born May 18, 
1935 in Tacoma, Wash. He received the 
B.A. degree in mathematics from Reed Col- 
lege, Portland, Ore., in 1955. 

From 1955 to 1956, he was at Convair 
Aircraft Company, San Diego, Calif., in the 
Digital Computing Laboratory. From 1957 
to 1958, he was at the Astronautics Division 
of Convair, where he worked on the digital 
simulation of missile guidance systems. He 
now is a member of the technical staff of the 
Pilotless Aircraft Division, Boeing Airplane 
Company, Seattle, Wash., where he is en- 
gaged in work on systems evaluation and 
simulation. 

He is a member of the American Mathe- 
matical Society, the Association for Com- 
puting Machinery, and the Society for In- 
dustrial and Applied Mathematics. 
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William F. Gunning (S’41—-A’43-—M’55) 
was born on August 11, 1916, in New York, 
N. Y. He received the A.B. degree in physics 
from the University of California at Los 
Angeles in 1941. His professional experience 
started at Douglas Aircraft Company’s 
Santa Monica Electrical Research Labora- 
tory, where he worked on flight and static 
test instrumentation from 1941 to 1948. He 
transferred to Project Rand and stayed 
with this organization when it became the 
Rand Corporation in 1948. He was in charge 
of rebuilding an early model REAC analog 
computer, adding chopper stabilizer ampli- 
fiers, and the first centralized removable 
program plugboard and operator’s console 
to be used on a large analog computer in- 
stallation. After a tour on part-time loan to 
the group constructing the SWAC at 
U.C.L.A., he became project engineer on 
Rand’s Johnniac digital computer which is 
a member of the Princeton I.A.S. family of 
machines. 

From 1953 to 1956, he was with Inter- 
national Telemeter Corporation working on 
problems associated with pay-as-you-see 
TV and magnetic core storage systems. 

He josned the Systems Division of Beck- 
man Instruments in 1956 where he is Chief 
Project Engineer working on data logging 
and data acquisition equipment for indus- 
trial application. 

Mr. Gunning has served on several Joint 
Computer Conference committees. 
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Ladis D. Kovach (A’53) was born on 
November 21, 1914 in Budapest, Hungary. 
He received the B.S. degree in physics in 
1936 and the M.S. degree in mathematics in 
1948 from Case Institute of Technology, 
Cleveland, Ohio; the M.A. degree in educa- 
tion from Western Reserve University, 
Cleveland, Ohio, in 1940; and the Ph.D. 
degree in mathematics from Purdue Uni- 
versity, Lafayette, Ind., in 1951. 

From 1936 to 1948, he was associated 
with Picker X-Ray Corporation, American 
Shipbuilding Company, and Ohio Crank- 
shaft Company as electrical designer. He was 
an instructor in mathematics from 1946 to 
1951, first at Case and then at Purdue. Since 
1951, he has been with the Douglas Aircraft 
Company as a design specialist. 

Dr. Kovach is a member of the American 
Physical Society, Sigma Xi, and the Ameri- 
can Mathematical Society. 
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Morton W. Marcovitz (S’50-—A’52—M’57) 
was born on September 20, 1927, in Phila- 
delphia, Pa. He received the B.S.E.E. de- 
gree from Drexel Institute of Technology in 
1951, and the M.S.E.E. degree from the 
University of Pennsylvania in 1957. 

He has been associated with the re- 
search activity of the Burroughs Corpora- 
tion since 1951, working on applications of 
digital computer techniques to data process- 
ing systems. 

Mr. Marcovitz is a member of Eta 
Kappa Nu, Tau Beta Pi, AIEE, ACM, and 
RESA. 
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Douglas B. Netherwood (M’56) was born 
in Dallas, Texas, on May 8, 1920. He re- 
ceived the B.S. degree from the U. S. Mili- 
tary Academy in 1943 and the M.S. degree 
in electronics from the Air Force Institute of 
Technology in 1957. He served in India and 
Burma during World War II and was Ad- 
visor on Research and Development to the 
Nationalist Government of China from 1946 
to 1947. He was a Communications Officer 
in the Strategic Air Command from 1947 to 
1955. At present he is Chief of the Electron 
Tube and Circuit Research Section, Elec- 
tronic Components Laboratory, Wright Air 
Development Center, Ohio. His field of 
special interest is the simplification of elec- 
tronic components and networks. 

Major Netherwood is a member of the 
ACM and Association for Symbolic Logic. 
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Palakodety Venkata Rao was born Au- 
gust 18, 1922. He received the Master’s de- 
gree in physics from Andhra University, 
Waltair, India, in 1943, and then studied 
electrical engineering at the Indian Institute 
of Science. He later worked as a design elec- 
trical engineer with the India Electrical 
Works in Calcutta, and since 1948, has been 
a member of the faculty of the Indian Insti- 
tute of Science at Bangalore. He will shortly 
complete his work for the D.Sc. degree in 
electrical engineering at Andhra University. 

Mr. Rao is an associate member of the 
AIEE. 
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Joseph H. Rawlings was born in Omaha, 
Neb. on April 16, 1924. He received the B.S. 
degree from Western Reserve University, 
Cleveland, O., in 1952. 
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He has been with the National Advisory 
Committee for Aeronautics since 1946 and is 
now head of Computer Engineering at the ~ 
Lewis Flight Propulsion Laboratory, Cleve- 
land. 
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Eric Seif (S’52—A’53) was born on Sep- 
tember 25, 1925, in Vienna, Austria. He re- 
ceived the B.S. degree in electrical engineer- 
ing from Drexel Institute of Technology, 
Philadelphia, Pa., in 1952 and the M.S. de- 
gree from the University of Pennsylvania, 
Philadelphia, Pa., in 1957. 

Since 1952 he has been with Burroughs 
Corporation Research Center, Paoli, Pa. 

Mr. Seif is a member of Tau Beta Pi, 
Eta Kappa Nu, Phi Kappa Phi, and RESA. 


‘7 
“e 


Alfred K. Susskind (S’47—A’50-M’56) 
was born on October 4, 1923, in Coblenz, 
Germany. He received a B.E.E. degree from 
the Polytechnic Institute of Brooklyn in 
1948, and the S.M. degree from M.I.T., 
Cambridge, Mass., in 1950. 

From 1948 to 1950, he was a research 
assistant in the Digital Computer Labo- 
ratory at M.I.T., working on the electronic 
design for the Whirlwind Computer. From 
1950 to 1955, he was a research engineer in 
the Servomechanisms Laboratory, working 
on numerical control of machine-tools. 

Since 1955, he has been an assistant pro- 
fessor in the Department of Electrical En- 
gineering at M.I.T., giving courses in digital 
data processing engineering, and doing re- 
search in airborne data-processing problems. 

Mr. Susskind is a member of Eta Kappa 
Nu, Tau Beta Pi, and Sigma Xi. 
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L. R. Turner was born on August 19, 
1910, in Williamsport, Pa. He received the 
B.S. degree in physics in 1931 and the M.S. 
degree in 1932 from Pennsylvania State Uni- 
versity, University Park. 

He has been with the National Advisory 
Committee for Aeronautics since 1938, work- 
ing in engineering until 1949 and then in 
applied mathematics. He is Chief of the 
Mechanized Computation and Analysis 
Branch at Lewis Flight Propulsion Lab- 
oratory, Cleveland, Ohio. 

He is a member of the ACM, Pi Mu 
Epsilon, and Sigma Rho Sigma. 
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JOINT COMPUTER COMMITTEE 


SENEWS 


SCIENCE EDUCATION SUBCOMMITTEE NEWSLETTER 


Vols 1,-No..1 


Wuat Is SENEWS? 


A newsletter addressed to computer-oriented mem- 
bers of IRE, ACM, and AIEE to help them promote 
computer interest and knowledge among high-school 
age students. 


Wuo Is RESPONSIBLE? 


The Science Education Subcommittee under sponsor- 
ship of the Joint Computer Committee. C. W. Farr of 
Lincoln Laboratory, M.I.T., is Chairman of the Sub- 
committee and Acting Editor of SENEWS. 


Way Is Ir? 


Computers are here to stay. High-school students are 
already showing a budding curiosity about computers. 
Professional computer people are showing wide interest 
in helping the kids to get started. All over the country 
computer educational activity is springing up. A rally- 
ing point is needed for sharing ideas and results. 


SEND Us A COMPUTER STORY 


... of your local high-school activity. If there is no 
news story, guess you'll have to go out and make some 
news. Let’s challenge the youngsters to reach their po- 
tential strength in the data processing world of to- 
morrow. Watch SENEW’S for ideas, and send us your 
success stories! 

Write to: C. W. Farr, JCC Science Education Sub- 
committee, M.I.T. Lincoln Laboratory, Lexington 73, 
Mass. 


Have DeEsicn, WILL BUILD 


A high-school junior and perennial regional science 
fair winner from Cedar Falls, Iowa, named David 
Ecklein, has designed and is building a thirty-five 
hundred vacuum tube digital computer, using surplus 
market tubes and sockets and adapting office shelving 
for use as chassis and cabinet material. He has designed 
the mathematical logic and circuitry with essentially 
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no assistance other than normal instruction in high- 
school science, what he has been able to learn from 
books, and his own experiments. Part of the circuitry 
has been assembled, with the help of a high-school soph- 
omore friend whom he, in turn, has helped on another 
science project. 

The point of the story is that the ability and motiva- 
tion for a budding computer engineer can spring up in 
the middle of Iowa or in the suburbs of Boston or Los 
Angeles, for example. There is one big difference: in the 
coastal industrialized areas, computers and computer 
engineers abound, and a youngster finds it is easy to 
obtain local adult professional help in the pursuit of his 
project. 

The JCC Science Education Subcommittee heard of 
this boy by way of Washington, D. C. (the President’s 
Committee on Scientists and Engineers), and is proud 
to have brought him to the attention of JCC members 
who have set the wheels in motion to give him well de- 
served counsel and encouragement. (Ed. note: David 
has sold over 500 Bibles to raise money for his science 
project. At last report he was flying to New York ona 
Saturday—Tuesday visit as guest of IBM.) 


Low-BupGET HIGH-SCHOOL COMPUTER 


In Buffalo, N. Y., Aaron Buchman, mathematics in- 
structor at the Hutchinson Central Technical High 
School, has been conducting a course in computer pro- 
gramming and coding open to college preparatory sen- 
iors. Under his direction, members of the Tech Math 
Club constructed a tape activated automatic digital 
computer for use in the course. 

The one-semester course consists of theory and lab- 
oratory work. After 30 classes in theory, the students 
get two laboratory and three theory periods per week. 
Included in the theory are number systems, drill in 
binary system arithmetic operations, elements of Boo- 
lean algebra, and computer logic. 

The computer is an electromagnetic relay machine 
which operates in the binary number system. Tapes are 
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ON ELECTRONIC 


COMPUTERS June 


A. Vanderburgh, R. Enticknap, N. Levy, and R. Mayer (back row, left to right) conduct a visit of Concord 


High School students to M.I.T. Lincoln Laboratory’s experimental computer, TX-2. 


punched out of pressboard nine inches wide, and carry 
ten channels. Each of the six registers of the storage 
section holds a six-digit binary number. Panel lights in- 
dicate the contents of each register at all times. With 
two operator registers and one output register, the 
arithmetic unit performs direct addition, subtraction, 
and multiplication. Division is accomplished by a taped 
routine. 

Before completing the course, each pupil must pro- 
gram and code an approved program which requires 
branch points, cycling routines, a minimum number of 
orders, and which must run without trouble on the com- 
puter. 

Among the problems successfully programmed and 
now a part of the library of routines are those to: 


1) Find the total surface area of a rectangular solid’ 

2) Find the average of two signed numbers, 

3) Evaluate a second degree polynomial using suc- 
cessive integer values of the variable, 

4) Solve simultaneous linear equations, 

5) Evaluate a two-rowed determinant, 

6) Evaluate various formulas from mathematics and 
physics for sets of values of the variables involved 
by having the computer order the machine to take 
a new value of one of the parameters with each 
pass of the tape. 


A description of the machine and the mathematics 
course in which it is used was presented by Mr. Buch- 
man at the ACM meeting in August, 1956. Mathematics 
Teacher carried the story in March, 1958. The course 
now includes three weeks of IBM 650 programming. 

Here, a teacher’s initiative and a forward-looking 
school system have combined to challenge the abilities 
of bright high-school students with a solid introduction 
to the technological revolution taking place within their 
time. 


STUDENTS STUDY 


Concord (Mass.) High School’s Mathematics Depart- 
ment experimented in the 1958 spring term with a six- 
week curriculum unit on computers. Fifteen students 
participated in the unit, which was designed by Norton 
Levy of the Concord High School Mathematics De- 
partment, assisted by M.I.T. Lincoln Laboratory staff 
members, who acted as consultants, arranged computer 
facilities tours, and furnished a demonstration-lecture 
team for a major part of the course. Lectures included 
number systems, computer principles (slanted away 
from electronic circuitry), simple programming theory, 
study of actual coding of two programs for the IBM 704 
computer, a follow-up visit to a computer facility while 
the programs were being run, and computer applica- 
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tions. Course evaluation was by written examination. 
This curricular experiment, guided by a teacher-engi- 
neer team, was sponsored by Ford Foundation’s School 
and University Program for Research and Develop- 
ment, under grant to Harvard’s School of Education. 

An interesting sidelight of this story is the “card- 
board computer,” a mechanical model made from card- 
board and common pins by R. P. Mayer in one evening 
and used in a high-school lecture-demonstration. This 
is an outgrowth of the effort described by Mayer in his 
paper, “A Proposal for Training Youngsters in Digital 
Computing Techniques,” presented to the ACM in 
Los Angeles in August, 1956. The “cardboard computer” 
was a binary unit 83 by 11 inches, comprising a man- 
ually activated pulse input capable of producing either 
a “zero” or a “one” pulse, a counter, control element se- 
lecting between two one-bit storage registers, and a one- 
bit adder. 


So Do TEACHERS 


A high-school teacher computer course consisting of 
seven lecture and workshop meetings is described in 
“Tntroduction to Arithmetic Computers,” a report to 
the Hughes Aircraft Company by ten participating 
(high-school science and mathematics) teachers, Los 
Angeles School-Industry Science Program, December, 
1956. (Copies are available from Dr. L. C. Van Atta, 
Head, Technical Information and Education, Hughes 
Aircraft Company, Culver City, Calif.) The report is in- 
tended for use as a syllabus for a teachers’ 40-hour 
workship. The lectures reviewed computer history, 
binary arithmetic, elementary computer logic, and in- 
cluded visits to a mechanical differential analyzer, to 
the Hughes REAC (analog) and IBM 650 (digital), 
and to the Bureau of Standards SWAC (digital) com- 
puter facilities. One lecture was devoted to SWAC pro- 
gramming; homework called for completion of the 
coding of a fifth degree polynomial problem (appendix 
to the report contains detailed flow diagram and coded 
program). 

This computer activity is only one phase of a well- 
organized community program! in which the school 
system employs a full-time school-industry coordinator 
who works with a Council. The Council’s board of di- 
rectors includes representatives of secondary schools, 
universities, industry, and the professional societies; 
the Council publishes a newsletter, organizes industrial 
demonstration-lecture teams for high schools, coordi- 
nates an industrial summer work program for high- 
school teachers, and in July, 1957 participated in an 
industry-education conference financed by Hughes Air- 
craft Company and held at the Lake Arrowhead facili- 
ties of the University of California. For information 


1L. C. Van Atta, “The Los Angeles school-industry science pro- 
gram,” Phys. Today, vol. 10; September, 1957. 
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write to J. H. Cooper, Executive Secretary, Southern 
California Industry-Education Council, 954 West 37th 
Street, Los Angeles 7, Calif. 


DID-Y OU Ha ART eee. 


... the one about the computer programmer who 
couldn’t balance his checking account. He knew that 
177+1=200; but the bank took a narrower view. 


IRE, ACM, AIEE 


Richard W. Melville (Computer Development Group, 
Stanford Research Institute), is Chairman of the PGEC 
Student Relations Committee, and IRE representative 
of the SENEWS Editorial Board. The original Chair- 
man of the JCC Science Education Subcommittee, he 
was appointed when its charter was a nebula of swirling 
hot gases; he explored and charted the unknown and 
has been a valued consultant in establishing the present 
subcommittee policy and program and in bringing forth 
this newsletter. 

Professor George E. Forsythe (Mathematics De- 
partment, Stanford University) is Chairman of the 
ACM Secondary School Educational Committee (ac- 
tually ACM has two other education committees, Uni- 
versity-Education and Industry-Education). Professor 
Forsythe is ACM representative on the SENEWS 
Editorial Board and has furnished valued criticism and 
encouragement. Those interested should write him for 
his Bibliography on High-School Mathematics Educa- 
tion, which includes ten books on digital computers. 
Professor Forsythe has earnestly questioned our edito- 
rial policy of confining SENEWS to computers, a ques- 
tion that has probably arisen also in the minds of many 
of our readers. 

While it is true that computer engineers as individuals 
can be broadly helpful in high-school science educa- 
tion, the JCC Science Education Subcommittee finds 
many other groups at work in the field and has focused 
on the task of providing “a newsletter addressed to 
computer oriented members of IRE, ACM, and AIEE, 
to help them promote computer interest and knowledge 
among high-school age students.” 

Dr. Morris Rubinoff (Government and Industrial 
Division, Philco) is Chairman of the AIEE Computing 
Devices Committee, and has given the JCC Science 
Education Subcommittee the benefit of his views on our 
newsletter program. In the temporary overseas absence 
of Dr. Rubinoff we are grateful to Vice-Chairman R. A. 
Imm and others in AIEE for counsel and criticism in 
launching SENEWS. 


Book AND FILM 


“Where can I find publications on computers under- 
standable to beginners?” You have probably been asked 
this question. Here is a short-short list of suggestions: 

“Faster Than Thought,” edited by Bowden, pub- 
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lished in 1953 by Pitman, London (available in U. S. at 
$9.50 and worth it). Subtitled “A Symposium on Digital 
Computing Machines,” it covers history and theory, 
surveys British and American computers, and reviews 
computer applications. The book is thorough and is 
written with “sparkle” (e.g., “...an expert is some- 
times defined as one who solves minor problems and 
avoids small errors as he sweeps forward to the ‘Grand 
Fallacy.) 

“The Electronic Brain and What It Can Do,” by 
Gorn, Manheimer, and Brandwein, published in 1956 by 
Science Research Associates (paperback, 50 cents). 
We have shown this to a range of people from high- 
school students and teachers to physicists, and they all 
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liked it. It uses several card tricks to Jemonstrate yes-no 
logic. 

“Making Electrons Count” is a 16-mm sound film 
telling the story of programming for the Whirlwind I 
digital computer; old (1954), but good for beginning 
audiences. This half-hour color film is available for 
mailing charges from M.I.T. Computation Center. 


FEEDBACK 


If you have read this far, you are identified as a person 
interested in our program. How about a story of high- 
school computer activity from your locality, or a word 
of constructive criticism, or just tell us you liked our 
newsletter, or just tell us! 


PGEC News 


PATENTS RELATING TO COMPUTERS 


Is there enough interest among TRANSAC- 
TIONS readers to warrant publishing some- 
thing in each issue about computer patents? 
The simplest and easiest thing to do is to 
publish a list giving number, date, in- 
ventor’s name, assignee, and title. Many of 
these would be from classification 235-61. 
More elaborate information might include 
an abstract of the disclosure. Readers in- 
terested should let the editor know of their 
interest. They should also realize that some 
help from one or more of them will have to 
be solicited if a decision is made to proceed 
with the publishing of patent information. 


New AFFILIATE 


The Institution of Electrical Engineers 
of Great Britain has been approved as an 
affiliate of the PGEC. 


Dr. WILLIAM G. FULLER MEMORIAL 
AWARD 


The Professional Group on Component 
Parts has announced that several awards of 
$250 will be made to seniors and graduate 
students for outstanding papers on the sub- 
ject of component parts submitted by De- 
cember 31, 1958, to any IRE publication or 
presented at any national meeting or sym- 
posium sponsored or supported by the IRE. 
Interested students are urged to consult 
their IRE Faculty Advisor. 


SUMMARY OF PGEC ApDMINISTRA- 
TIVE COMMITTEE MEETING OF 
Marcu 26, 1958 


The final meeting of the 1957-1958 Ad- 
ministrative Committee was held at IRE 
Headquarters, New York, N. Y. on March 
26, 1958, during the IRE National Conven- 
tion. A total of 15 members were present or 
represented by proxy. 

New officers were elected for the year 
starting April 1, 1958: 


Chairman—Willis H. Ware, Santa Mon- 
ica, Calif. 

Vice-Chairman—Richard O. Endres, 
Philadelphia, Pa. 


Five new Administrative Committee 
members were elected for a three-year term, 
starting April 1, 1958: 


Walter L. Anderson, Arlington, Va. 

Arnold A. Cohen, St. Paul, Minn. 

Daniel Haagens, Long Island City, N.Y. 

Frank B. Heart, Cambridge, Mass. 

Keith W. Uncapher, Santa Monica, 
Calif. 


The First International Conference on 
Information Processing, as it is officially 
called, is now definitely set under UNESCO 
sponsorship for June 15-20, 1959, according 
to I. L. Auerbach’s report. The place will be 
either Paris or Rome, to be decided this 
summer, and the method of selecting papers 
is being worked out. 

The conference subjects are listed under 
these headings: 


Methods of digital computing, 

Logical design and common symbolic 
language for machines, 

Automatic translation of languages, 

Collection, storage, and retrieval of in- 
formation, 

Pattern recognition and machine learn- 
ing. 


A lively discussion followed the report of 
the Bibliography Committee. Committee 
Chairman L. F. Jones presented several 
ways in which the series of abstracts in the 
TRANSACTIONS, formerly prepared under 
the editorship of H. D. Huskey, could be 
continued, A distinction must be made be- 
tween the regular publication of abstracts 
or reviews of current literature and the occa- 
sional publication of a bibliography of com- 
puter literature up to a certain date, includ- 
ing careful classification and cross-indexing 
of the material. It was concluded that the 


abstracts were the more urgent at this time 
to avoid too long a gap after the termination 
of Huskey’s series. A bibliography, however, 
presents an equally important task, as well 
as a more difficult one, and the Committee 
will work on it concurrently. 

The interest at the meeting indicates 
that the preparation of a good bibliography 
and continuation of the abstract series could 
be a most important service for the PGEC 
to give its members. In recognition of this, 
the Bibliography Committee was elevated 
to the rank of a permanent committee. 

At the conclusion of another successful 
year for the PGEC, I wish to thank mem- 
bers of the Administrative Committee and 
the various committee chairmen for their 
efforts and cooperation. The PGEC is the 
largest Professional Group in the IRE and 
it is still growing. This reflects not only the 
continuing growth and importance of the 
computer field, but also the contributions of 
the active members which have made the 
PGEC an important institution in the tech- 
nical field. 

In particular, I would like to mention 
Ike Auerbach who initiated the organization 
of the First International Conference on 
Information Processing (which somebody, 
following a venerable international custom, 
will surely label ICIP some day). Ike is now 
PGEC representative and chairman of the 
NJCC Committee which organizes the 
United States’ end of the conference. He 
has put in a great deal of work during the 
past year, contacting many organizations 
in this and other countries, and making a 
trip to Paris which culminated in the con- 
ference becoming firmly established. Ike’s 
current unofficial title is Minister of Foreign 
Affairs. 

I hope that this example will inspire 
other PGEC members to seek a field of 
activity in which they can make an effective 
and unique contribution. 


W. BucHHOLz 
Chairman 
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2° INFORMATION FOR AUTHORS 


es 
The PGEC TRANSACTIONS is published quarterly and will bear date- 
lines of March, June, September, and December. Abstracts of papers 
appearing in the TRANSACTIONS will appear also in IRE PROCEEDINGS. 
The PGEC publication schedule requires about one month for review 
and correction of all accepted manuscripts. The professional IRE 
Editorial Staff requires an additional two months’ production time from 
receipt of manuscripts to completion of the printed journal. Y 


MANUSCRIPTS: Three copies of the manuscript should be sub- 
mitted. They should be typewritten (original and two carbon copies), 
and double spaced on only one side of each sheet. References should 
appear as footnotes, numbered consecutively, and include in the fol- 
lowing order the author’s name (including initials), title of reference 
work, journal name, volume, initial and final page numbers, and date 
of publication. Footnotes should be listed on a separate sheet and not 
inserted in text. Each paper must be accompanied by two copies of a 
summary not more than 200 words in length. Reviewing normally will 
require about four weeks from receipt of manuscript. 


ILLUSTRATIONS: Only original illustrations should be submitted: 
they will be returned if desired. Photostatic copies of originals are not 
ac, aptable, except where they are exceptionally clear, with sharp black 
and white contrasts. All line drawings (graphs, charts, diagrams, etc.) 
should be prepared on drafting cloth or white drawing paper in India ink. 
It is preferable that only the coordinate lines show in graphs. All 
lettering must be large enough to be legible when reduced 50 to 75 per 
cent in size. Photographs should be glossy prints. Figure numbers 
should be indicated on the back of each illustration. Figure numbers 
and captions should be listed on a separate sheet accompanying manu- 
script. All drawings, photographs, and other manuscript material should 
_ be not larger than 83 by 11 inches for ease in handling. 


Please submit all manuscripts to 
J. P. Nash, PGEC Editor 
Missile Systems Division, Lockheed Aircraft Corp., 
3251 Hanover Street, Palo Alto, California 


